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Preface 

The  purpose  of  this  study  was  to  characterize  the  performance  of  two  types  of  nonlinear 
optical  crystals  in  their  use  for  producing  the  second  harmonic  of  10.6  (im  radiation.  The 
crystals  and  laboratory  facilities  were  provided  by  the  Wright  Research  and  Development  Cen¬ 
ter’s  Electro-Optics  Sources  Branch,  Wright  Patterson  AFB,  Ohio. 

The  computational  work  in  preparation  for  the  experiments  led  to  the  optimal  designs  for 
characterizing  each  of  the  crystals.  The  results  of  this  analysis  showed  that  AgGaSe2  has  prom¬ 
ise  as  a  high  efficiency  source  of  doubled  10.6  |im  radiation,  while  ZnGeP2  may  be  more 
suited  for  other  applications  such  as  optical  parametric  oscillation.  The  work  should  be 
continued  with  AgGaSc2  in  order  to  more  fully  understand  the  damage  mechanisms  near  the 
points  of  peak  performance. 

In  performing  this  work,  1  have  had  a  great  deal  of  help  from  several  individuals.  I  wish  to 
thank  my  faculty  advisor,  Dr.  Won  B.  Roll,  for  his  guidance  and  assistance  throughout  the 
effort.  1  am  very  grateful  to  my  sponsor,  Dr.  Kenneth  L.  Schepler,  Wright  Research  Develop¬ 
ment  Center  Electro-Optic  Sources  Branch,  for  his  direction  and  enthusiastic  support  as  the 
work  progressed.  To  Dr.  Paul  Schrieber,  I  owe  a  debt  of  gratitude  for  his  help  in  providing 
insightful  comments  and  instructing  me  on  the  use  of  the  instrumentation.  Thanks  go  to  Gail 
Brown  and  Capt.  Ges  Seiger  of  WRDC/MLPJ  for  their  response  in  spectrally  characterizing 
the  crystals  and  other  samples.  A  word  of  thanks  is  owed  to  Robert  Wade  and  Art  Becraft  for 
their  efforts  and  patience  in  fabricating  a  number  of  components,  and  to  Virginia  MacMillan 
for  running  several  FTIR  scans  as  well.  Finally,  1  wish  to  thank  my  family  and  friends  for 
helping  me  to  maintain  a  bit  of  sanity  over  the  course  of  these  pas!  few  months. 

Monte  D.  Turner 
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Abstract 

This  experiment  was  to  characterize  the  nonlinear  performance  of  AgGaSe2  and  ZnGeP2  for 
frequency  doubling  C02  laser  output  at  10.6  pm.  The  experiments  included  a  series  of  theo¬ 
retical  estimates  of  the  source  and  crystal  performance  to  define  the  measurements  to  be  per¬ 
formed.  The  parameters  of  interest  were  the  phase  matching  angles,  angular  acceptance, 
conversion  efficiency  and  damage  thresholds.  The  9mm  x  9mm  x  29mm  antireflection  coted 
sample  of  AgGaSe:,  with  a  phase  matching  angle  of  0  =  55.02°  and  an  internal  acceptance 
bandwidth  A9  =  0.298°,  achieved  a  maximum  second  harmonic  conversion  efficiency  of 
20.7%.  Surface  damage  limited  the  improvement  of  the  efficiency  and  the  threshold  for  the 
coated  surface  was  near  4.5 kWlcm2.  Further  studies  were  conducted  using  small  witness  sam¬ 
ples  from  the  crystal  boule.  For  the  case  of  the  uncoated  sample  of  ZnGeP2,  the  theoretical 
computations  point  out  that  the  performance  will  be  poor  at  best  due  to  an  extremely  small 
effective  nonlinear  coefficient  at  the  temperature  tuned  phase  matching  angles.  These  calcula 
tions  are  presented  though  no  second  harmonic  was  observ  ed  during  these  experiments. 
Improvements  to  the  experimental  designs  and  theoretical  estimates  are  cited  for  continued 
work  with  both  of  these  crystals. 


IX 


CHARACTERIZATION  OF  AgGaSe2  AND  ZnGeP2  FOR 
FREQUENCY  DOUBLING  C02  LASER  OUTPUT 


I.  Introduction 

High  performance  and  frequency  agile  laser  sources  are  required  for  a  variety  of  uses  in 
optical  systems  and  applications.  Non-linear  optical  materials  provide  a  means  to  extend  the 
utility  of  these  sources  by  converting  the  output  frequencies  from  known  laser  systems  to  other 
regions  of  the  optical  spectrum.  The  carbon  dioxide  laser  is  a  particularly  mature  technology 
whose  high  output  efficiency  makes  it  an  attractive  candidate  for  frequency  conversion.  There 
is  a  significant  interest  in  the  mid-infrared  region  near  5  pm  due  to  a  natural  transmission  win¬ 
dow  in  the  earth's  atmosphere  and  applications  in  spectroscopy  and  remote  sensing. 

In  the  development  of  non-linear  optical  materials,  their  utility  in  the  mid  and  far  infrared 
regions  have  been  limited  by  transmission  losses  due  to  absorption  and  scattering.  Select 
materials  from  a  class  known  as  the  chalcopyrites  have  been  used  to  demonstrate  second  har¬ 
monic  generation  across  a  wide  band  of  the  infrared  spectrum  (19:1).  With  a  relatively  low 
intrinsic  absorption  in  this  region,  scattering  losses  resulting  from  the  crystal  growth 
techniques  have  limited  their  ultimate  usefulness.  Recently,  silver  gallium  selenide  and  zinc 
germanium  phosphide,  have  been  grown  in  a  manner  which  significantly  reduces  these  losses 
by  Cleveland  Crystals  under  contract  to  the  Wright  Research  Development  Center’s  Electro- 
Optic  Sources  Branch,  Wright  Patterson  AFB,  Ohio.  Characterization  of  second  harmonic 
generation  using  these  materials  is  required  for  comparison  of  their  relative  figures  of  merit. 

These  figures  of  merit  include  conversion  efficiency,  angular  acceptance,  phase  matching 
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angle  and  damage  threshold.  The  measurements  taken  will  provide  an  estimate  of  their  useful¬ 
ness  in  high  averag ;  power  frequency  conversion  from  repetitively  pulsed  laser  sources. 

In  this  research,  theoretical  estimates  of  the  conversion  efficiency  will  be  made  using 
published  values  of  the  nonlinear  coefficients,  phase  matching  angles  and  damage  limits  along 
with  measured  characteristics  of  the  pump  source  in  order  to  specify  the  physical  parameters  of 
the  experiments.  Deviations  in  the  actual  measurements  will  highlight  the  differences  between 
the  present  and  published  values  as  well  as  the  shortcomings  in  the  mathematical  models 
employed. 

The  sample  of  AgGaSe2  (9mm  x  9mm  x  29mm),  anti-reflection  coated  for  10.6  and  5.3  pm, 
will  be  characterized  in  terms  of  the  phase-matching  angle,  angular  acceptance  and  conversion 
efficiency  for  the  focussed  and  non-focussed  cases.  The  uncoated  sample  of  ZnGeP2 
(6.2mm  x  6.2mm  x  1 1mm)  will  be  characterized  in  terms  of  the  phase-matching  angle  at  a 
given  temperature,  angular  acceptance  and  conversion  efficiency  for  the  focussed  case  only. 
The  repetitively  pulsed  C02  laser  source  provides  a  very  stable  TEM m  output  where  the  tempo¬ 
ral  profile  of  the  pulses  are  Gaussian.  The  measured  efficiencies  will  be  compared  with  calcu¬ 
lations  involving  the  marsured  crystal  properties  and  literature  values  of  the  nonlinear 
coefficient.  Finally,  damage  thresholds  for  AgGaSe2  will  be  investigated  using  two 
9mm  x  9mm  x  2mm  witness  samples  which  are  coated  on  one  side  only.  The  threshold  for 
damage  will  be  related  in  terms  of  peak  and  average  intensities  for  both  surfaces  of  these  sam¬ 
ples. 

This  document  is  organized  to  provide  the  reader  with  sufficient  background  to  appreciate 
the  mathematical  developments  as  well  as  the  intricacies  of  the  experiment  procedures  and 
results.  Chapter  II  introduces  the  theory  of  nonlinear  optics  and  the  application  to  the  phase 
matched  conditions  of  second  harmonic  generation.  The  next  chapter  involves  the  character¬ 
ization  of  the  C02  laser  source  and  both  nonlinear  crystals.  Theoretical  estimates  of  the  crys¬ 
tals’  performance  follows  in  Chapter  IV.  These  results  are  used  to  specify  the  experimental 
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set  ups  for  the  non-focussed  and  focussed  efficiency  measurements  and  are  enumerated  in 
Chapter  V.  The  results  and  discussion  are  presented  in  Chapter  VI  with  a  comprehensive  Sum¬ 
mary  and  Suggesrions/Recommendations  in  the  concluding  chapters.  Background  theory  of 
linear  optics  in  anisotropic  materials  is  found  in  Appendix  A.  Appendix  B  describes  a  program 
developed  to  numerically  evaluate  one  of  the  theoretical  estimates  of  second  harmonic  genera¬ 
tion  conversion  efficiency. 
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In  the  theory  of  nonlinear  optics,  the  material’s  response  to  an  incident  electromagnetic 
field,  or  the  instantaneous  polarization  of  the  medium,  will  be  expanded  to  include  both  the 
linear  and  nonlinear  effects.  In  considering  the  linear  theory  of  propagation  of  light  through  a 
medium,  its  speed  is  defined  by  the  index  of  refraction  of  the  material.  In  an  isotropic  media 
this  value  is  independent  of  both  the  polarization  of  the  light  and  its  direction  of  travel.  In  an 
anisotropic  medium  however,  the  incident  radiation  is  broken  into  two  orthogonally  polarized 
components  whose  speeds  (or  indices  of  refraction)  will  vary  with  direction.  The  magnitude  of 
the  variation  is  determined  by  the  material  properties  and  will  correspond  to  the  two  solutions 
of  the  Fresnel  equations  in  the  medium.  A  review  of  linear  anisotropic  theory  is  found  in 
Appendix  A. 

The  instantaneous  polarization  of  an  anisotropic  medium  may  be  expressed  in  a  Taylor 
scries  which  includes  the  linear  susceptibility,  .  In  tensor  notation,  this  appears  as, 

P,  ~  X<jEj  +  2£0xyi£/£4  +  +  ....  (2.1) 

where  P,  is  the  ith  component  of  the  instantaneous  polarization  and  £,  is  the  i'th  component  of 
the  instantaneous  field  (23:504).  The  linear  susceptibility  is  shown  as  the  first  term  in  Eq.(2.2), 
while  the  second  is  the  nonlinear  term  responsible  for  second  harmonic  generation  and  is 
generally  written  as  dIJt  =  EoX.,*-  In  this  development,  the  third  term  in  the  expansion, 
which  gives  rise  to  such  interesting  effects  as  phase  conjugation  and  third  harmonic  generation, 
can  be  considered  negligible. 

Note  the  second  order  term  is  dependent  on  the  presence  of  the  two  field  terms  which  will 
both  be  of  the  general  form  in  Eq.(2.2).  We  represent  the  first  field  at  (0,  as 


where  1,  j  and  k  refer  to  the  field  components  along  the  principal  axes  x,  y  and  z  respectively. 
The  field  E^1  may  be  written  similarly.  If  the  medium  is  nonlinear  in  response,  the  fields  will 
give  rise  to  polarizations  at  the  sum  and  difference  frequencies  of  the  input  fields,  />,C0,  '<D|±t0j(r). 
Restricting  our  attention  to  the  case  of  second-harmonic  generation  where  (tii  =  Ci^,  the  result¬ 
ing  cross  terms  of  the  complex  conjugates  will  vanish  end  the  second  order  nonlinear  polariza¬ 
tion  will  be 


P2*  =  d  kE*E?  (2-3) 

for  Type  I  phase  matching  (to  be  covered  later).  Note  that  the  degeneracy  due  to  the 
frequencies  being  equal  is  compensated  by  a  factor  of  one  half  in  Eq.(2.2)  above.  The  second 
order  nonlinear  susceptibility  tensor  considered  here  is  not  a  frequency  dependent  function  and 
is  a  material  dependent  property  often  measured  by  second  harmonic  generation  experiments. 
Since  no  physical  significance  can  be  assigned  to  changing  the  order  of  j  and  k  in  Eq.(2.3) 
(23:507),  the  second  order  tensor  may  be  contracted  to  1 8  independent  elements  by  assuming 
the  convention  below.  Note  in  Eq.(2.3)  that  although  the  contraction  of  the  indices  does  not 
change  the  polarization,  the  order  to  the  field  terms  themselves  in  the  tensor  notation  is  very 
important. 


xx  =  1,  yy  =  2,  zz  =  3,  (2.4) 

yz  -  zy  -  4,  xz  -  zx  =  5,  xy  =  yx  =  6. 

The  resulting  dv  forms  a  3  x  6  matrix  that  operates  on  the  crossed  terms  of  the  incident  field  for 


Type  I  phase  matching  as 


du 

dx  3 

d\4 

d\i 

d> : 

Py 

= 

dn 

di  3 

du 

dv 

du 

^31 

dn 

dyy 

dy 4 

dys 

dy*; 

z 

z] 

E) 

2E,Ey 
2  E,EX 


(2.5) 
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Both  of  the  crystals  considered  in  this  analysis  fall  into  the  crystal  point  group  42 m  where 

the  only  nonzero  components  are  du,  d^,  and  d*,  where  d14  =  We  now  have  a 
description  of  the  physical  properties  of  the  material  given  the  direction  of  propagation  and 
polarization  of  the  input  field. 

Interaction  of  Light  and  the  Nonlinear  Medium 

Below,  we  will  develop  a  set  of  expressions  which  quantify  the  interaction  between  the 
input  radiation  frequencies  and  the  medium  in  setting  up  time  dependent  polarizations  along 
the  three  coordinate  axes.  Consider  Maxwell’s  equation  (Eq.(A.8  -  1 1))  written  in  the  form 
that  incorporates  the  polarization  explicitly 

—  a  -  _  (2.6) 

V  XJH  =  J  +  t(£o £  +  P) 

at 


V  x£ 


(2.7) 


where  P  =  +  Pnl  and  the  linear  susceptibility  is  treated  as  a  scalar  quantity  and 

the  nonlinear  polarization  term  is  taken  as  the  second  term  in  Eq.(2. 1).  Following  the 
development  in  Yariv  (23:516)  we  rewrite  Eq.(2.7)  as 

a  -  dPNL  (2.8) 

V  x  H  =  o£  +  t-zE  +  — r~~ 

dt  dt 


where  the  nonlinear  polarization  term  may  be  written  =  2dvkEjEk.  The  prime  on  the 


second  order  susceptibility  indicates  that  this  is  a  variable  transformed  to  coincide  with  the 
direction  of  propagation  with  respect  to  the  principal  coordinate  system.  It  is  known  as  the 
effective  nonlinear  coefficient,  d ,0.  The  current  density  has  been  rewritten  in  terms  of  the 
conductivity,  a,  and  may  be  regarded  as  a  frequency  dependent  loss  term.  It  may  or  may  not 
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be  included  in  developments  to  predict  overall  second  harmonic  generation  efficiency.  Taking 
the  curl  of  both  sides  of  Eq.(2.8)  and  replacing  VxW  with  Eq.(A.l  1)  along  with  the  vector 
identity  V  x  V  x  £  =  VV  •  E  -  V2£  ,  yields 


V2£ 


BE  ¥E  d2- 

*°ar  +  +  M(v/V 


(2.9) 


If  we  restrict  the  problem  to  the  one  dimensional  case  of  propagation  along  the  z-axis,  the 
partial  derivatives  with  respect  to  the  x  and  y  axes  may  be  neglected.  Considering  the 
solutions  along  the  three  coordinate  axes,  where  0),,  and  0)3  represent  the  frequencies, 
the  corresponding  fields  in  the  form  of  travelling  waves  appear 

£,  ’(z,0  =  ^[Eu(z)e‘<a'‘  v>  +  c.c.] 

fiiVO  «  5[£„(i)e'lv'v,  +  c.c.]  (2-10) 

£,”’(*■')  -  |[£l/(2)<l"">'  V’  +  C.C.] 


where  the  subscripts  /,  j,  k,  refer  to  the  cartesian  coordinate  system  (and  may  take  on  the  values 
x  and  y)  whose  orientation  with  respect  to  the  principal  axes  is  accounted  for  by  dlIf.  The 
specific  form  of  the  /th  component  of  the  nonlinear  polarization  term  at  frequency  0)!  given  the 
form  of  the  travelling  waves  in  Eq.(2.2)  is 


d'vlEh(z)El(z)e 


(2.11) 


The  conservation  of  energy  condition,  (O3  =  to,  +  a>2 ,  is  valid  for  all  time  t,  and  the  operations 


specified  in  Eq.(2.9)  may  be  carried  out  noting  that  the  second  derivative  terms  with  respect  to 
the  spatial  coordinate  z  are  very  small.  The  wave  equation  above  must  be  satisfied  by  each  of 
the  three  field  components  separately,  where  the  form  of  the  nonlinear  polarization  is 
Eq.(2. 11),  the  result  is 
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dEu 

dz 


-  /Mo  r  /  Mo  ,■  _  -<<*3  “  ~  *i> 

■  “W  i>£j'£af 

r^,-  ,  -%  /  Mo  .•  r  „• 

2ye2^2A  +  l(^y  ^d^Eyfi^e 
°3  /Mol,  /  Mo  .•  _  r 

"  2  V  £3 £3y  1  ^  \  ej  d^E"E^e 


(2.12) 


This  set  of  equations  accounts  for  the  interaction  between  the  three  fields  within  the  media  and 
covers  all  second  order  interactions.  An  alternate  development  would  attack  the  problem  from 
the  aspect  of  the  amount  of  energy  delivered  to  the  medium  by  each  field  as  a  function  of  time 
and  space.  Second  harmonic  generation  is  the  special  case  where  the  two  input  frequencies  are 
identical.  These  set  up  a  nonlinear  polarization  at  twice  their  value  leading  to  a  dipole 
radiation  term  at  CO3. 


Ltj  •!  IM’  iKlJMiliU  Ll4i'  1 LM 


Consider  the  last  of  Eqs.(2.38)  for  the  case  of  no  absorption  at  the  third  frequency,  o3  =  0, 


and  identical  input  fields  such  that  and  (03  =  20),.  In  the  initial  development,  we  can 

also  assume  that  the  fields  at  the  input  frequencies  are  not  significantly  diminished  by  conver¬ 
sion  to  the  second  harmonic  and  thereby  drop  their  dependence  on  z  and  write 


\  “h'yf^d'vkEuEli  exp[/ -  k\,}  -  Oz } 


(2.13) 


The  factor  of  one  half  comes  from  the  development  of  d ^  previously  for  the  case  of  second 

harmonic  generation.  is  the  propagation  constant  for  the  beam  at  fa),  which  is  polarized 
along  the  direction  i.  Recall  that  the  anisotropic  medium  will  separate  the  incident  radiation 
into  two  eigenpo'arization  modes.  In  general,  the  input  frequency  will  be  a  linear  superposi¬ 
tion  of  the  two  where  the  superscripts  on  k  indicate  that  the  magnitudes  will  depend  on  these 
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polarizations  states. 

Assuming  a  crystal  of  finite  length,  L,  integrating  Eq,(2.13),  for  the  case  of  no  initial  second 
harmonic,  £3,( 0)  =  0 ,  the  result  is 

.1  (214) 

E»  L>  “  ~l2<vy  iMc 

where  Ak  =  k^  -  kl']  -  k[k)  and  is  often  called  the  phase  mismatch.  The  repeated  indices 


indicate  summations  and  in  general  there  will  be  four  separate  terms  in  Eq.(2.14).  In  order  to 
eliminate  the  complex  terms  in  this  equation  and  reduce  it  to  a  more  tractable  form,  we  can 
generally  need  only  consider  the  two  cross  terms,  and  dju  0  *  k).  They  will  domi¬ 
nate  the  nonlinear  optical  tensor  due  to  either  to  the  phase  matching  conditions  or  large 
nonlinear  coefficients.  Multiplying  both  sides  by  the  complex  conjugate,  recognizing  the 
symmetry  of  the  nonlinear  tensor  and  the  trigonometric  form  of  the  last  term,  the  equation  may 
be  rewritten  in  terms  of  real  variables 


where  i  =  x. 


£‘(L)£Jy(L) 
k  =  y  (or  i  =  y. 


4po 

£3 


^{djElE^2 


sin2^A kL 


(2.15) 


k  =  x  ),  and  there  is  no  summation  over  the  repeated  indices. 


Also  note  the  fundamental  frequency  has  been  included  in  this  equation. 

In  terms  of  power  over  a  specified  area,  we  may  use  the  expression  derived  from  the  time 
averaged  Poynting  theorem 


Inserting  Eq.(2.14),  this  expression  becomes 


(2.16) 


T  _  /if!! 

rea  ~  2  Ve,“'M")£"£:'*L  a 


P 

J 

area 


sin2^  A kL 


(2.17) 


(W 


Finally,  in  terms  of  intensity,  for  the  case  when  input  power  is  divided  equally  between  the  two 
input  fields  at  the  same  fiequency,  the  result  is 
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where  e, 


p( 2a) 


£3  =  Eo^2  ■ 


Mo 

Eo 


v>]{d'jL‘ 


n 


[pfayy  ^  sin^~ AkL 

area  J  QmlJ 


(2.18) 


This  expression  yields  the  amount  of  power  developed  in  the  second  harmonic  for  a 
continuous  input  beam  at  the  fundamental  frequency  with  no  depletion  or  absorption.  It 
inherently  assumes  the  intensity  level  is  constant  across  the  beam  and  that  the  there  is  no 
deviation  of  the  two  fields  over  the  length  of  the  crystal  tno  walkoff).  In  Chapter  IV,  more 
rigorous  developments  will  incorporate  some  or  all  of  these  additional  considerations. 

In  all  cases,  the  amount  of  second  harmonic  generated  is  limited  by  the  finite  coupling  of 
the  input  field  to  the  dipole  nature  of  the  crystalline  material,  and  to  the  manner  in  which  the 
microscopic  fields  add  to  create  the  output  beam. 

Phase  Matched  Conditions 

For  the  maximum  power  to  be  generated,  the  phase  mismatch  in  Eq.(2.18)  must  be  zero,  i.e. 
Ak  s  0.  If  Ak  *  0,  the  field  generated  at  one  point  within  the  crystal  will  be  out  of  phase  with 
that  of  another,  resulting  in  an  interference  between  the  wavefronts  described  by  the  sinci2  term. 
This  may  arise  for  a  number  reasons  including  deviations  in  angle,  temperature  fluctuations, 
divergence  of  the  input  beam  and  a  finite  wavelength  bandwidth  of  the  source.  Recalling  the 
definition  of  Ak  in  terms  of  the  indices  of  refraction  (assuming  co,  =  0^  and  dropping  the  sub¬ 
scripts) 


Ak  =  k™  -  2k™  =  -(n^-n™) 

c 


(2.19) 


Since  the  index  of  refraction  generally  increases  with  frequency,  we  would  ideally  seek  a 
situation  where  the  difference  between  the  two  indices  is  a  minimum. 

If  we  note  the  representation  of  the  index  surfaces  for  the  fundamental  and  the  second 
harmonic  in  a  uniaxial  crystal,  Figure  1,  we  see  that  two  of  the  ellispes  intersect  at  a  certain 
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angle  denoted  as  0^.  These  ellipses  are  plotted  for  the  case  of  a  negative  uniaxial  crystal  ( n,  < 
n0)  and  it  can  be  seen  that  the  ordinary  index  of  refraction  for  the  fundamental  is  equal  to  the 
extraordinary  index  of  refraction  of  the  second  harmonic  at  this  angle. 


Figure  1 .  Phase  matching  condition  in  a  negative  uniaxial  crystal  shown  in  the  xrx3  plane 
(22:213).  0„  is  the  phase  matched  angle. 

If  we  were  to  propagate  a  beam  at  this  angle  with  respect  to  the  principal  axis  x3  inside  the 
crystal,  polarized  perpendicular  to  this  axis  (recall,  the  polarizations  are  susceptible  to  the 
index  of  refraction  along  the  axis  parallel  to  the  electric  field,  here  the  ordinary  index),  the  ort¬ 
hogonal  polarization  at  the  second  harmonic  would  have  an  index  of  equal  magnitude.  This  is 
known  as  the  phase  matched  angle. 

Figure  1  is  an  example  of  Type  I  phase  matching  for  second  harmonic  generation  in  a  nega¬ 
tive  uniaxial  crystal.  The  two  input  fields  possess  the  same  polarization,  along  x,  (or  x2),  and 
the  output  is  orthogonal  to  both.  If  the  two  fields  were  input  at  orthogonal  polarizations,  and 
the  output  was  parallel  to  either  of  them,  it  would  be  called  Type  II.  Type  III  is  very  similar, 
where  the  output  would  fall  on  the  other  polarization.  Both  of  the  crystals  characterized  is  in 
this  analysis  are  capable  of  Type  I  phase  matching  only  given  the  fundamental  wavelength  of 
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10.6  ^tm- 


We  may  rewrite  the  equation  for  the  extraordinary  index  of  refraction  in  the  uniaxial  crystal 

in  a  parametric  form  related  to  tire  angle  of  propagation  with  respect  to  the  optic  axis  (x}),  6, 

for  either  of  the  two  frequencies  as  follows 

1  _  cos2  8  sin2  8  (2.20) 

n.(  8)2  n]  n) 


For  the  case  of  Type  I  phase  matching  in  a  negative  uniaxial  crystal,  we  can  equate  Eq.(2.20), 
written  for  2 (0,  with  to  yield  (22:213) 


sin  8.  = 


(O'2  -  (O'2 
(Of2  -  (Of2 


(2.21) 


where  0m  is  the  phase  matched  angle.  For  the  case  of  positive  uniaxial  crystals  this  becomes 


(Of2  -  (Of2 

sin  e_  =  - - - 


(2.22) 


(Of  -  (Of 


The  coupled  wave  equations  bear  out  the  fact  that  the  input  at  the  fundamental  and  the 
output  at  the  second  harmonic  are  orthogonally  polarized  for  Type  I  phase  matching,  but  a 
mechanical  argument  related  to  the  dipole  oscillators  is  in  order  (Note  Figure  2).  The 
calculated  solutions  for  the  specific  cases  of  the  materials  to  be  characterized  will  be  carried 
out  in  the  next  chapter. 

The  fundamental  polarization  is  parallel  to  the  dipole  polarization  and  induces  an  oscillation  in 
this  plane  at  twice  the  natural  frequency.  Since  dipole  oscillators  radiate  a  polarization  perpen¬ 
dicular  to  their  axis  and  since  the  two  frequencies  travel  at  the  same  speed  (same  index) 
through  the  medium,  a  coherent  constructive  field  is  built  up.  In  all  other  directions  the 
oscillations  produce  fields  that  destructively  interfere  and  are  not  observed. 
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Opote  Oecilelor 


Figure  2.  Dipole  oscillator  diagram,  showing  excitation  parallel  to  the  input  polarization 
direction  and  the  induced  oscillation  producing  the  output  at  the  orthogonal  polarization. 

Figures  of  Merit 

With  a  basic  understanding  of  the  theory  of  second  harmonic  generation  in  hand,  it  is  now 
possible  to  begin  a  comparison  of  various  materials  in  terms  of  their  relative  performance  char¬ 
acteristics. 

The  effective  nonlinear  coefficient  is  dependent  upon  the  nonlinear  optical  tensor  and  the 
direction  of  propagation  through  the  medium.  The  specifics  of  its  derivation  arc  covered  in 
Chapter  IV. 

The  angular  sensitivity,  p#,  relates  the  rate  with  which  the  second  harmonic  rolls  off  with  a 
change  in  angle.  These  are  found  by  retaining  the  First  two  terms  in  Taylor  series  expansion  of 
Eq.(2.19)  with  Eq.(2.20)  inserted  for  n,(Q)  given  Type  I  phase  matching. 


M(0  )L  =  2Pe(0-9J 

(2.23) 

-oi L  ,  (O'2-  (O'2 

- sin(2(L) - - - ,  negative  uniaxial 

c  2(0 

(2.24) 

-oiL  .  ,  (O’3 "  (O'  . 

- — -sin(20„) - - - .positive  uniaxial 

c  2(0 

(2.25) 
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The  angular  acceptance  at  the  phase  matching  condition  can  be  related  to  the  full  width  at 
half  the  maximum  values  of  the  sine 2  term  in  Eq.(2.18).  The  anguiar  bandwidth  can  be  found 


as  twice  the  angle  in  radians  which  satisfies  the  parametric  relation 

sin2  Pe(0  ~  9»i)  1 


p(2u)  x 


(2.26) 


(Pe(6-ej]2  2 

The  walkoff  angle  between  the  fundamental  and  second  harmonic  in  a  uniaxial  crystal  may 
be  expressed  as 


tanp  =  ^K(0)]2 


_1 _ 1_ 

n\  nl 


sin  20 


(2.27) 


where  n,(6)  is  found  from  Eq.(2.20)  and  nE  and  nQ  are  found  from  the  Sellmeier  equations  for 


the  material.  In  a  negative  uniaxial  crystal,  the  second  harmonic  will  deviate  from  the 
propagation  direction  by  this  angle  and  for  a  positive  uniaxial  crystal,  it  will  be  the  fundamen¬ 
tal. 

The  measured  conversion  efficiency,  f>t2“>/PM,  will  be  dependent  upon  these  parameters  as 
well  as  the  physical  arrangement  of  the  experiment.  It  is  a  sensitive  function  of  the  different 
terms  contributing  to  the  phase  mismatch,  intensity,  focussing  and  the  crystal  properties  such 
as  absorption. 

This  review  of  the  basic  nonlinear  theory  provides  a  basis  upon  which  we  can  develop  an 
understanding  of  the  experimental  results. 
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III.  Experimental  Apparatus  and  Source  Characterization 


It  is  imperative  in  the  theoretical  predictions  of  the  conversion  efficiency,  as  well  as  the 
other  figures  of  merit,  that  the  values  used  to  describe  the  physical  nature  of  the  laser  beam  at 
the  fundamental  wavelength  and  the  crystal  be  as  accurate  as  possible.  With  these  parameters 
in  hand,  we  are  able  to  specify  the  optimum  experimental  arrangements  for  the  various  mea¬ 
surements.  In  addition,  deviations  between  the  experimental  and  theoretical  values  will  point 
out  shortcomings  in  the  calculations  or  errors  in  the  experimental  procedure.  For  this  reason, 
the  experimental  results  of  the  source  and  crystal  measurements  which  are  used  to  derive  the 
estimates  of  efficiency  in  the  following  chapters  are  presented  below. 

Detectors  and  Stages 

A  number  of  different  detectors  were  used  in  this  series  of  experiments  which  will  be 
reviewed  here  for  the  sake  of  reference  in  later  sections. 

Continuous  wave  power  measurements  and  short  duration  repetitively  pulsed  energy  mea¬ 
surements  were  conducted  using  Scientech  Model  365  one  inch  diameter  surface  absorbing 
disc  calorimeters.  The  analog  meters  provided  a  0.001  to  10  Watt  (or  Joule)  operating  range. 
At  the  lower  power  or  energy  ranges,  isothermal  enclosures  were  used  to  minimize  the  fluctua¬ 
tions  in  the  measurements  due  to  air  currents  in  the  laboratory.  The  surface  absorbers  were 
coated  with  Nextel  black  paint  which  provided  a  flat  spectral  response  for  all  wavelengths 
under  study  in  the  experiment. 

Temporal  measurements  of  the  repetitively  pulsed  fundamental  and  second  harmonic  were 
made  using  a  pair  of  liquid  nitrogen  cooled  resonant  cavity  HgCdTe  detectors.  The  detectors 
and  preamplifiers  were  originally  designed  for  heterodyned  measurements  and  have  a  band¬ 
width  on  the  order  to  100  MHz.  With  the  resonant  cavity  design,  the  detectors  were  relatively 
insensitive  to  wavelengths  other  than  those  of  interest  The  outputs  were  captured  with  a  200 
MHz  bandwidth  digital  oscilloscope. 
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An  Infrared  Systems  Model  700  room  temperature  pyroelectric  vidicon  camera  was  used  to 
record  the  spatial  profile  of  the  fundamental  and  second  harmonic  beams.  The  vidicon  tube 
(210x150  effect  elements)  is  sensitive  from  2pm  to  20pm  and  has  a  dynamic  range  of  0  to 
256  (8  bits,  7  useful).  The  camera  employs  a  25  Hz  chopper  to  modulate  the  input  beam  for 
the  pyroelectric  array  and  an  electronic  suppression  of  this  frequency.  The  camera  is  con¬ 
trolled  via  a  Big  Sky  Software,  Beamcode,  IBM  computer  based  system.  Data  acquired  with 
this  camera  was  useful  for  relative  measurements  only.  It  proved  notoriously  unreliable  in 
performing  absolute  measurements  on  the  continuous  or  repetitively  pulsed  laser  source.  It 
was  able  to  provide  a  basis  forjudging  beam  and  mode  quality  of  the  laser  in  any  case. 

The  mechanical  stages  used  for  angular  adjustments  in  these  experiments  provided  sub¬ 
minute  of  arc  accuracy  and  could  be  arranged  to  control  any  desired  axis  of  rotation.  Relative 
adjustments  with  linear  translation  and  elevation  stages  were  accurate  to  within  0.1  mm  along 
all  axes  of  interest. 

£Q2  laseLSpurfg_&  Output  Characterization 

The  carbon-dioxide  laser  source  used  in  these  experiments  provided  continuous  wave  and 
repetitively  pulsed  output  from  an  actively  controlled  frequency  stabilized  local  oscillator/ 
main  amplifier  design.  A  schematic  diagram  of  the  optical  system  is  shown  below  in  Figure  3. 
The  laser  is  actually  a  comprised  of  two  separate  lasers  for  frequency  control  and  power  stabil¬ 
ity  followed  by  a  single  pass  four  channel  waveguide  amplifier.  Together  they  provide  a  high 
average  power  repetitively  pulsed  output  beam  with  good  mode  quality  and  frequency  stability. 

On  the  lower  left  side  of  Figure  3,  the  local  oscillator  is  a  very  short  (14.5  cm)  grating  tuned 
cavity,  135  //mm,  which  is  tuned  for  the  center  of  the  P(20)  line  of  C02.  The  cavity  length  of 
the  oscillator  is  varied  to  maintain  a  resonance  at  60  Mhz  off  of  line  center.  At  this  offset,  the 
output  power  is  a  very  sensitive  function  of  cavity  length.  The  error  signal  generated  to  drive 
the  piezoelectrically  controlled  mirror  mount  comes  from  the  comparisoi  of  the  detector  out¬ 
put  near  the  local  oscillator’s  output  port  to  a  set  reference  level. 


I 


16 


Main  Amplifier 
Output 


Figure  3.  Schematic  diagram  of  the  C02  laser  source.  There  are  two  separate  lasers  present; 
the  grating  tuned  local  oscillator  and  the  electro-optically  modulated  amplifier  /  coupled  cavity 
laser.  The  output  from  the  EO-Modulated  laser  then  passes  through  a  four  channel  waveguide 
amplifier  to  the  beam  telescope  and  exits  the  system. 

Pan  of  the  local  oscillator’s  output  is  also  split  off  to  act  as  the  frequency  standard  for  the 
rest  of  the  system  while  the  rest  exits  the  housing  as  a  horizontally  polarized  beam. 

The  second  laser  provides  the  basis  for  die  electro-optically  modulated  output  via  the  intra¬ 
cavity  CdTe  modulator.  The  folded  amplifier  is  approximately  70  cm  in  length  where  the 
intermediate  mirror  is  mounted  on  a  piezoclectrically  controlled  stage.  The  output  passes  into 
the  coupled  cavity  through  the  45%  reflective  ZnSe  window  prior  to  the  CdTe  modulator.  At 
this  point  the  beam  is  horizontally  polarized  and  after  passing  through  the  modulator  twice,  the 
beam  is  vertically  polarized  and  is  sheared  off  by  the  germanium  polarizer/coupler. 

A  small  amount  of  the  beam  leaks  through  the  high  reflector  in  the  coupled  cavity  to  the 
mixed  with  the  output  from  the  local  oscillator.  The  two  beams  are  detected  via  heterodyning 
where  a  constant  offset  frequency  of  the  local  oscillator  provides  the  error  signal  to  the  inter- 
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mediate  cavity  mirror’s  stage.  The  long  term  output  frequency  stability  is  maintained  to  within 
30  MHz  of  the  P(20)  line.  The  optimum  output  power  is  maintained  by  the  output  of  the 
pyroelectric  detector  after  the  four  channel  amplifier.  A  small  dither  is  placed  on  the  rear  cou¬ 
pled  cavity  mirror  and  the  error  signal  provides  the  signal  for  the  piezoelectric  stage. 

The  beam  goes  through  a  series  of  mirrors  prior  to  passing  through  the  four  channel 
waveguide  amplifier  and  existing  the  system  through  a  Cassegrain  beam  expanding  telescope. 
Note  in  Figure  3,  the  dotted  line  traces  the  path  of  the  main  output  beam  from  its  origin  within 
the  electro-optically  modulated  cavity  tlirough  the  output  port. 

The  vertically  polarized  output  beam  has  an  average  power  of  approximately  8  Watts  CW 
and  7  W  when  modulated  at  30  kHz.  The  output  mode  is  TEMqq  except  for  a  slightly  elliptic 
nature  caused  by  the  beam  combining  elements  within  the  cavity.  Two  other  modulation  fre¬ 
quencies  were  available,  50  and  100  kHz,  but  not  often  used  due  to  lower  peak  powers  (though 
higher  average  powers). 

Polarization  of  the  Laser  Source.  The  polarization  of  the  main  amplifier  and  the  local 
oscillator  were  measured  by  rotating  an  infrared  polarizer  of  known  polarization  through  a  180° 
arc  and  noting  the  variation  in  power.  According  to  Malus’  Law,  the  variation  can  be  fit  to  a 
cosine  of  the  angle  of  orientation  of  the  polarizer  with  respect  to  input  beam.  Both  of  the 
beams  were  attenuated  to  approximately  10  mW  in  order  to  avoid  damage  to  the  very  fine  grid 
of  gold  wires  held  between  a  pair  of  germanium  windows  that  made  up  the  polarizer. 

According  to  the  design  of  the  overall  laser  system,  the  two  output  beams  should  be  ortho¬ 
gonally  polarized.  As  can  be  seen  in  Figure  4,  the  polarization  of  the  local  oscillator  is  approx¬ 
imately  15°  off  of  horizontal.  This  does  not  pose  a  problem  in  the  heterodyning  detection 
inside  the  cavity  because  the  two  beams  at  that  point  remain  in  the  same  plane  of  polarization. 

The  output  from  the  main  amplifier  is  nearly  vertically  polarized  where  any  small  deviation 
should  not  significantly  effect  the  amount  of  properly  polarized  power  delivered  to  the  crystal. 
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Figure  4.  Polarization  of  the  two  beams  f't  to  a  cosine. 

Characterization  of  the  Beam  Divergence  and  the  Focussed  Soot  Sizes.  As  noted  in 
the  development  of  the  phase-matching  requirements  for  second  harmonic  generation,  the 
phase  mismatch,  A k,  is  a  function  of  the  deviation  of  the  input  beam  from  the  optimum  angle. 
This  can  be  due  to  mechanical  position  or  to  the  divergence  of  the  laser  as  it  passes  through  the 
crystal.  At  the  phase-matched  angle,  the  divergence  will  dominate  the  phase-mismatch  (espe¬ 
cially  in  the  focussed  case). 

Since  all  of  these  experiments  were  conducted  near  the  output  port  of  the  laser, 
measurements  of  the  divergence  would  have  to  take  into  account  the  fact  that  we  are  in  the 
near- field.  For  this  wavelength,  the  far-field  is  a  least  10  meters  from  the  output  port  of  the 
laser.  Thus,  measurements  of  the  beam  profile  were  taken  at  various  distances  from  the  laser 
and  numerically  fit  to  yield  an  analytic  description  of  the  divergence  of  the  beam. 

The  complex  beam  parameter,  which  describes  the  Gaussian  beam  as  function  of  wave¬ 
front’s  radius  of  curvature  and  the  spot  size  is  given  as, 
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(3.1) 


1  _J_  Ao 

q(z)  R(z)  nnu?(z) 

This  paiameter  can  be  used  to  describe  the  propagation  of  a  Gaussian  beam  through  an  optical 
system  with  the  appropriate  ray  transfer  matrices. 

In  this  particular  arrangement,  we  assume  the  location  of  the  beam  waist  to  lie  inside  the 
cavity  (near  the  output  port  of  the  waveguide  amplifier),  this  is  where  the  radius  of  curvature 
goes  to  infinity,  and  the  complex  beam  parameter  reduces  to 

]_  _  (3.2) 

<7o  nnu$ 

We  need  only  to  build  the  appropriate  set  of  ray  transfer  matrices  which  describe  the  propaga¬ 
tion  of  the  beam  from  the  waist  to  points  outside  the  laser.  All  of  the  intervening  optics 
between  the  beam  waist  and  the  output  port  can  be  disregarded  because  we  will  fit  the 
experimental  data  with  an  analytic  description.  The  spot  size  as  a  function  of  distance  from  the 
laser  is  described  by 


where  z  =  z2  +  Offset  The  Offset  is  the  constant  distance  from  the  output  window  of  the 

laser  back  to  the  location  of  the  beam  waist,  iv0,  and  z2  is  the  distance  from  the  laser  to  the 
point  of  interest. 

Measurements  of  the  continuous  wave  beam  profiles  were  taken  at  different  distances  from 
the  laser,  and  the  waist  radius  was  determined  to  be  u>0  =  1.019mm  where  Offset  =  0.31m. 
For  the  case  or  ’he  repetitively  pulsed  beam,  the  waist  radius  was  extracted  from  the  measure¬ 
ments  of  the  focussed  beam  spot  sizes. 

The  spot  size  formed  by  the  focussing  element  must  be  characterized  accurately  in  order  to 
calculate  the  intensity  (irradiance)  levels  of  the  beam  as  a  function  of  the  distance  from  the 
lens.  In  this  set  of  measurements,  the  change  in  power  observed  at  a  fixed  positior  was 
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recorded  as  a  knife  edge  is  scanned  across  the  horizontal  axis  of  the  beam.  Having  accurately 
calibrated  the  motion  of  the  knife  edge  as  a  function  of  time,  the  spot  size  at  that  distance  from 
the  lens  may  be  deduced  from  rate  of  fall  in  the  observed  power.  The  experimental  layout  for 
this  method  is  shown  in  Figure  5. 

Assuming  the  focussed  beam  is  a  two  dimensional  Gaussian  beam,  the  spatial  intensity 
distribution  may  be  described  by 


2P0  ( 

I(x,y)  =  — ^expl- 


2x2  +  2 y 

1F~ 


A 


(3.4) 


where  x  and  y  are  the  axes  chosen  perpendicular  to  the  beam,  ons  the  spot  size,  and  P0  is  the 

total  laser  power  (20:2809).  If  we  scan  the  beam  horizontally  with  the  razor  blade,  the  power 
transmitted  will  be  given  by 
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(3.5) 


Integrating  out  the  dependence  along  the  vertical  axis  of  the  beam,  assumed  here  to  be  the 
y-axis,  the  equation  becomes 


P(x) 


(3.6) 


which  takes  on  the  form  of  the  error  function. 

As  the  knife  edge  scans  across  the  beam,  the  power  recorded  will  drop  slowly  at  first  as  the 
blade  intersects  the  edge  of  the  pattern  and  rapidly  decrease  as  it  approaches  the  peak  of  the 
pattern.  If  we  are  able  to  easily  locate  the  90  %  and  10  %  of  total  average  power  points  in  the 
decay  of  the  signal  recorded,  Eq.(2.6)  above  may  be  written  in  terms  of  this  difference  in  time 
(or  distance).  From  this  distance  we  can  easily  extract  the  1/e  radius  of  the  beam  by 

1  (3.7) 

a/  =  0.7803  — 

’  6U4e 
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M.  is  the  distance  on  the  chart  paper  it  took  for  the  power  to  drop  from  the  90%  to  the  10% 

level,  is  the  chart  paper  speed  and  vbUde  is  the  speed  of  the  blade  as  it  cuts  through  the 
beam. 

Note  in  the  diagram  of  the  experimental  apparatus,  the  razor  blade  is  mounted  on  the  stem 
of  an  electrically  driven  micrometer  stage  with  0.001  mm  incremental  distance  readout  set  to  a 
speed  of  0. 1 1  mm/sec. 


Figure  5.  Experimental  setup  for  recording  the  focussed  spot  sizes. 

The  power  meter  used  to  record  the  average  power  of  the  beam  was  a  Scientech  1 "  disk 
calorimeter  with  a  0  to  100  mV  full  range  voltage  output.  The  voltage  source  had  a  much 
faster  response  time  than  the  analog  readout  (on  the  order  of  100  Hz)  and  thus  allowed  a 
Hewlett-Packard  1 175A  strip  chan  recorder  to  accurately  depict  the  fall  in  power  as  the  blade 
passed  through  the  beam  and  was  set  to  a  speed  of  15.0  ±  O.Olcm/sec. 

The  separation  between  the  lens  and  the  micrometer  stage  was  varied  in  1  min  intervals  and 
the  blade  was  scanned  through  the  beam  at  each  point.  With  a  series  of  spot  sizes  in  hand  it 
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was  only  a  matter  of  fitting  them  with  an  analytic  representation. 

Beginning  once  again  with  the  waist  at  a  known  position  inside  the  laser,  the  spot  size  as  a 
function  of  distance  from  the  focussing  lens  is  given  by 

8) 

where  A,  B,  C  and  D  refer  to  the  components  of  the  ABCD  ray  transfer  matrix  for  the  Gaussian 
beam.  It  is  the  result  of  the  multiplication  of  the  matrices  shown  below 

fA  B\  _  (\  Zj'l  1  ®  (\  Zj  +  Offset'j 

[c  »)  "  [o  lj  [-J  lj  lo  1  J 

This  last  matrix  in  Eq.(2.9)  accounts  for  the  propagation  over  a  distance  z2  from  the  front  of  the 
laser  to  the  lens,  through  a  focussing  element  of  focal  length/,  and  finally  a  propagation  over  a 
distance  z}  to  the  point  of  interest  (where  the  razor  cut  through  the  focussing  beam). 

The  fit  to  the  measured  points  is  found  by  varying  the/in  Eq.(2.9)  and  w*,  in  Eq.(2.10)  using 
the  known  values  of  the  distances  z2  and  z3.  The  results  of  this  analysis  are  shown  in  Figure  6. 
Using  this  same  beam  waist  radius,  similar  results  were  achieved  for  a  different  lens  of  a 
shorter  focal  length  as  shown  in  Figure  7. 
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Figure  6.  Focussed  spot  size  where  u%  =  0.7019mm  and  /=  12.5cm  as  a  function  of  dis¬ 
tance,  Zj,  from  the  lens. 


Distance  from  lent  |mj 

■  Measured  values  - Caloiated 

Figure  7.  Focussed  spot  size  where  =  0,7019mm  and /=  6.05cm  as  a  function  of  distance, 
Zj,  from  the  lens. 
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This  leads  us  to  consider  that  the  divergence  of  the  laser  in  the  repetitively  pulsed  mode  is 
larger  than  the  CW  mode.  Since  the  measurements  of  the  focussed  spot  sizes  were  conducted 
at  a  pulse  repetition  frequency  of  30  KHz,  and  the  calculated  spot  sizes  agreed  well  with  the 
data,  the  beam  waist  derived  in  this  analysis  will  be  used  in  the  theoretical  estimates  of  the  effi 
ciency.  The  change  in  divergence  is  most  likely  attributable  to  the  electro-optic  modulator  in 
operation  (in  the  CW  mode  the  electric  field  applied  is  static). 

Pulse  Repetition  Formats.  The  CdTe  electro-optic  modulator  preceding  the  four  chan¬ 
nel  waveguide  amplifier  allows  for  three  high  repetition  rate  pulsed  formats  at  30,  50  and  100 
kHz.  The  modulator  provides  for  a  smooth  Gaussian  envelope  in  time  as  opposed  to  the  pulse 
shape  of  most  Q-switched  lasers.  The  temporal  profile  for  the  three  different  pulse  repetition 
frequencies  is  shown  in  Figure  8. 


Figure  8.  Temporal  format  of  the  repetitively  pulsed  source. 

The  temporal  profiles  were  recorded  with  the  liquid  nitrogen  cooled  resonant  cavity  HgCdTe 
detectors.  The  average  power  increases  with  the  higher  modulation  frequencies,  but  the  peak 
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powers  decrease.  Since  the  theoretical  estimates  of  efficiency  for  pulsed  sources  are  propor¬ 
tional  to  the  peak  irradiance  values,  most  of  the  efficiency  measurements  were  completed  using 
the  30  kHz  modulation  frequency. 

Frequency  Stability.  The  frequency  stability  was  characterized  by  two  methods.  The 
first  involved  the  observation  of  the  main  amplifier’s  output  through  an  infrared  grating  spec¬ 
trometer  designed  for  carbon-dioxide  laser  wavelengths.  Only  the  P(20)  line  at  10.6  pm  was 
observed.  The  second  method  involved  heterodyning  the  main  amplifier  and  local  oscillator 
beams.  The  constant  offset  frequency  of  60  MHz  was  observed  after  careful  alignment  of  the 
optical  system. 

Spatial  Profiles.  The  two-dimensional  spatial  profile  of  the  repetitively  pulsed  formats 
were  observed  using  the  pyroelectric  vidicon  camera  system  once  the  beam  had  been  attenu¬ 
ated  to  less  by  10  mW/cm2;  approximately  three  orders  of  magnitude  of  attenuation  were 
required  to  observe  the  unfocussed  beam  at  30  cm  from  the  output  port  of  the  laser.  A 
representative  spatial  profile  is  shown  in  Figure  9. 


Figure  9.  Spatial  profile  of  the  repetitively  pulsed  main  amplifier  at  30  kHz.  The  diagram 
represents  the  entire  field  of  view  of  the  infrared  vidicon  camera  where  the  effective  area  is  6 
mm  square. 
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The  beam  appears  to  be  Gaussian  ( TEMoo )  with  only  minor  local  deviations.  This  provides 
confidence  in  the  calculations  of  the  focussed  spot  sizes  since  the  theory  of  optical  matrix 
transformations  as  outlined  above  is  applicable  to  the  Gaussian  mode. 

Characterization  of  the  Crystals 

The  sample  of  AgGaSe2  was  provided  under  contract  to  the  Wright  Research  Development 
Center’s  Electro-Optic  Sources  Branch  by  Cleveland  Crystals,  Cleveland  OH.  The 
9mm  x  9mm  x  29mm  crystal  was  provided  with  a  two  layer  1/4  (PbF2  and  ThF2  based)  anti¬ 
reflection  coating  for  10.6pm  and  5.3pm  and  was  cut  from  the  boule  at  57.17°  with  respect  to 
the  optic  axis  (their  calculated  phase  matching  angle);  note  Figure  10. 

The  6.2mm  x  6.2mm  x  1 1mm  sample  of  ZnGeP2  was  an  uncoated  sample  which  was  cut  for 
non-critical  phase  matching  with  the  face  perpendicular  to  the  x-axis  (also  the  [010]  direction); 
note  Figure  1 1 . 


Figure  10.  AgGaSej  crystal  cut  for  Type  I 

phase  matching  at  57.17°.  Figure  11.  ZnGeP2  crystal  cut  for  Type  I 

phase  matching  at  90°. 

The  crystals  were  characterized  by  spectral  scans  as  well  as  direct  absorptivity  measure¬ 
ments  using  the  smaller,  low  power  beam  of  the  local  oscillator. 

Absorptivity  Measurements.  In  order  to  avoid  errors  that  would  be  introduced  under 
the  assumption  that  the  laser  source  possessed  a  long  term  amplitude  stability,  this  measure- 
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ment  was  conducted  in  three  separate  parts  (26:5).  Note  the  definitions  established  in  Figure 
12  for  the  components  of  the  incident  beam  of  intensity,  /,  and  the  various  reflected  and 
transmitted  portions. 


Figure  12.  Calculation  of  the  absorption/scattering  loss.  R  is  the  reflectivity  of  the  crystal 
surface  (assumed  to  be  equal  on  both  ends),  I  is  the  incident  intensity,  L  is  the  length  of  me 
crystal  and  a  is  the  loss  per  unit  length  of  the  material. 

We  seek  to  establish  an  experimental  method  which  will  remove  /  from  the  calculations  of 
the  absorptivity,  a.  From  the  three  measurements  taken  in  the  manner  as  snown  in  Figure  13, 
note  that  the  ratios  are  not  functions  of  the  intensity  of  the  incident  beam  and  will  thus  remove 
the  uncertainty  in  the  long  term  amplitude  stability  of  the  source. 

If  we  establish  the  ratio 
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Ratio  1  -  B/A 


Figure  13.  Three  sets  of  measurements  to  be  taken  in  determining  the  absorptivity. 


„  Ratio-1  ,,  _ .  2  , 

T  =  Raiio-2  =  0  '  S)eXP(-<Ii) 


(3.10) 


and  note  the  specific  definition  of  the  terms  that  make  up  Ratio-3  as 


Ratio- 3  = 


R  +  R(  1  -  /?)2exp(-2ctL) 

(1  -  R)2exp(-aL) 


(3.11) 


we  note  that  Eq.(2.10)  can  be  substituted  into  Eq.(2.1 1)  to  yield  a  parametric  equation  with 

respect  to  the  reflectivity  of  the  surface  of  the  crystal  R. 

_  (Ratio-3)T  (3.12) 

i  +  r2/(i  -  R)2 


It  is  a  straight  forward  matter  to  iteratively  solve  for  the  R  and  subsequently  solve  of  the 
absorptivity  by 


(3.13) 
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An  alternative  method  of  measuring  the  absorptivity  simply  involves  the  measurement 
shown  in  Figure  1 3  for  Ratio-3  given  that  we  have  a  reasonably  good  estimate  of  the  surface 
reflectivity.  In  general,  we  can  use  that  defined  by  the  previous  analysis  as  a  good  starting 
point.  This  second  method  involves  measuring  the  transmittance  of  the  crystal  directly  and 
applying  Eq. (2. 13). 

The  low  power  continuous  wave  output  of  the  local  oscillator  was  used  to  conduct  these 
measurements,  where  the  beam  was  apertured  down  such  that  the  physical  dimensions  of  the 
crystal  did  not  obscure  the  beam  and  lead  to  erroneously  high  values  of  the  absorptivity. 

The  results  of  the  measurements  of  the  two  crystals  is  shown  in  the  table  below.  Note  that 
the  measured  values  for  both  of  the  methods  agree  very  well.  The  error  is  rather  large  due  to 
the  uncertainty  in  the  surface  reflectivity. 


Table  1. 

Results  of  the  absorptivity  measurements 


_ A&GaSez _ 

Absorptivity  measured  by  method  1  0.0267  ±  0.0270  /cm 

Surface  reflectivity  of  the  sample  0.0055  ±  0.0005 

Absorptivity  measured  by  method  2  0.0271  ±  0.0038  /cm 


_ ZnQsPj _ 

0.68 1  ±  0.0650  /cm 
0.248  ±0.025 
0.661  ±0.0623 /cm 


Spectrographic  Measurements.  The  samples  were  scanned  with  a  Fourier  transform 
infrared  spectrometer  in  order  to  establish  a  relative  value  of  the  absorptivity  at  the  second  har¬ 
monic's  frequency  as  well  as  provide  an  additional  check  for  the  measurements  above. 
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Figure  14.  Infrared  transmission  spectrum  of  the  AgGaSe2  crystal. 

In  both  crystals,  the  transmission  (which  can  be  related  to  the  absorptivity)  is  relatively 
constant  across  the  region  of  interest  from  5.3pm  to  10.6pm.  With  the  laser  source  and  the 
crystals  to  be  studied  well  characterized,  we  are  able  to  continue  with  the  theoretical  estimates 
of  the  crystal  performance  and  the  experimental  designs. 
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To  begin  the  experiments  to  characterize,  and  hopefully  optimize,  the  performance  of  the 
crystals  used  in  the  second  harmonic  generation  of  10.6  pm  radiation,  a  number  of  calculations 
will  prove  useful  in  establishing  the  experimental  design.  Given  the  constraints  of  available 
equipment,  a  primary  concern  was  establishing  the  intensity  levels  of  the  non-focussed  and 
focussed  beam  at  the  face  and  center  of  the  crystal.  With  this,  we  are  able  to  estimate  the  effi¬ 
ciencies  and  compare  the  values  of  peak  and  average  power  to  published  damage  thresholds. 
They  also  provide  a  means  to  approximate  the  temperature  profiles  generated  during  exposure. 
This  would  provide  an  estimate  of  the  stresses  being  placed  on  the  crystal. 

Presented  below  are  the  various  mathematical  developments  for  the  calculation  of  conver¬ 
sion  efficiency  under  different  experimental  conditions.  Comparison  to  the  experimental  data 
are  carried  out  in  Chapter  VI. 

NpndeBlflgd plans  Wave  Approximation 

As  an  initial  estimate  of  the  conversion  efficiencies,  we  begin  by  analyzing  the  dependences 
noted  in  Eq.(2.18).  Written  in  terms  of  the  angular  sensitivity,  pe,  the  conversion  efficiency, 
p  (2«>)/ p  (“),  becomes 


(P(<0))V.  2sin2p,(8w-8) 
EoJ  (area  /  ^  (Po(e„  -  8))2 


n  =  2 


(4.1) 


where  d ^  —  dvi  is  the  effective  nonlinear  coefficient  from  the  previous  development. 


This  solution  was  derived  under  the  assumption  that  the  plane  wave  incident  on  the  crystal 
is  not  affected  by  the  development  of  the  second  harmonic  (non-depleted  input  approximation). 
The  wavefront  is  assumed  to  be  infinite  in  extent  and  have  a  constant  amplitude  in  both  time 
and  space.  The  efficiency  is  proportional  to  the  intensity,  P^/area,  and  to  the  square  of  the 
crystal  length,  L.  The  dependence  upon  the  direction  of  propagation  is  accounted  for  in  two 
separate  terms.  The  deviation  from  the  phase-matching  angle,  0m  ,  produces  the  characteristic 


33 


sine2  drop  off,  and  ^possesses  an  angular  dependence  determined  by  the  crystal  symmetry 
and  the  corresponding  nonlinear  optical  tensor. 

In  this  experiment,  we  have  control  over  some  of  these  variables,  such  as  the  intensity  levels 
and  the  angular  orientation  of  the  crystals,  but  little  control  over  the  constants  such  as  crystal 
length,  phase  matching  angle  and  For  this  reason,  we  will  have  to  determine  the  values  for 
these  constants  in  order  to  specify  the  appropriate  magnitudes  for  the  variables. 

Calculation  of  Beam  Parameters  ■  Peak  and  Average  Intensities 


Given  the  spot  size  of  the  laser  at  some  point,  we  are  able  to  estimate  the  peak  intensity 
(fluence,  or  irradiance)  of  the  spatial  profile  at  the  temporal  peak.  From  the  characterization  of 
the  laser  source,  we  are  able  to  approximate  the  pulses  to  be  Gaussian  in  space  and  time.  For 
the  case  of  the  30  kHz  repetition  frequency,  the  measured  FWHM  (full  width  at  half  of  the 
maximum  value)  is  36.4  ns  (see  Figure  8)  which  leads  to  a  characteristic  time  at  the  1/e  point 
of  the  pulse  of  t  =  21.86/u.  The  power  of  a  single  pulse  as  a  function  of  time  is 

Pit)  =  (4‘2) 

where  is  the  peak  power.  Given  a  value  for  the  average  power,  PAVC,  and  the  pulse 

repetition  frequency,  v,ip,  the  peak  power  may  then  be  written  (9:663) 


This  is  the  total  power  at  the  temporal  peak  of  the  pulse.  It  can  be  related  to  the  intensity 
given  a  specified  area  through  which  the  power  will  pass.  Given  the  equation  for  the  radial 
intensity  distribution  of  a  Gaussian  beam 

/(r,0)  =  (4'4) 

The  total  power  is  given  as  the  integral  of  this  function  over  all  space  and  is  equal  to  the  peak 

power  derived  above. 
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Thus,  the  peak  intensity  may  be  written  (15:302) 

,  P AVG 


(4.6) 


‘Ptak 


where,  in  general,  the  spot  size,  u>,  will  be  a  function  of  distance  from  the  laser  or  the  focussing 
element. 

Eq.(2.8)  yields  the  spot  size  as  a  function  of  the  distance  from  the  laser  with  or  without  a 
focussing  element  and  can  be  used  to  find  the  intensity  at  the  face  of  the  crystal.  In  order  to 
find  the  values  at  the  center,  we  must  use  the  transformation  matrix  resulting  from 
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where  the  ratio  of  the  indices  accounts  for  the  interface  of  the  crystal,  n2,  of  length  L.  Figures 
16  and  17  show  the  spot  sizes  and  the  peak  and  average  powers  for  the  unfocussed  case. 


Figure  16.  Spot  size  as  a  function  of  distance  from  the  output  port  of  the  laser  at  the  face  and 
center  of  the  AgGaSe2  crystal. 


Figure  17.  Peak  and  average  intensity  at  the  center  of  the  beam  as  a  function  of  distance  from 
the  output  port  of  the  laser  given  for  the  face  and  center  of  the  AgGaSe2  crystal. 

The  calculations  for  the  focussed  case  are  summarized  in  Figures  18  through  22  for  the  face 
and  center  of  the  two  crystals  for  two  different  distances  from  the  laser.  These  plots  are  nor¬ 
malized  for  an  average  power  of  1  Watt  at  a  pulse  repetition  frequency  of  30  kHz  (FWHM  = 
36.4  ns)  such  that  the  values  can  be  interpreted  given  the  experimental  setup. 
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Distance  from  lens,  z 3  [cm] 


Figure  18.  For  a  focal  length  of  12.5  cm,  the  plot  shows  the  spot  size  as  a  function  of  the 
distance  from  the  focussing  element,  z3,  at  the  face  and  center  of  the  AgGaSe2  crystal  for  values 
of  z2  (distance  from  the  laser  to  the  lens)  of  0.5  m  and  1  m.  Note  the  distance  z3  is  to  the  face 
of  the  crystal,  thus  as  the  crystal  is  moved  towards  the  focus  (away  from  the  lens)  the  beam 
waist  passes  through  the  center  of  the  crystal  first. 


Distant*  from  lens,  z3  [cm] 


Figure  19.  For  a  focal  length  of  6.05  cm,  the  plot  shows  the  spot  size  as  a  function  of  the 
distance  from  the  focussing  elemenf,  z3,  at  the  face  and  center  of  the  AgGaSe2  crystal  for  values 
of  z2  (distance  from  the  laser  to  the  lens)  of  0.5  m  and  1  m, 


Figure  20.  For  a  focal  length  of  12.5  cm,  the  plot  shows  the  peak  and  average  intensity  as  a 
function  of  the  distance  from  the  focussing  element,  z3,  at  the  face  and  center  of  the  AgGaSe2 
crystal  for  values  of  z2  (distance  from  the  laser  to  the  lens)  of  0.5  m  and  1  m.  These  values 
have  been  normalized  to  an  average  power  1  Watt  at  a  pulse  repetition  frequency  of  30  kHz. 


Distance  from  Ians,  z3  (cm) 


Figure  21.  For  a  focal  length  of  6.05  cm,  the  plot  shows  the  peak  and  average  intensity  as  a 
function  of  the  distance  from  the  focussing  element,  z} ,  at  the  face  and  center  of  the  AgGaSe2 
crystal  for  values  of  z2  (distance  from  the  laser  to  the  lens)  of  0.5  m  and  1  m.  These  values 
have  been  normalized  to  an  average  power  of  1  Watt  at  a  pulse  repetition  frequency  of  30  kHz. 
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Figure  22.  For  a  focal  length  of  12.5  cm,  the  plot  shows  the  peak  intensity  and  spot  size  as  a 
function  of  the  distance  from  the  focussing  element,  z3,  at  the  center  of  the  ZnGeP2  crystal  for 
values  of  z2  =  0.5  m  (distance  from  the  laser  to  the  lens).  These  values  have  been  normalized 
to  an  average  power  of  1  Wan  at  a  pulse  repetition  frequency  of  30  kHz. 

These  for  the  spot  sizes,  peak  and  average  intensities  may  then  be  applied  to  the  theoretical 
estimates  of  the  conversion  efficiency.  Note  in  Figure  18,  that  the  focussed  waist  radius  is  not 
changed  by  the  introduction  of  the  crystal,  only  its  position  is  moved  in  accordance  with  the 
change  in  slope  n,/n2. 

Phasematched  Condition  and  the  Effective  Nonlinear  Coefficient 

The  phase  matched  condition  is  a  value  which  can  be  derived  as  the  numerical  solution  to 
Eqs.(A.24)  in  Appendix  A.  We  will  restrict  our  attention  to  reference  the  principal  coordinate 
system  and  redefine  the  subscripts  1,2  and  3  in  this  set  to  refer  to  the  components  of  the  field 
vectors  parallel  to  the  principal  axes  x,  y  and  z  respectively.  If  we  consider  the  direction  of 
propagation  in  a  spherical  coordinate  system,  the  directions  cosines  may  be  written 
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ka 

=  -  =  sin  0  cos  6 

1*1 


=  -=r-  =  sin  0  sin  < 

1*1 


(4.8) 


s,  =  ~  =  cos  0 

1*1 

where  0  is  the  angle  between  the  propagation  vector  and  the  z  axis  and  4>  is  the  angle  between 
the  projection  of  the  propagation  vector  onto  the  x-y  plane  and  the  x  axis. 

Recall  that  for  the  non-trivial  solution,  the  determinant  of  the  matrix  for  this  set  must  be 
equal  to  zero.  When  multiplied  out,  it  reduces  to  a  form  which  is  quadratic  in  n2.  If  we 
normalize  the  first  term  and  consider  the  cases  of  the  fundamental  and  second  harmonic 
separately,  the  two  separate  solutions  take  on  the  forms 
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(4.14) 

The  indices  of  refraction  refer  to  their  values  along  the  respective  principal  axes  for  a  given 
frequency,  to  or  2(0,  and  temperature. 

The  solutions  to  Eq.(4.9)  and  Eq.(4.10)  are 

V2  (4.15) 
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where  i  =  1  or  2,  refers  to  the  solution  for  the  case  of  the  square  root  term  being  positive  or 
negative  respectively  (24:65). 

Note  for  all  classes  of  crystals  (biaxial,  positive  and  negative  uniaxial),  n ^  and 

n2<*,2  >  n2a.\-  For  Type  I  phase  matching,  n^2  =  n^,,  resulting  in 

y[2  =  >/ 2  (4-17) 

For  the  case  of  a  uniaxial  crystal,  the  ordinary  and  extraordinary  indices  of  refraction  at  the 
fundamental  relate  to  the  indices  in  the  solution  to  the  phase-matched  condition  by 
n*,a> =  ny,<a  -  no.a  and  ni  b>  =  n4  a .  This  is  true  for  2co  as  well,  and  by  inserting  Eqs.(4. 11-14) 
into  Eq.(4.17),  the  phase-matched  angle  for  AgGaSe2  which  satisfies  this  relation  is 
0  =  55.022°.  The  graphical  representation  of  this  solution  is  shown  in  Figure  23.  The 
phase  matching  angle  corresponds  to  the  intersection  of  the  ordinary  index  of  refraction  for  co 
with  the  extraordinary  index  for  2co.  Note  that  there  is  no  dependence  of  the  solution  on  $  for 
uniaxial  crystals  because  ni  at  =  nyM.  The  optimum  value  for  4>  is  determined  by  the  effective 
nonlinear  coefficient.  The  general  solution  to  the  Type  I  phase-matched  condition  in  negative 
uniaxial  crystals  reduces  to  Eq.(2.2l)  and  for  positive  uniaxial  crystals  to  Eq.(2.22). 
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which  may  be  satisfied  for  some  particular  direction  of  propagation  defined  by  0  and  4>. 


The  indices  of  refraction  to  the  fundamental  and  second  harmonic  are  generally  derived 
from  the  appropriate  Sellmeier  equations  and  material  dependent  constants  which  are  listed  in 
the  table  below.  The  nominal  indices  of  refraction  at  10.6  j xm  for  the  two  crystals  considered 
in  this  study  are  2.59  for  AgGaSe2  and  3.1 1  for  ZnGeP2. 


Table  2. 

Sellmeier  coefficients  for  AgGaSe2  and  the  temperature  dependent  ZnGeP2 
(Note :  n2  =  A  +  5/(1  -C/X2)  +  Dl(\-E/X2)) 

X  is  in  micrometers,  and  E  (the  square  of  the  fundamental 
lattice  vibration  frequency)  is  generally  a  constant, 
and  temperature  is  in  degrees  Centigrade. 


Crystal  &  E 
Value  (pm2) 

Temperature 

(°C) 

Polarization 

A 

B 

C 

D 

AgGaSe2 

25 

0 

3.9362 

2.9113 

0.1507 

1.7954 

(1600.00) 

e 

3.3132 

3.3616 

0.1459 

1.7677 

ZnGeP? 

400 

o 

4.6171 

5.4709 

0.1495 

1.4912 

(662.55) 

e 

4.6791 

5.6826 

0.1567 

1.4577 

200 

o 

4.5654 

5.3382 

0.1419 

1.4913 

e 

4.6732 

5.4842 

0.1499 

1.4581 

100 

0 

4.5209 

5.2917 

0.1376 

1.4911 

e 

4.6559 

5.4001 

0.1460 

1.4580 

-100 

0 

4.3761 

5.2540 

0.1275 

1.4903 

e 

4.5801 

5.2747 

0.1368 

1.4576 

Adapted  from  Kildal  (16:317)  and  Bhar  (3:839) 

i  - ; - 

For  the  case  of  ZnGeP2,  whose  indices  of  refraction  are  given  as  functions  of  the  tempera¬ 
ture,  the  phase-matched  angle  is  found  in  the  same  manner.  The  only  difference  lies  in  that 
this  is  a  positive  uniaxial  cry  stal  and  the  solution  will  be  the  intersection  of  extraordinary  index 
of  refraction  at  the  fundamental  and  the  ordinary  at  the  second  harmonic.  Knowing  that  the 
temperature  required  for  non-critical  phase-matching  at  10.6pm  is  near  140°C  (3:541),  we  can 
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Figure  23.  Graphical  solution  to  the  phase-matching  angle  in  AgGaSe7.  Note  the  intersection 
of  the  extraordinary  index  of  refraction  at  the  second  harmonic  frequency  with  the  ordinary 
index  of  the  fundamental.  Note,  the  lines  are  staggered  due  to  the  limited  accuracy  of  the  pub¬ 
lished  indices  of  refraction. 

fit  a  second  order  polynomial  to  the  Sellmeier  coefficients  as  cited  in  Table  3  for  the  three 
higher  temperatures.  Table  4  shows  the  coefficients  of  these  fits  over  the  range  from  100°C  to 
400°C. 

These  temperature  dependent  coefficients  can  be  substituted  into  the  same  form  of  the 
Sellmeier  equation  cited  with  Table  3  leading  to  a  temperature  dependent  phase  matching 
angle.  The  results  of  this  analysis  are  shown  in  Figure  24.  The  analysis  says  the  material  is 
non-critically  phase  matched  over  a  large  temperature  range,  from  130°C  to  190°C.  The  dis¬ 
continuity  in  the  solution  can  be  attributed  to  the  limited  accuracy  of  the  computational  routine 
employed  to  find  the  intersections  of  the  temperature  dependent  indices  of  refraction. 
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Table  3. 


Constants  of  the  second  order  polynomials 
describing  the  temperature  dependence  of  the 
Sellmeier  coefficients  ioxZnGeP2  over  the  range  from  100°C  to  400°C 

Note:  (  Coefficient(7")  =  XT2  +  HT  +  C  ) 

T,  temperature,  is  in  degrees  Centigrade. 


Coefficient 

Polarization 

3 

C 

A(T) 

0 

-6.2167  x  10'7 

6.315  x  104 

4.46397 

e 

-4.7833  x  10 7 

3.165  x  10^ 

4.62903 

B  (T) 

0 

6.6167  x  10'7 

2.665  x  10^ 

5.25843 

e 

5.0333  x  107 

6.900  x  10^ 

5.32607 

C  (T) 

o 

-1.6667  x  10" 

4.800  x  10 5 

0.132967 

e 

-1.6667  x  10 8 

4.400  x  10  s 

0.141767 

D  (T) 

o 

1.4917  x  10 6 

-7.455  x  10-4 

1.55073 

e 

1.0000  x  10 8 

4.000  x  10  6 

1.45770 

E(T) 

o 

0.0 

662.55 

e 

0.0 

662.55 

Temperature  | Degrees  Centigrade] 


Figure  24.  Phase  matched  angle  as  a  function  of  temperature  for  ZnGeP2.  Note  that  the 
phase-matched  angle  is  nearly  7t/2  over  a  relatively  wide  angular  range.  This  is  the  case  of 
non-critical  phase  matching  in  this  positive  uniaxial  crystal. 
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These  solutions  for  both  types  of  crystals  tell  us  the  angle  of  propagation  required  to  mini¬ 
mize  the  phase  mismatch  M  and  thus  maximize  the  sine 2  term.  There  is  an  angular  depen¬ 
dence  included  in  the  effective  nonlinear  coefficient  which  has  yet  to  be  determined. 

Returning  to  the  set  of  simultaneous  equations  (Eqs.(2.24))  from  which  the  phase-matching 
conditions  were  derived,  we  are  now  interested  in  the  electric  field  components  of  the  funda¬ 
mental  and  second  harmonic  waves.  We  can  express  these  components  in  terms  of  the  cosines 
of  the  angles,  a,  p,  and  y,  with  respect  to  each  of  the  principal  axes,  x,  y  and  z  respectively,  as 

(cos  amo>1,  cos  PMOil  cos  ymatii)  (4>1 8) 

where  m  =  1  and  2  indicate  fundamental  and  second  harmonic,  and  i  =  1  or  2,  indicate  the  two 


possible  values  of  the  refractive  indices  (24:56).  For  second  harmonic  generation,  two  of  these 
field  components  will  be  the  same.  Given  a  direction  of  propagation  in  the  principal  coordi¬ 
nate  system,  the  direction  cosines  may  be  found  for  this  set  of  equations  under  the  restriction 

(cosam{DI)2  +  (cospmM1)2  +  (cosymioi)2  =  1  (419) 

Recall  Eq.(2.3),  the  second  order  polarization,  which  may  be  written  in  terms  of  these 


direction  cosines  and  the  two  applied  electric  fields  as  (24:66) 


P*°  =  a,dllka]ai  E /  ((a)Ei  (to) 


For  Type  I  phase  matching  we  have 


a,ak 


and  for  the  42m  class  crystals. 


cos  a, 
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cos2  y^2 
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0  0  0  0 


(4.20; 
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(4.23) 


a,  =  (coso^lCos^lCosy^,) 

Completing  the  matrix  multiplications  indicated,  we  find 

4/r  =  2d14cosa2£0ilcosP<n>2cosYa,2 

+  2d14  cos  ,  cos  a^2cos  7^2 

+  2d*  cos  72a,  1  COS  2  cos  p^2 

For  AgGaSe2,  a  negative  uniaxial  crystal,  the  solution  to  the  phase  matching  conditions 

requires  the  fundamental  to  be  perpendicular  to  the  optic  axis  (an  ordinary  wave).  This  makes 

7^2  =  tU2  and  all  terms  containing  cos  7^2  go  to  zero.  Only  the  third  term  in  Eq.(4.24) 

survives,  which  is  dependent  upon  d36.  We  can  relate  the  direction  cosines  of  the  electric  field 

to  the  angles  describing  the  propagation  direction  by  the  following 


Y2u>.i  ~ 


a, 


,  2  =  --<j> 

40. 2  2  Y 


(4.25) 


under  the  restriction  that  the  electric  field  of  the  fundamental  lies  in  the  x-y  plane,  d.^then 
becomes 


dcir 


(=-el 

|cos 

{2  TJ 

cos(tc  -  4>) 


(4.26) 


which  can  be  related  by  trigonometric  identities  to  the  more  familiar  form  quoted  for  negative 
uniaxial  crystals  of  this  class  (25:64) 

dtn  =  -d*sin0sin2(}>  (4.27) 

where  *36  =  43xl0'12  pmlV  for  AgGaSe2  (1  1:788).  This  is  a  relatively  insensitive 


function  of  4>  near  <J>  =  nJ4. 

Similar  arguments  may  be  applied  to  the  case  for  ZnGeP2,  which  is  a  positive  uniaxial 
crystal  and  the  electric  field  of  the  fundamental  is  restricted  to  lie  perpendicular  to  the  x-y 
plane.  Only  the  first  term  of  Eq.(4.24)  survives  and  it  may  be  related  to  the  angles  describing 
the  propagation  direction  in  a  similar  manner  to  obtain 

dtn  ~  d|4  sin  20  cos  2<J>  (4.28) 
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where  d,4  =  111  x  10‘12  pmlV  for  ZnGeP2  (4:302). 


An  immediate  source  of  concern  with  the  solution  for  ZnGeP2  lies  in  that  at  relatively 
moderate  temperatures,  less  that  200°C,  the  phase  matching  angle  is  approximately  n/2,  which 
makes  dfS\t ry  small.  In  fact,  in  order  to  take  advantage  of  the  relatively  high  nonlinear 
coefficient  dl4,  we  would  like  to  operate  at  an  angle  close  to  6  =  nJ4  and  0  =  0.  If  we  note  the 
orientation  of  the  principal  axes  in  Figure  1 1,  it  will  very  difficult  to  achieve  angles  approach  ¬ 
ing  these  values  in  the  actual  experiment.  The  specifics  of  these  limitations  will  be  discussed 
in  the  next  chapter.  Plots  of  the  d#  for  both  of  the  crystals  are  shown  in  Figures  25  and  26  as  a 
function  of  0  and  0.  Note  that  Type  II  phase  matching  is  not  possible  with  either  of  these 
crystals. 


Figure  25.  First  quadrant  of  the  dj  surface 
for  AgGaSe2.  The  plot  ranges  from 
0  $  0  S  itandO  S  $  ^  1^2, 
where  the  peak  value  is  d*  =  (43pm/V)z. 


Figure  26.  First  quadrant  of  the  surface 
for  ZnGeP2.  The  plot  ranges  from 
0  £  0  £  vJ2  and  x  to  0  £  0  <;  nJ2 
,  where  the  peak  value  is  dl,  =  (lllpm/V)2. 
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Bringing  the  developments  of  the  effective  nonlinear  coefficients  and  the  angular  sensitivi¬ 
ties  together,  the  first  order  approximation  to  the  conversion  efficiency,  Eq.(4.1),  may  be  writ¬ 
ten  in  terms  of  the  direction  of  propagation  in  the  principal  coordinate  system,  the  physical 
dimensions  of  the  crystal  and  average  intensities  (CW  input).  A  summary  of  the  calculated 
material  parameters  which  will  be  used  in  the  following  sections  is  shown  in  Table  5. 

Table  4. 


Calculated  Material  Parameters  for  AgGaSe2  and  ZnGeP2.  Values 
for  ZnGeP2  calculated  for  the  phase  matched  condition  at  21 1°C. 


AgGaSe2 

ZnGeP2 
(at  21 1°C) 

Ordinary  index  of  refraction 
at  10.6pm 
at  5.3pm 

2.5915 

2.6136 

3.1020 

3.1426 

Extraordinary  index  of  refraction 
at  lO.fipm 
at  5.3pm 

2.5585 

2.5809 

3.1430 

3.1833 

Phase  matching  angle,  0„ 

55.02° 

85.16° 

Angular  acceptance,  A0 
(for  a  29mm  AgGaSe2  crystal 
and  a  1 1mm  ZnGeP2  crystal) 

0.304° 

0.678° 

Angular  sensitivity,  pe 
(length  of  the  crystals,  L, 
not  included,  note  Eqs.(2.24-25)) 

-18089.8  rad/m 

4166.2  rad/m 

Birefringent  walkoff  angle,  p 

0.67° 

-0.13° 

Eq.(4.1)  will  be  valid  for  relatively  low  conversion  efficiencies  in  the  non-focussed  case 
(plane  wave).  The  solution  is  misleading  in  that  n  prompts  the  use  of  longer  crystals  for  higher 
efficiencies.  In  the  next  section  it  will  be  shown  that  as  the  power  in  the  second  harmonic 
grows,  it  gives  rise  to  a  dipole  oscillation  at  the  fundamental  frequency  which  reradiates  in  the 
phase-matched  direction. 
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Eq.(4.1)  provides  an  initial  estimate  of  the  overall  conversion  efficiency,  but  does  not 
include  several  important  parameters  of  interest.  If  we  retain  z  dependence  of  the  fundamental 
field  in  the  solution,  we  arc  able  to  track  the  amplitude  of  the  second  harmonic  as  a  function  of 
distance  as  well  as  the  depletion  of  the  fundamental.  The  development  in  Yariv  (23:528) 
yields  the  form 


>(2w. 


p(a) 


=  tanh 


V  2co, 


£,(0)kz 


(4.29) 


where  £,(0)  is  the  amplitude  of  the  fundamental  at  2  =  0,  and  tc  is  defined  by 


K 


*123^" 


e0 


(4.30) 


In  the  limit  that  the  argument  of  the  hyperbolic  tangent  approaches  infinity,  the  solution 
approaches  1/2.  Conservation  of  energy  requires  that  no  more  than  one  half  as  many  photons 
can  be  produced  at  the  second  harmonic  than  those  input  at  the  fundamental.  It  concludes  that 
with  longer  crystals,  the  closer  the  final  product  will  come  to  this  value. 

In  a  more  exhaustive  development,  the  divergence  of  the  pump  beam,  its  degradation 
through  the  crystal  and  the  angular  acceptance  of  the  material  play  major  roles  in  accurately 
predicting  the  efficiency  and  limiting  values. 

If  we  considered  the  case  of  a  single  pulse  at  the  frequency  of  the  fundamental  propagating 
through  the  crystal,  it  gives  rise  to  a  dipole  oscillation  at  the  second  harmonic  frequency.  This 
oscillation  radiates  a  field  whose  amplitude  grows  in  the  phase  matched  direction  as  function 
of  distance.  Given  that  the  energy  coupling  of  the  fundamental  pulse  to  the  dipole  radiators  is 
proportional  to  the  instantaneous  power,  the  temporal  peak  of  the  pulse  will  be  degraded  as  it 
gives  up  its  energy  to  the  crystal.  The  spatial  peak  will  be  degraded  as  well,  but  we  will 
assume  that  the  incident  monochromatic  plane  wave  possesses  a  constant  spatial  intensity 
profile.  If  we  disregard  the  local  variations  in  the  temporal  shape  of  the  pulse  due  to  this 
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factor,  we  can  expect  the  tariff  argument  in  Eq.(4.29)  to  describe  the  limit  of  this  coupling. 

The  exact  solution  to  the  coupled  wave  equations  (Eqs.(2.12)),  however,  describes  the 
secondary  reaction  of  the  dipole  oscillation  in  the  media.  The  field  at  the  second  harmonic  will 
promote  an  oscillation  at  the  fundamental  frequency  which  will  in  turn  cause  reradiation  at  this 
wavelength.  This  rise  and  fall  is  caused  by  the  competing  wave  interaction  of  the  two  fields 
can  be  described  by  a  general  class  of  solutions  known  as  solitons  (solitary  wave  solutions). 
This  back  conversion  to  the  fundamental  frequency  brings  about  a  limitation  on  the  useful 
length  of  material  in  producing  the  second  harmonic. 

For  Type  I  phase  matching,  the  equations  describing  the  interaction  of  the  first  and  second 
harmonics  for  a  single  pulse  with  an  input  intensity,  I,  and  a  crystal  thickness  L  have  been 
solved  (1:1925)  to  give  the  form  (13:576), 

Jl  .  1/2  I  (4.31) 

r|  =  tanhl  -tanh"'(5rt[2Tii/2,  1  +  82/4ti0]) 

where  the  nonlinear  drive,  ri0,  and  die  dephasing  terms,  8,  are  given  by 

T|0  =  C2IL2  (4.32) 

8  =  (1/2)A*L  .  (4.33) 

The  coupling,  C,  given  in  units  of  GW1'2  is  shown  below,  and  sn  in  Eq.(4.31)  is  the  Jacobi 
elliptic  function. 

C  =  5.46  d.jji  pm/V)  /  X1(n/n)(fl1rc2rtJ)''2  (4.34) 

Depletion  of  the  input  beam  due  to  conversion  to  the  second  harmonic  is  included  in  this 

development,  but  losses  due  to  absorption  have  been  omitted.  In  the  limit  of  a  small  nonlinear 
drive  term,  ti0,  this  function  reduces  to  the  sine2  dependence  of  Eq.(4.1)  and  is  a  direct  measure 
of  the  phase-matched  conversion  efficiency  (13:576).  For  larger  drives,  the  fundamental  pulse 
is  more  strongly  depleted,  and  the  point  where  T)0  =  1  can  be  considered  as  the  point  where 
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saturation  begins.  This  non-periodic  function  shows  the  strong  dependence  of  this  process  on 
the  beam  parameters  used.  The  contour  plot  shown  in  Figure  27  is  very  instructive  in  tracking 
the  predicted  conversion  efficiency  for  different  experimental  arrangements. 


Dephasing 
8  = 


Drive  %  =  C2L2/ 

Figure  27.  Contour  plot  of  conversion  efficiency  as  a  function  of  the  dimensionless  drive  and 
dephasing  terms  derived  using  the  Jacobi  elliptic.  There  are  nine  contours  corresponding  to 
efficiencies  from  0.1  to  0.9  over  the  range  of  the  plot. 

The  Jacobi  elliptic  function  in  Eq.(4.31)  may  be  written  explicitly  as  a  series  expansion 
(15:575)  of  the  nome  q  =  cxpi-nK’/K)  and  the  argument  v  -  nu/(2K ) 
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sn(u  |  m) 


(4.35) 


sin(2n  +  l)t> 


The  variables  K  and  K'  arc  functions  of  the  real  parameter  m  and  may  thus  be  written  without 
the  complex  notation  in  K'  as 


In  general,  the  m  in  Eq.(4.35)  must  be  less  than  one  in  order  for  the  series  to  converge.  Note 
the  specific  forms  of  the  parameters  of  the  Jacobi  elliptic  in  Eqs.(4.32-33),  where  1  +  52/4ri0 
will  always  be  greater  than  one.  By  Jacobi’s  real  transformation  (15:573) 

sn(u\m)  =  pir2sn(t/|p)  where  m>0,  v=umv 2  (4.37) 

Eq.(4.31)  may  be  rewritten  in  the  form 

q  =  tanhf -tanh-(  2^(1  +  52/4q0) 

[2  (Vl+62/4ri0  L 

where  the  nonlinear  drive  and  dephasing  terms  are  defined  above. 

Figure  27  describes  the  dependences  of  the  efficiency  on  the  crystal  length,  peak  intensity 
(of  the  temporal  and  spatial  profile)  and  the  phase  mismatch.  If  we  consider  the  case  of  second 
harmonic  generation  at  the  phase  matched  angle,  we  will  still  have  finite  contributions  to  the 
dephasing  from  the  far-field  divergence  half  angle  of  the  beam  (in  air),  50, 

Ak  =  p08©  +  p"602  <4-39) 

For  critically  phase-matched  materials,  such  as  AgGaSe2  whose  the  phase  matching  angle 

0  *  k/2,  the  first  term  dominates  and  A k  is  linear  in  50.  For  non-critically  phase  matched 
materials,  ZnGeP:  below  200°C,  the  first  term  goes  to  zero  and  the  dephasing  is  quadratic  in 
divergence. 

In  Figure  27,  one  may  be  misled  to  conclude  that  higher  intensity  levels  will  lead  to  better 
conversion  efficiencies.  Telescoping  the  beam  down  would  provide  this,  but  the  increase  is 


accompanied  by  a  corresponding  increase  in  the  divergence.  As  an  estimate  of  the  efficiencies 
expected  with  the  focussing  lenses  and  intensity  levels  proposed,  note  the  relations  between  the 
beam  and  crystal  dimensions  in  the  nonlinear  drive  term 

Tio  =  C\lI  =  C2(L1iarea)PFtti  <4-40) 

where  the  area  may  be  approximated  by  Ku?  at  the  center  of  the  crystal,  wmay  be  interpreted 

as  the  minimum  spot  size  in  the  focussed  case.  Doubling  the  dimensions  of  the  beam  and 

crystal  has  no  affect  on  the  drive  (13:577).  The  peak  intensity  must,  however,  be  chosen  to 

avoid  the  damage  thresholds  of  the  material. 

The  dephasing  term  may  be  expressed  as 

1  1  (4.41) 

6  =  ±AJcL  = 

where  Q  is  the  beam  quality  and  is  the  ratio  of  the  far-field  divergence  half  angle  to  the 
diffraction  limited  divergence  (both  taken  inside  the  crystal) 

50  58  (4.42) 

5q  X/(Kn  i a) 


By  eliminating  the  crystal  thickness  from  Eqs.(4.33-34),  the  condition  for  high  efficiency  (in 
excess  of  50%),  independent  of  the  crystal  length,  is  given  by  (13:577) 

P?tak  *  (P ^JCf  (n</462)  •  Q1  <4-43> 

In  AgGaSe2  we  can  estimate  this  value  with  a  few  minor  computations.  The  only 

assumption  with  regard  to  the  spatial  intensity  (Gaussian)  pattern  of  the  beam  will  be  in 
calculating  the  peak  intensity  of  the  pulse.  If  we  consider  the  case  of  focussed  second 
harmonic  generation  using  a  lens  of  focal  length  12.5  cm,  with  a  minimum  spot  size  (peak 
intensity)  as  given  in  Figure  18,  we  can  find  the  far-field  divergence  as  one  half  of  the 
numerical  (or  analytic)  derivative  of  Eq.(3.8). 


Q 


56  _  10.39mrad/2 

X/nnu)  ~  10.6p/n/(7t(2  5915)(125.3p/n)) 


0.500  rad 


(4.44) 
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which  says  that  the  peak  power  in  the  pulse  must  be  greater  than  90.7kW  in  order  to  achieve 

high  conversion  efficiency  in  this  configuration.  In  our  case,  the  peak  powers  range  for 
4kW  to  8kW  and  we  will  not  observe  conversion  efficiencies  quite  this  high. 

In  order  to  maximize  the  efficiency,  we  wish  to  operate  with  as  small  of  a  dephasing  term  as 


54 


I 

I 

I 


possible  while  maintaining  a  value  for  the  nonlinear  drive  term  as  near  to  optimum  as  possible, 
but  below  the  damage  threshold.  In  actuality,  this  configuration  for  AgGaSe2,  given  the 
available  power  and  the  12.5  cm  focal  length  lens,  is  near  optimum. 

Improvements  to  Plane  Wave  Solutions  -  Radiai  Intensity  Profile 

In  both  of  the  solutions  above,  the  spatial  intensity  profile  of  the  CW  beam  or  pulse  was 
assumed  to  be  constant  In  the  initial  derivation  of  the  efficiency  calculations,  we  are  able  to 
include  this  variation  in  a  rather  simple  fashion  by  rewriting  Eq.(2.15)  as  (23:528) 

(4.50) 

.  Vl__? _ 1 

ltjk> 


I E 


where 


K  .=  Ei  =  £'»  =  Eee 
If  we  allow  the  power  in  the  fundamental  to  be  described  by 


(W 


>(<*) 


-  iVT  JI 


E^Ydxdy 


.  ./I>i 

VMo  °4 


Kill 


(4.51) 


(4.52) 


Mo  J  J  cron  uctioo  V  Mo 

■ufo  is  the  minimum  beam  waist  for  the  focussed  case.  Combining  Eq.(4.50)  and  Eq.(4.52) 


p( 2“>) 

1  =  Jm  =  *1 


'MoYV(4)*L2  pw  sin2;A kL 
\eoJ  n 3  7t4  (j/ttL? 

which  is  identical  to  Eq.(4. 1)  except  that  it  was  derived  for  a  Gaussian  beam  input  with  z0»L 


(4.53) 


and  is  accounted  for  in  the  factor  area  =  Jn (13:530).  z0  is  the  Rayleigh  range  defined  by 
z0  •-  n ufafk.  This  solution  is  applicable  to  focussed  second  harmonic  generation  as  suggested 
in  Figure  28. 

The  efficiency  will  increase  with  the  length  of  the  crystal  up  to  the  point  where  z0  =  L ,  at 

which  point  the  divergence  of  the  beam  over  the  length  of  the  crystal  becomes  appreciatable. 
The  optimum  condition  exists  for  z0  =  LI 2  (22:217).  For  the  case  of  AgGaSe2,  with  a  lens  of 
focal  length  12.5  cm  placed  50  cm  from  the  laser 
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Figure  28.  Gaussian  beam  focussed  inside  a  nonlinear  optical  crystal. 


L  =  2[n(125.3pjn)2  (2.5915)]/I0.6jxm  =  2.412  cm  (4-54) 

which  suggests  that  a  slightly  weaker  focussing  case  would  bring  us  closer  to  the  optimum 

efficiency.  This  2grees  with  the  exact  analysis  in  the  sense  that  with  a  large  focussed  spot  size, 
we  would  have  a  lower  divergence  (and  dephasing  term)  leading  to  higher  overall  efficiencies. 
Eq.(4.53)  is  applicable  for  relatively  short  crystal  lengths  at  moderate  intensity  levels  in  the 
CW  harmonic  generation. 

For  pulsed  second  harmonic  generation,  the  exact  solution  (Eq.  (4.38))  does  not  easily  lend 
itself  easily  to  folding  in  the  radial  intensity  distribution  of  the  beam.  A  rough  approximation 
can  be  achieved  by  segmenting  the  radial  pattern  into  discrete  rings  at  different  intensity  levels, 
weighting  their  contribution  to  the  total  efficiency  by  the  amount  of  energy  in  each  and 
summing  the  components.  We  can  describe  the  Gaussian  radial  intensity  distribution  by 


/*V)  -  w"'-*  <455) 

We  can  then  sum  a  series  of  discrete  evaluations  of  the  efficiency  using  the  exact  solution  and 


weight  each  by  the  integral  cf  the  total  power  in  A r  by 
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The  steps,  A r,  and  final  value  for  the  summation  in  r  are  determined  to  yield  an  accurate  result 
in  a  reasonable  amount  of  computation  time. 

Numerical  Approximation  to  the  Eocussed  Beam  Sfllution 


Implicit  in  the  previous  derivations  was  the  assumption  that  the  focussed  spots  at  the  funda¬ 
mental  and  second  harmonic  frequencies  overlap  over  the  entire  length  of  the  crystal.  Recall 
from  the  discussion  of  the  anisotropic  properties  of  these  crystals,  one  of  the  two  orthogonally 
polarized  beams  may  walkoff  from  the  other.  This  deviation  is  accounted  for  in  the  solution 
presented  below  which  numerically  tracks  the  spatial  development  of  the  field  at  2wover  the 
length  of  the  crystal.  It  takes  into  account  diffraction,  double  refraction  and  absorption  exactly 


within  the  limit  of  the  paraxial  approximation. 

Following  the  development  in  Boyd  (5:3600;  and  Eckardt  (12:5),  we  define  the  conversion 
efficiency  in  Eq.(4.57)  in  terms  of  the  integral  expression  shown  in  Eq.(4.61).  This  theory  was 
originally  developed  for  CW  second  harmonic  generation  (5:3600)  and  has  been  modified  for 


the  pulsed  case  with  two  different  approximations.  Assuming  a  Gaussian  shaped  pulse  in  '  me, 

I — 

a  factor  of  l/\2  was  introduced  to  derive  an  effective  average  power  over  the  length  of  the 


exposure.  A  second  factor,  RSaJ,  was  used  to  account  for  saturation  of  the  medium  (10:925). 


1 


J/?s^'Lfc/V0)Mo,(U,U) 


(4.57) 


where 


O),  Fundamental  frequency,  radians/sec 

d ^  Effective  nonlinear  coefficient,  m/V 

Xo  Free  space  wavelength,  m 

L  Length  of  the  crystal,  m 

n  Nominal  index  of  refraction 

The  fundamental  power  at  the  peak  of  the  temporal  profile  is  defined  by, 


P,(  0) 


P 

*  Average 


(4.58) 
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where  x  is  the  He  point  of  the  pulse  (characteristic  time)  and  v,v  is  the  pulse  repetition 
frequency.  Saturation  is  accounted  by 

„  -  ^ 

KSl  — 


(4.59) 


H**+l 

This  factor  is  ge  erally  included  once  the  rest  of  the  computations  have  been  completed  and 
leads  to  a  favorable  approximation  of  the  pulsed  efficiency  up  to  about  50%.  We  go  on  to 
define 


h{c,  (3,  k,  fa)  -  (^/^^(o.p.K,^) 

_2_ 

Vre- 


f(a,p,K,^p)  =  |//(o  +  4(ir,K4,p)|2exp(-452kii 


We  can  define  o’  =  o  +  4(U  such  that  (5:3602) 

I  //(o’,  K,^,p)  |2  =  ^  (/?  (o’,  tc, %, p)2  +  /  (O’,  >c, |i)2] 
The  real  and  imaginary  pans  (17:379)  may  be  written  as 


(4.60) 
(4.6’ ) 


(4.62) 


/?(a’,tc,£,p) 

J 

U'*M) 

dx  „ 

- ;(COSO’T  +  X  sin  o’T)<T 

1  +r 

(4.63) 

/(o’.tc.^.M) 

-  J 

+|i) 

-V'-M) 

dx  „ 

- r(sino’T-Tcoso’x)<? 

1  +r 

(4.64) 

The  parameters  are  defined  as; 

o  = 

\b*k 

,  phase  mismatch 

(4.65) 

P  - 

Bt;'a 

,  ratio  of  p/5o 

(4.66) 

K  = 

,  absorption 

(4.67) 

4  = 

Lib 

,  focussing  parameter 

(4.68) 

b  = 

,  confocal  parameter 

(4.69) 
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H  =  (L-2f)/L 


,  focai  position 


(4.70) 

(4.71) 


1  1 
a  =  a, --a,  and  a  =  a,  +-o^ 

a,  and  CC2  are  the  absorptivities  at  the  fundamental  and  second  harmonic  respectively  and 

14,  is  the  focussed  minimum  spot  size.  The  integral  generally  yields  a  single  peak  for  a  given 
value  of  the  double  refraction  parameter, 

B  =  p(L/k,)1'2/  2  (4-72) 

where  p  is  the  walkoff  angle  in  radians,  and  k,  is  the  propagation  constant  inside  the  media. 

This  theory  is  written  for  the  case  of  negative  uniaxial  crystals  and  may  be  extended  to  the 
positive  uniaxial  case.  These  computations  may  be  completed  using  Romberg  integration 
routines  for  both  of  the  expressions  above  (See  Appendix  B). 

For  a  specified  crystal  type  and  physical  dimensions,  this  expression  yields  a  single 
maximum  for  the  optimum  focussing  conditions.  Our  AgGaSe2  crystal,  which  is  2.9  cm  in 
length,  has  a  walkoff  angle  of  1 1.69  mrad  and  B  =  1.234.  The  corresponding  optimum 
focussing  parameter,  extracted  from  Figure  4  of  Boyd,  is  \  =  1.775  which  leads  to  an 
optimum  focussed  spot  size  of 

rrr 

W°  \  2TlnfU 

This  value  is  smaller  than  that  predicted  earlier  by  the  confocal  focussing  condition,  and  is 
smaller  than  the  set-up  involving  the  12.5  cm  focal  length  lens  placed  50  cm  from  the  laser. 

We  can  easily  reach  this  value  by  moving  the  lens  further  from  the  C02  source  (z2  =  65cm,  and 
this  was  accomplished  in  the  latter  portions  of  the  experiment. 

For  the  1.1  cm  sample  of  ZnGeP}  heated  to  21 1°C,  at  the  critically  phase  matched  angle  of 
35.i  p  =  0.13°  and  B  =  0.155,  leads  to  an  opti  um  focussing  parameter  of  ^  =  2.82.  This 
yields  an  optimum  spot  size  of  74.8jim  which  can  be  achieved  with  the  available  experimental 
setup. 


1  0.029m  10.6pm  , 

- — -  =  103.1  i im 

2  ti(2. 5915)(1.775)  w 


(4.73) 
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Heating  Estimates  and  Damage  Thresholds 


As  a  precursor  to  the  actual  experiments,  a  first  order  study  of  the  heating  in  the  crystal 
caused  by  the  absorption  of  the  fundamental  as  a  function  of  the  average  input  power  was  con¬ 
ducted.  The  temperature  rise  can  be  estimated  by  the  amount  of  energy  absorbed  at  a  particular 
point  over  a  specified  time  duration  via  the  heat  capacity,  C,  and  absorptivity  of  the  material  of 
mass  m  by 

AQ  (4.74) 

AT  =  — 
mC 


AQ  is  the  amount  of  energy  absorbed  and  we  will  relate  it  to  the  power  by  considering  a  1  sec 

exposure  to  the  beam.  We  use  the  time  averaged  energy,  or  power,  because  of  the  very  high 
repetition  rates  of  the  C02  source.  If  we  consider  the  absorptivity  to  be  constant,  we  can 
integrate  the  total  energy  absorbed  in  1  mm  thick  slices  along  the  length  of  the  crystal  by 

(4.75) 

AQ(l)  =  PAhAi  J  [e\p(-a/ )  -  exp(-oc(/  +  A//10))]  dl 


where  /  can  range  from  0  to  L  =  29  mm,  AI  =  0. 1  mm  ,  At  =  1  sec  and  the  absorptivity  is 


defined  per  millimeter.  If  wc  were  to  integrate  this  expression  over  the  length  of  the  crystal, 
we  would  get  the  total  energy  absorbed. 

To  fold  in  the  radial  dependence  of  the  energy  distribution,  we  considered  the  energy 
absorbed  by  a  particular  segment  of  the  crystal  and  write  in  the  from  of  a  radial  Gaussian 
distribution  as 


Q(r,D 


Q(  D  -rW 

nu? 


(4.76) 


Finally  the  temperature  profile  at  a  given  point  within  the  crystal  as  a  function  of  radius  may  be 
written 


AT{r,l) 


Qd ) 

nu?(mC) 


(4.77) 
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The  heat  capacity  is  generally  given  in  J/mole/’C  and  the  mass,  m,  can  be  found  from  the 
density  of  the  material  and  the  given  area  of  integration  (consider  a  long  cylinder  where  the 
cross-sectional  area  cancels  out  in  the  computation).  For  the  case  of  the  AgGaSe2  crystal,  with 
C  =  104J/mole/°C  and  p  =  5.7g/cm3 ,  the  temperature  profile  for  three  spot  sizes,  w,  at  the 
center  of  the  crystal  with  a  =  0.0027  /mm  are  shown  in  Figure  29. 


Figure  29.  Temperature  profiles  for  a  one  second  exposure  at  30  kHz,  nominal  average  power 
of  6.0  Watts  at  the  center  of  the  AgGaSe2  crystal.  Note  that  conduction  of  the  heat  was  not 
considered  in  this  estimate. 

These  profiles  appear  very  steep,  as  they  are  shown  for  the  center  0.2  mm  radius  of  the 
beam,  and  although  no  thermal  conductivity  was  taken  into  account,  we  can  expect  reasonably 
high  stress  within  the  material  when  exposed  to  these  kinds  of  power.  A  more  rigorous  explo¬ 
ration  would  relate  the  energy  absorbed  to  the  thermal  conductivity  and  expansion  coefficients 
to  deduce  the  stresses.  For  our  purposes,  this  gives  us  a  feel  for  limiting  the  exposure  durations 
at  high  powers  in  light  of  the  anticipated  temperature  build  up. 

From  the  variety  of  analytic  and  numerical  approximations  to  the  expected  conversion  effi¬ 
ciency,  we  can  now  estimate  the  optimum  experimental  parameters  under  which  to  test  and 
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characterise  the  materials.  The  heating  studies  only  give  us  a  crude  idea  of  the  stresses  on  the 
material  but  are  useful  in  comparing  to  known  damage  thresholds  of  the  two  crystals.  These 
will  be  reviewed  in  the  next  chapter  with  the  experimental  method. 
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With  the  review  of  the  theoretical  calculations  in  hand,  we  are  able  to  describe  the  exper¬ 
imental  methods  employed  in  determining  each  of  the  figures  of  merit.  The  various  depen¬ 
dences  of  the  conversion  efficiency  on  the  peak  and  average  powers,  direction  of  propagation 
and  focussing  were  further  explored. 

As  noted  in  the  derivation  of  the  effective  nonlinear  coefficients,  the  conversion  efficiency 
of  both  crystals  will  be  heavily  dependent  upon  the  direction  of  propagation  with  respect  to  the 
optic  axis,  0,  and  relatively  insensitive  to  the  angle  <j>.  The  critically  phase  matched  AgGaSe2 
will  have  a  very  small  angular  acceptance  and  demand  precise  movement  about  the  crystal’s 
central  axis  to  accurately  measure  the  conversion  efficiency’s  dependence  on  this  angle.  We 
anticipate  being  able  to  measure  SHG  in  the  non-focussed  and  focussed  cases  as  functions  of 
angle  of  incidence  with  respect  to  the  principal  axes,  intensity  levels  (position  of  the  crystal 
with  respect  to  the  lens)  and  dephasing  (different  lenses). 

The  experiments  with  ZnGeP2  will  prove  a  bit  more  involved  due  to  the  requirement  for 
temperature  tuning.  Physically  mounting  the  heating  oven  will  restrict  the  limits  of  angular 
deviation  from  perpendicular  to  the  face.  Recall,  this  is  a  non-critically  phase-matched  mate¬ 
rial  over  a  reasonably  large  temperature  range,  but  we  anticipate  extremely  low  conversion 
efficiencies  due  the  magnitude  of  d^near  0  -  71/ 2.  We  do  not  anticipate  observing  any  second 
harmonic  with  the  low  intensity  levels  of  the  non-focussed  case,  thus  these  experiments  will  be 
restricted  to  the  focussed  case  only.  Additionally,  in  order  to  boost  d &  we  must  elevate  the 
temperature  of  the  crystal  so  that  0„  <  nl 2.  The  major  disadvantage  this  produces  is  that  the 
material  becomes  critically  phase-matched  and  these  oblique  angles  of  incidence  will  prove 
more  difficult  to  align  with  the  oven  in  place. 
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The  initial  experiments  were  conducted  without  the  use  of  a  focussing  element  in  order  to 
compare  the  results  to  the  theoretical  calculations  involving  the  plane  wave  solutions  as  well  as 
to  avoid  any  unforeseen  heating  or  damage  mechanisms  that  may  have  been  encountered. 

From  the  calculated  spot  size  and  peak  intensity  values  for  the  unfocussed  laser  beam,  at 
these  relatively  low  intensity  values,  the  closer  we  operate  to  the  output  port,  the  higher  the 
conversion  efficiency.  We  must  also  consider  the  physical  size  of  the  crystal  face  in  order  to 
ensure  that  most  of  energy  traverses  the  crystal’s  length  without  being  clipped.  If  the  spot  size 
is  less  than  one  quarter  of  the  width  of  the  crystal  face,  2.25  mm,  approximately  98%  of  the 
energy  will  pass  through  the  crystal.  This  restricted  the  operating  distance  to  within  15  cm  of 
the  output  port  of  the  laser.  In  Figures  30  and  31,  the  measurements  of  the  unfocussed  conver¬ 
sion  efficiency  were  taken  at  this  point. 

An  electrically  operated  mechanical  shutter  was  incorporated  into  the  design  in  order  to 
limit  the  exposure  durations.  It  was  generally  repeatable  to  within  ±0.05sec  for  exposures 
longer  than  0.5  seconds.  The  timing  and  actuator  circuit  was  foot  activated  and  had  to  be  iso¬ 
lated  from  the  power  sources  within  the  room  in  order  to  reduce  the  EMI  produced.  In  later 
experiments,  the  exposure  duration  was  directly  measured  with  a  small  HeNe  laser  and 
photovoltaic  detector. 

In  the  non-focussed  experiments,  two  methods  were  used,  the  difference  being  that  the  first 
involved  a  series  of  energy  measurements  using  a  single  detector  and  the  second  involved 
power  measurements  using  two  detectors.  Initially,  we  were  not  able  to  acquire  two  energy 
meters  and  were  forced  to  conduct  the  energy  experiment  with  a  single  1 "  disk  calorimeter. 
Note  Figure  30,  where  the  energy  in  the  fundamental  was  measured  for  the  0.5  sec  exposure 
before  and  after  the  crystal  was  placed  on  the  stage.  The  energy  at  the  second  harmonic  was 
recorded  and  compared  to  these  fixed  values  to  derive  an  efficiency  as  a  function  of  the  angle 
6. 
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Figure  30.  Experimental  set-up  for  measuring  the  non-focussed  conversion  efficiency  as  a 
function  of  angle  of  incidence,  0,  using  energy  measurements  of  a  single  meter. 

The  second  set  of  measurements  at  this  same  distance  were  taken  in  terms  of  power.  Figure 
31  shows  the  second  calorimeter,  a  1"  diameter  Coherent  Radiation  laser  power  head  and 
meter,  was  used  to  measure  the  fundamental  reflected  from  the  10.6  pm  attenuator  (discussed 
later)  after  it  passed  through  the  crystal.  By  taking  into  account  the  absorptivity  of  the  crystal, 
this  can  be  used  to  relate  the  power  incident  to  the  second  harmonic  power  produced.  At  these 
very  low  efficiencies,  depletion  of  the  fundamental  was  be  negligible. 

Note  the  small  laser  used  for  alignment  purposes  in  all  of  these  experiments  was  a  1.0  mW 
HeNe  operating  at  547  nm.  Since  the  C02  beam  is  difficult  to  track  through  the  system  with¬ 
out  disturbing  the  optics,  the  visible  alignment  laser  was  an  essential  tool. 

After  the  beam  passed  through  the  crystal,  the  second  harmonic  and  fundamental  were 
separated  by  the  filter  stack  composed  of  a  30.6pm  reflective  interference  coating  (reflects 
99%  of  the  fundamental)  on  a  ZnSe  substrate,  followed  by  a  1/8"  thick  sapphire  plate.  The 
sapphire  plate  blocks  all  wavelengths  longer  than  approximately  6p.m.  The  transmission  spec¬ 
trum  of  this  filter  stack  was  measured  using  a  FTIR  and  is  shown  in  Figure  32.  In  recording 


65 


Figure  31.  Experimental  set-up  for  measuring  the  non-focussed  conversion  efficiency  as  a 
function  of  angle  of  incidence,  0  using  power  measurements  with  two  meters. 

the  background  and  sample  scan  with  the  FTIR,  a  polarizer  aligned  parallel  to  the  anticipated 
output  polarization  of  the  second  harmonic  (parallel  to  the  table  top)  was  placed  in  the  path  of 
the  beam.  This  helped  to  remove  the  ambiguities  in  the  transmission  introduced  by  the  use  of 
an  interference  filter. 

The  mounting  stage  provided  control  of  the  angular  orientation  0  to  within  ±3  min  and  <J> 

to  ±0.5  min.  An  absolute  measure  of  the  distance  between  the  laser  and  the  crystal  face  could 
be  made  with  an  accuracy  of  ±1  mm  while  relative  measurements  made  with  a  micrometer 
stage  were  ±0.005  mm.  Lateral  movement  of  the  crystal  was  included  to  compensate  for  the 
movement  of  the  crystal  face  .  i  it  was  slewed  about  its  central  axis.  As  a  reference  ooint,  note 
Figure  33  which  depicts  the  various  axes  and  directions  of  movement  over  which  we  had  con¬ 
trol.  They  are  denoted  by  letters  for  reference  in  Chapter  IV. 
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Figure  32.  Infrared  transmission  spectrum  of  the  pair  of  filters  used  to  block  the  fundamental 
wavelength  at  10.6  p.m  and  pass  the  second  harmonic  at  5.3  p.m  where  the  transmission  is 
49.8%. 


Figure  33.  Diagram  showing  the  axes  of  rotation  and  translation  for  the  AgGaSe2  crystal. 
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At  the  point  of  peak  conversion  efficiency,  the  temporal  waveforms  of  the  fundamental  and 
second  harmonic  were  recorded  using  the  high  speed  LN2  cooled  HgCdTe  detectors.  The 
detectors  were  simply  placed  where  the  calorimeter  heads  resided  and  at  these  moderate  con¬ 
version  efficiencies,  the  signal  strength  of  the  second  harmonic  was  at  the  appropriate  level  for 
detection.  It  is  interesting  to  note  the  resonant  cavity  detectors  cut  for  the  fundamental 
wavelength  were  very  sensitive  to  the  5.3  pm  radiation  as  well. 

The  procedure  for  taking  each  of  the  measurements  included  several  steps.  The  energy 
meters  tended  to  drift  with  the  thermal  fluctuations  in  the  room  caused  by  air  flows  especially 
when  taking  measurements  on  the  more  sensitive  ranges  (0.01  J  or  less).  They  required  contin¬ 
uous  observation  despite  their  isothermal  enclosures.  Once  they  had  settled  on  a  reliable  zero 
reading,  the  next  step  was  to  ensure  the  laser  was  running  in  a  stable  mode.  This  was  done  by  ^ 
observing  the  loop  error  meters  monitoring  the  piezoelectrically  controlled  mirrors  of  the  three 
oscillators.  Once  both  the  meters  and  laser  were  stabilized,  the  crystal,  in  the  proper  orienta¬ 
tion,  was  exposed  by  depressing  a  foot  switch  to  activate  the  mechanical  shutter.  The  energy 
in  the  fundamental  and  second  harmonic  were  then  recorded.  The  meters  were  reset  and  the 
zero  points  observed  to  ensure  no  fluctuations  had  taken  place  during  the  exposure. 

Focussed  efficiency  measurements  of  A%GaSez 

In  the  design  of  the  focussed  beam  experiments,  we  are  given  greater  flexibility  in  arranging 
the  components  since  the  constraints  on  the  placement  of  the  crystal  with  respect  to  the  laser 
have  been  lifted.  The  initial  bulk  of  the  experiments  were  conducted  with  the  focussing  ele¬ 
ment  50  cm  from  the  laser  since  at  this  distance,  the  spot  size  of  the  laser  is  just  less  than  1/4  of 
the  diameter  of  the  1"  diameter  ZnSe  positive  meniscus  focussing  lens.  Also,  the  calculated 
values  for  the  peak  and  average  intensities  would  be  below  published  damage  thresholds  for 
AgGaSe2.  In  previous  studies,  the  surface  damage  for  an  uncoated  sample  occurred  for  a  peak 
intensity  level  of  12  MW/cm2  (1  1:788)  in  a  single  pulsed  mode  and  14.7  MW/cm2  (12:10)  in  a 
repetitively  pulsed  format.  The  average  power  damage  thresholds  at  the  high  repetition  rates 
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was  given  as  19.4  kW/cm2  given  a  spot  size  a  the  crystal  face  of  164pm  (12: 10).  In  the  case  of 
the  12.5  cm  focal  length  lens,  with  a  minimum  focussed  spot  size  of  125.3  \im  (placed  50  cm 
from  the  laser),  the  average  power  seen  at  the  crystal  face  would  have  to  exceed  9.6  W  in  order 
to  reach  the  average  power  damage  threshold.  The  average  power  at  30kHz  seen  through  the 
focussing  lens  is  below  6  Watts. 

Energy  measurements  using  1  second  exposures  were  used  in  the  initial  investigations.  In 
order  to  investigate  any  transient  heating  effects  during  exposure,  a  copper-constantine  thermo¬ 
couple  was  incorporated  into  the  crystal  mount  and  the  temperature  was  recorded  during  vari¬ 
ous  exposures.  The  thermocouple  was  sensitive  to  iO-l^C  changes  in  temperature  and  there 
was  no  obse  vable  rise  during  exposures  ranging  up  to  5  seconds  in  duration.  For  this  reason, 
average  power  measurements  of  the  conversion  efficiency  were  conducted  as  well.  Here, 
exposures  up  to  20  seconds  were  taken  as  the  meters  at  {joints  I  and  II  (set  to  measure  power) 
settled  on  the  average  power  readings. 

Since  the  average  energy  from  the  laser  tended  to  fluctuate  slowly  over  time,  the  energy,  or 
power,  incident  on  the  crystal  was  measured  as  a  function  of  that  reflected  from  the  beam  split¬ 
ter  placed  near  the  shutter  (Note  Figure  34).  Without  the  crystal  in  place,  the  energy  detected 
at  point  II  was  calibrated  against  that  detected  at  point  I  over  a  given  exposure  duration.  The 
energy  could  be  related  to  the  average  power  using  this  exposure  time  and  a  similar  calibration 
between  the  two  meters  with  the  shutter  held  open.  This  calibration  was  conducted  with  each 
set  of  experiments  in  order  to  ensure  the  accuracy  of  the  results. 

The  physical  arrangement  of  the  experimental  apparatus  is  shown  in  Figure  34.  With  the 
additional  room  provided  by  placing  the  lens  far  from  the  laser,  we  were  able  to  incorporate 
several  other  elements  to  control  the  beam.  The  ZnSe  beamsplitters  were  used  to  vary  the  peak 
and  average  intensity  levels  while  the  chopper,  shown  near  the  mechanical  shutter,  allowed  for 
adjustment  of  the  average  intensity  without  affecting  the  peak  intensity  of  the  transmitted 
power. 
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Figure  34.  Experimental  arrangement  for  the  focussed  beam  conversion  efficiency  measure¬ 
ments  in  AgGaSe2. 

The  second  HeNe  laser  was  used  to  measure  the  exposure  duration.  When  the  shutter  was 
open,  the  signal  detected  by  the  silicon  photovoltaic  detector  was  viewed  on  an  oscilloscope 
and  the  FWHM  measured. 

The  HgCdTe  detectors  were  used  to  measure  the  relative  temporal  profiles  of  the  fundamen¬ 
tal  and  second  harmonic  at  various  intensity  levels  and  crystal  orientations.  The  sapphire  plate 
in  front  of  one  of  the  detectors  indicates  which  was  used  to  record  the  5.3  pm  signal. 

Recording  the  spatial  profiles  of  the  fundamental  and  second  harmonic  was  a  straight  for¬ 
ward  procedure  of  placing  the  pyroelectric  vidicon  camera  in  the  position  of  the  energy  meters 
and  adjusting  the  intensity  of  the  team  to  the  appropriate  level  for  a  non-saturated  response. 

The  procedure  for  exposure  was  the  same  as  that  for  the  non-focussed  case.  Variations  in 
the  chopper  frequency  and  duty  cycle  were  investigated  as  well. 

Focussed  Efficiency  Measurements  of  ZnGeP , 

In  the  case  of  Z.nGeF2 ,  the  temperature  controlled  oven  was  incorporated  into  the  design. 
From  the  calculations  of  d ^  the  direction  of  propagation  with  respect  to  the  optic  axis,  0, 
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would  be  the  crucial  determinant  of  the  conversion  efficiency.  Near  the  noncriticai  phase- 
matching  angle  and  temperature,  we  expect  very  poor  results  and  are  forced  to  go  to  higher 
temperatures.  Although  at  these  temperatures  the  material  is  critically  phase-matched,  we 
hoped  to  observe  some  output  as  the  temperature  passes  through  the  phase-matching  point  for  a 
given  orientation.  The  feedback  loop  for  the  oven  control  circuit  is  shown  in  Figure  35. 


Figure  35.  Experimental  arrangement  for  the  focussed  beam  conversion  efficiency  measure¬ 
ments  in  ZnGePj. 

The  hot  plate  is  operated  in  a  full  crVoff  manner  by  a  single  set  point  thermistor  gauge 
monitor  reading  the  temperature  with  a  copper-constantine  thermocouple  mounted  near  the 
crystal  in  a  copper  heat  sink.  During  initial  warm-up  there  is  a  fair  amount  of  overshoot  in 
temperature  because  of  the  thermal  mass  of  the  hot  plate  and  copper  block.  Once  settled,  the 
thermistor  is  able  to  maintain  a  set  temperature  to  with  ±5°C.  A  second  thermocouple 
mounted  on  the  opposite  side  of  the  crystal  monitors  the  temperature  at  this  point  with  an  out¬ 
put  connected  to  a  strip  chan  recorder. 

The  crystal  orientation  may  be  varied  in  0  (an  external  angle  in  this  case)  up  to 
approximately  15°  where  the  hot  plate  begins  to  obscure  the  beam  impinging  on  the  crystal. 
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The  internal  angle,  0,  at  this  point  is  8nMidr.Z0ntp  =  O.O&lrad  ,  or  5°  of  arc.  The  phase  match¬ 
ing  temperature  at  this  point  is  21 1°C  and  is  easily  attainable  with  this  hot  plate.  There  were  no 
windows  around  the  crystal  to  help  maintain  a  low  temperature  gradient,  so  the  overshoot  in 
the  temperature  control  may  have  been  to  our  advantage. 

Once  the  orientation  and  position  of  the  crystal  behind  the  lens  was  set  and  the  temperature 
was  near  the  critical  value  the  exposure  procedure  was  the  same  as  the  previous  cases. 

Damage,.. St  wflies.Qf  AgGaS.?i 

Damage  point  data  for  AgGaSe2  is  limited  in  the  open  literature  due  to  the  expense  and 
rar  f  the  material  at  this  point.  The  threshold  values  may  vary  widely  with  crystal  quality, 
surfar  coatings  :md  in  some  instances,  interplay  between  the  average  and  peak  power  damage 
mechanisms.  There  may  be  instances  where  it  is  i  fiuenced  by  long  term  fatigue  due  to 
repeated  exposures  ;tt  the  same  point.  Although  damage  occurred  with  the  original  sample  at 
two  points  in  the  efficiency  measurements,  small  witness  samples  were  provided  with  coated 
and  non-coated  surfaces  for  this  specific  purpose  of  determining  these  thresholds. 

The  published  threshold  limits  may  be  reached  or  exceeded  by  using  the  6.05  cm  focal 
length  Ge  lens  or  by  moving  the  12.5  cm  focal  length  ZnSe  lens  to  a  point  further  from  the 
laser.  For  each  of  the  samples  it  is  then  a  matter  of  translating  their  position  with  respect  to  the 
position  of  the  focus  until  damage  occurs.  The  threshold  can  be  characterized  in  terms  of  peak 
and  average  intensities,  spot  size  and  bum  spot  diameters. 

The  experimental  procedures  outlined  here  provide  a  reference  for  the  results  and  discussion 
of  the  measurements  taken.  The  specifics  of  each  case  are  included  and  generally  involve  the 
adjustment  of  the  intensity  level  of  the  beam  or  the  crystal  orientation  in  various  manners. 
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With  these  developments,  we  are  able  to  compare  the  experimental  and  theoretical  results 


under  different  conditions. 

i -Focussed  Second  Harmonic  Generation  with  A^GaSez 
ie  first  exposures  of  the  AgGaSe2  crystal  were  conducted  in  the  non-focussed  case,  15.2 
cm  from  the  output  port  of  the  laser,  with  an  average  input  power  of  approximately  8.2  Watts. 
Energy  measurements  were  limited  to  one-half  second  exposures,  and  in  this  case,  there  was 
not  a  second  energy  calorimeter  available  to  monitor  the  energy  at  the  fundamental.  The  effi¬ 
ciencies  extracted  were  referred  to  measurements  of  the  energy  at  the  fundamental  before  and 
after  the  experiment  without  the  crystal  in  place.  This  added  uncertainty  is  reflected  in  the 
calculated  error  bars  as  well  as  some  odd  behavior  in  the  data  points  shown  in  Figure  36. 

With  these  results,  we  can  compare  the  theoretical  estimates  of  the  conversion  efficiency 
developed  in  Chapter  IV.  We  note  that  the  simple  theory,  Eq.(4.1),  is  directly  applicable  while 
the  exact  solution,  Eq.(4.38),  must  be  modified  slightly  to  incorporate  the  angular  phase  mis¬ 
match.  This  can  be  accomplished  by  rewriting  the  dephasing  term,  Eq.(4.33)  as 


6  =>  \AkL  =  ±(fcge  +  2pe(0-6J) 


(6.1) 


Using  this  form,  Eq.(4.1)  and  Eq.(4.38)  can  be  directly  compared.  As  we  note  in  Figures  36 
and  37,  a  constant  multiplier  is  used  with  each  to  account  for  deviations  between  the  measured 
and  calculated  values.  Although  the  direct  results  of  these  calculations  are  on  the  same  order 
of  magnitude  as  the  data,  they  are  difficult  to  compare  to  the  data  without  these  correction 
factors. 

The  drop  in  the  efficiency  at  the  angle  0.01°  can  most  likely  be  attributed  to  a  variation  in 
the  laser  power.  Note  that  the  data  points  in  Figure  36  are  shown  a  function  of  the  deviation 
from  the  phase-matched  angle  wKich  was  assumed  to  lie  at  the  peak  of  this  pattern.  This  angle 
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Figure  36.  Efficiency  as  a  function  of  the  internal  angle  of  propag  .<  0,  using  ..ugy  mea¬ 

surements  of  the  non-focussed  second  harmonic  generation  in  AgG  SV:.  Distance  f  rom  tfie 
laser  was  15.2  cm,  with  an  average  input  power  of  9.53  W.  Tne  dac  points  are  shown  as  the 

filled  triangles  along  with  the  calculated  error  bars  (2:59).  ( - )  Fit  to  the  data  using 

Eq.(4.1)  modified  by  a  factor  of  0.27,  ( . )  Fit  to  data  using  Eq.(4.38)  modified  by  a 

factor  of  0.17. 


was  5.25°  ±0.05"  away  from  perpendicular  to  the  face  in  the  direction  of  the  optic  axis.  The 
manufacturer  cut  the  crystal  for  a  phase-matched  angle  of  57.17°,  which  would  lead  to  a  phase- 
matched  input  at  10.805  ^tm  (Note  calculations  in  Chapter  IV).  At  10.6  4m,  if  we  divide  5.25° 
by  the  index  of  refraction,  2.5915,  the  peak  of  the  efficiency  pattern  (55.14°)  agrees  well  with 
die  phase  matched  angle  of  55.02°  calculated  from  the  Sellmeier  equations.  The  oblique  angle 
of  incidence  pays  off  in  the  sense  that  there  is  no  feedback  to  the  laser  by  the  energy  reflected 
from  the  crystal  surface. 

The  asymmetric  appearance  of  the  measured  pattern  about  the  phase  matched  angle  can  be 
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accounted  for  in  a  number  of  ways.  The  fact  that  the  phase-matched  angle  does  not  lie  perpen¬ 
dicular  to  the  face  may  account  for  some  this,  but  the  majority  is  due  the  variation  in  the  aver¬ 
age  output  power  of  the  laser  source.  This  was  not  directly  measured  in  this  set-up  due  to  the 
lack"  of  a  second  energy  meter. 

Noting  that  the  short  exposures  had  no  damaging  effects  on  the  crystal,  the  same  experiment 
was  repeated  taking  power  measurements  instead.  In  this  case  we  were  able  to  track  the  input 
power  at  the  fundamental  using  a  second  power  meter.  This  removed  the  uncertainty  in  the 
input  power  levels  and  led  to  a  more  symmetric  set  of  points  about  the  phase-matched  angle 
(Figure  37). 

The  average  power  of  7.0  W  is  lower  than  the  previous  set  of  measurements,  but  the  effi¬ 
ciencies  are  comparable.  We  note  the  width  of  this  pattern  is  a  bit  larger  than  the  calculated 
value.  This  may  be  due  to  crystal  inhomogeneity  or  heating  affects  during  these  long 
exposures. 

The  width  of  the  pattern  calculated  using  the  exact  solution  drops  off  with  a  steeper  slope 
about  the  center  peak  than  Eq.(4. 1),  but  both  display  the  sine 2  dependence  found  in  the  original 
developments. 

All  of  the  calculations  fall  in  the  same  order  of  magnitude  of  as  the  measured  data,  but  each 
suffers  different  combinations  of  assumptions  and  approximations.  These  may  attributed  to 
several  sources.  The  published  values  for  the  nonlinear  optical  coefficients  vary  by  nearly  a 
factor  of  two  (1 1 :788).  At  oblique  angles  of  incidence,  less  energy  than  calculated  may  be 
passing  into  the  crystal,  although  the  anti-reflection  coatings  should  significantly  reduce  this 
problem.  Finally,  the  beam  may  be  chpped  by  the  physical  extent  of  the  crystal  face  resulting 
in  a  lower  overall  conversion  efficiency. 

If  we  include  the  radially  dependent  intensity  o.  .ic  in  the  exact  solution  (Eq.(4.56)),  the 
efficiency  calculated  for  an  average  power  of  7.0  W,  using  the  far-field  half  angle  divergence 
at  the  phase-matched  angle  is  0.29%.  At  these  extremely  low  values  of  efficiency,  this  esti- 
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Figure  37.  Efficiency  as  a  function  of  the  internal  angle  of  propagation,  0,  using  power  mea¬ 
surements  of  the  non-focussed  second  harmonic  generation  in  AgGaSe2 ■  The  crystal  was 
placed  15.2  cm  from  the  laser,  with  an  average  input  power  of  7.0  W.  The  data  points  are 

shown  as  the  filled  triangles  with  the  calculated  error  bars.  ( - )  Fit  to  the  data  using 

Eq.(4.1)  modified  by  a  factor  of  0.34.  ( . )  Fit  to  data  using  Eq.(4.38)  modified  by  a 

factor  of  0.22. 

mate  is  reasonably  dose  to  the  peak  value  of  0.1 1%  for  the  same  conditions  considered  in  the 
power  measurements.  This  calculation  proves  to  be  very  tedious  and  is  summarized  in  Figures 
38  and  39. 

In  Figure  38,  the  radially  dependent  conversion  efficiency  is  the  exact  solution  calculated 
given  the  intensity  level  at  the  specified  values  of  the  radius,  r,  in  Eq.(4.56).  The  weighting 
factor  is  the  integral  expression  of  this  equation.  The  integral  of  the  area  under  the  curve  in 
Figure  39  is  the  total  efficiency 

The  temporal  waveforms  of  the  fundamental  and  second  harmonic  pulses  at  the  angle  of 
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Figure  38.  Plot  of  the  components  of 

Eq.(4.56)  as  a  function  of  the  radius.  Figure  39.  Result  of  multiplying  the  compo¬ 

nents  point  by  point. 

peak  efficiency  at  the  same  distance  as  the  measurements  above  were  recorded  in  the  three 
pulse  repetition  formats  available  from  the  source.  Each  is  shown  in  Figures  40,  41  and  42  for 
the  30,  50  and  100  kHz  pulse  repetition  formats  respectively.  In  all  the  cases,  the  second  har¬ 
monic  pulse  increases  efficiency  with  the  instantaneous  power  of  the  pulse  to  peak  at  the  tem¬ 
poral  peak  of  the  fundamental. 
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Figure  40.  Normalized  amplitude  of  the  fundamental  and  second  harmonic  pulses  as  a  func¬ 
tion  of  time  for  the  30  kHz  pulse  repetition  frequency.  Average  input  power  =  7.4  W.  The 
Gaussian  fit  to  the  fundamental  pulse  has  a  characteristic  time  of  21.86  ns. 


Figure  41.  Normalized  amplitude  of  the  fundamental  and  second  harmonic  pulses  as  a  func¬ 
tion  of  time  for  the  50  kHz  pulse  repetition  frequency.  Average  input  power  =  1 1.6  W. 
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Figure  42.  Normalized  amplitude  of  the  fundamental  and  second  harmonic  pulses  as  a  func¬ 
tion  of  time  for  the  100  kHz  pulse  repetition  frequency.  Average  input  power  =  17.0  W. 

The  spatial  profiles  of  the  fundamental  and  second  harmonic  beams  recorded  using  the 
pyroelectric  infrared  vidicon  camera  are  shown  in  Figures  43  and  44.  Since  the  conversion 
efficiency  is  so  low,  there  is  little  unusual  structure  in  either  of  the  surfaces.  At  higher  effi¬ 
ciencies,  we  may  would  expect  to  see  asymmetric  depletion  of  the  fundamental  as  the  two 
beams  deviate  from  one  another  c  ver  the  length  of  the  crystal. 


Figure  43.  Surface  plot  of  the  intensity  pro¬ 
file  of  the  fundamental  after  passing  through 
the  AgGaSe2  crystal. 


Figure  44.  Surface  plot  of  the  intensity  profile 
of  the  second  harmonic  produced  by  AgGaSe2 
crystal. 
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Using  the  12.5  cm  focal  length  lens,  the  initial  focussed  efficiency  measurements  were 
taken  at  a  distance  which  would  avoid  any  chance  of  early  damage  to  the  crystal.  Second  har¬ 
monic  energy  produced  over  the  course  of  a  series  of  1  second  exposures  was  recorded  as  a 
function  of  the  angle  of  propagation  0  in  the  media  and  the  results  are  shown  in  Figure  45. 


Internal  Angle  [Degrees] 


Figure  45.  Efficiency  as  a  function  of  the  internal  angle  0  for  the  focussed  case,  f  =  12.5  cm, 
distance  from  the  laser  to  the  lens,  z:  =  50  cm,  distance  from  the  lens  to  the  crystal  face,  z}  = 
11.1  cm  and  a  nominal  average  power  of  6.0  W.  The  solid  line  is  an  eight  order  polynomial  fit. 

In  these  measurements,  versus  0,  for  the  non-focussed  and  focussed  cases,  the  angular 
acceptance  can  be  extracted  as  the  internal  angular  bandwidth  over  which  the  efficiency  drops 
to  a  value  of  one  half  the  peak.  For  the  three  methods  of  measurement,  the  results  are  shown  in 
Table  6. 
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Table  5. 


Phase  matched  angles  and  angular  acceptance 
results  for  the  three  measurement  techniques 
at  the  phasematched  angle,  0  =  55.02°. 


Divergence  of 
the  input  beam 

Non-focussed  energy  measurements 

0.292°  ±0.03 

1.85  mrad 

Non-focussed  power  measurements 

0.333°  ±0.04 

1.85  mrad 

Focussed  energy  measurements, 

0.656°  ±0.03 

10.39  mrad 

/=  12.5  cm 

With  the  sensitivity  of  this  axis  characterized,  we  were  able  to  begin  studies  of  the  depen¬ 
dence  of  the  conversion  efficiency  on  the  intensity  of  the  focussed  beam.  This  was  accom¬ 
plished  by  fixing  the  angular  position  in  the  phase-matched  direction  and  translating  the  crystal 
along  the  axis  z3  with  respect  to  the  focal  point.  This  was  accomphshed  for  several 
combinations  of  average  power  levels  and  the  two  lenses  used  in  the  focussed  experiments. 

For  the  12.5  cm  focal  length  lens,  the  efficiency  as  a  function  of  distance  behind  the  lens  is 
shown  in  Figures  46,  47  and  48  for  with  an  average  power  at  the  crystal  of  0.30  W,  2.22  W  and 
5.45  W  respectively.  The  peak  and  average  powers  were  attenuated  using  various  filters.  Note 
the  plot  of  the  peak  intensity  at  the  center  of  the  crystal  shown  against  the  secondary  y-axis  in 
each  of  the  figures.  It  is  clear  that  in  this  arrangement,  the  peak  efficiency  occurs  when  the 
focus  lies  at  the  center  of  the  crystal. 
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Figure  46.  Efficiency  as  a  function  of  the  distance  from  the  focussing  lens,  zh  f  =  12.5  cm,  at 
the  phase  matched  angle,  Bm.  Distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with  a  nominal 
average  power  of  0.30  W.  The  data  points  are  shown  as  the  filled  triangles  with  the  calculated 

error  bars.  ( - )  Fit  to  the  data  using  Eq.(4.38)  modified  by  a  factor  of  0.385.  ( . ) 

Calculated  peak  intensity  at  the  center  of  the  crystal  used  in  theoretical  calculation  of  the  effi¬ 
ciency. 
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Figure  47.  Efficiency  as  a  function  of  the  distance  from  the  focussing  lens,  z3,  f  =  12.5  cm,  at 
the  phase  matched  angle,  0„.  Distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with  a  nominal 
average  power  of  2.22  W.  The  data  points  are  shown  as  the  filled  triangles  with  the  calculated 

error  bars.  ( - )  Fit  to  the  data  using  Eq.(4.38)  modified  by  a  factor  of  0.5.  ( . ) 

Calculated  peak  intensity  at  the  center  of  the  crystal  used  in  theoretical  calculation  of  the  effi¬ 
ciency. 

Restricting  our  attention  to  the  case  of  highest  efficiencies,  the  calculated  values  can  be 
improved  by  incorporating  the  radially  dependent  intensity  profile,  Eqs.(4.55-56),  with 
Eq.(4.38)  as  shown  in  Figures  38  and  37.  Using  the  far  field  half  angle  divergence  of  the  beam 
inside  the  crystal  in  the  dephasing  term,  6  =  1.363rad2  (also  note  Eq.(4.48)),  we  can  calculate 
the  efficiency  as  a  function  of  the  distance  between  the  crystal  and  the  lens,  z},  to  yield  an 
excellent  approximation  to  this  efficiency  (see  Figure  49). 
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Figure  48.  efficiency  as  a  function  of  the  distance  from  the  focussing  lens,  z},  f  =  12.5  cm,  at 
the  phase  matched  angle,  9„.  Distance  from  the  laser  to  the  lens,  :2  ~  50  cm,  with  a  nominal 
average  power  of  5.45  W.  The  data  points  are  shown  as  the  filled  triangles  with  the  calculated 

error  bars.  ( - )  Fit  to  the  data  using  Eq.(4.38)  modified  by  a  factor  of  0.72.  ( . ) 

Calculated  peak  intensity  at  the  center  of  the  crystal  used  in  theoretical  calculation  of  the  effi¬ 
ciency. 

The  6.05  cm  focal  length  Ge  positive  meniscus  lens  had  the  tightest  focus  of  the  elements 
used  in  these  measurements.  In  its  use,  we  must  be  particularly  careful  about  the  position  of 
the  focus  and  the  corresponding  intensity  levels  to  avoid  damage  to  the  crystal.  The  minimum 
spot  size  of  55.61pm  is  much  smaller  than  optimum,  as  described  by  Eq.(4.63),  and  we  expect 
relatively  poor  efficiency.  The  results  for  an  average  input  power  of  0.19  W  and  0.40  W  are 
shown  in  Figures  50  and  51  respectively.  The  short  focal  length  leads  to  a  large  divergence 
and  the  beam  is  actually  clipped  by  the  10.6  pm  attenuator  prior  to  striking  the  detector  at  Posi¬ 
tion  I  in  Figure  34.  This  aperaturing  is  taken  into  account  when  calculating  the  power  in  the 
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Figure  49.  Measured  and  calculated  conversion  efficiency  (Eq.(4.56))  as  a  function  of  the 
position  of  the  focussed  beam  waist  with  respect  to  the  center  of  the  crystal,  f  =  12.5  cm,  at  the 
phase  matched  angle,  0m.  Distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with  a  nominal  aver¬ 
age  power  of  5.45  W.  ( . — )  Peak  intensity  at  the  center  of  the  crystal  (note  relation  to  z}  in 

Figure  48). 

second  harmonic  as  the  crystal  is  moved  in  relation  to  the  focus.  Although  the  peak  of  the 
efficiencies  still  corresp>ond  to  the  focus  at  the  center  of  the  crystal,  the  values  vary  a  great  deal. 
This  can  mostly  be  attributed  to  changes  in  the  output  power  of  the  laser  over  time. 
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Figure  50,  Measured  conversion  efficiency  us  a  function  of  ilic  distance  from  llic  focussing 
lens,  z,,  f  ■  6,05  cm,  at  die  phase  mulched  angle,  0m.  Distance  front  the  laser  to  the  lens,  i,  *» 

50  cm,  with  u  nominal  average  jmwer  of  0.10  W,  limn  bars  for  each  of  the  data  nolnts  are 
included  und  the  calculated  peuk  Intensity  at  the  center  of  the  crystal  nji|>cnr h  mi  lire  dnihed 
line. 

In  the  mathematical  estimates  of  the  efficiency,  all  ol  the  Input  |  >wn  at  the  fundamental 
frequency  is  assumed  to  be  parallel  to  polarization  required  lo  excite  me  dipole  oscillators  at 
the  frequency  of  the  second  harmonic.  Deviations  In  the  Input  polarization  result  In  less  power 
being  applied  to  this  solution,  The  power  ptoduced  at  the  second  harmonic  frequency  will  lull 
off  as  the  cos*  of  this  angular  deviation.  Wr  wetr  able  to  observe  this  effect  with  the  slioitci 
local  length  lens  by  totaling  the  crystul  about  Its  central  uxls  (axis  I.  In  f  igure  53)  'I  lie  results 
for  relatively  low  Input  power  levels  in  shown  In  Figure  52,  Note  the  |*uk  of  this  function 
docs  not  lie  parallel  the  vertical  |>oluii/.aiioii  axis  ((f)  and  can  l*e  nlUibulcd  lo  a  slight  deviation 
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f iguie  ,M.  Mcn«urc J  cum vci nUm  efficiency  an  a  function  of  the  dliluncc  from  (lie  foaming 
Jena, I  *  <»,()>  cm,  ai  the  phaar  matched  angle,  0„  DMtancc  fiom  the  Inacr  to  the  Iciik,  » 
50  un,  with  a  nominal  average  laiwcr  of  0,40  W,  liiror  ban*  for  each  of  (lie  duiu  ooinia  ure 
Included  and  die  calculated  peak  hiteiiilty  at  die  center  of  die  cryila)  appear*  ai  the  daubed 
line,  Note  die  average  |K»wcr  from  the  laicr  wav  probably  the  main  aourcc  of  variability  in 
dioe  efficiencies. 


of  die  output  polail/aiinn  of  U*e  iaici  from  veitlcal  at  thli  jm ilnt  (tee  figure  4),  'flic  drop  off  is 
quite  a  bit  iiecpct  than  unft  but  due  to  the  udd  entry  anglek,  accurate  0  alignment  wua  not 
always  pokklblc. 

In  several  nil  of  die  then  meakurementk,  panic ular  concern  wai  paid  to  l.ccplng  the  peuk 
and  average  (miwcu  well  below  known  damage  dueiliulds  foi  the  material.  My  Incorporating  u 
clioppei  Into  dir  ck|*crlmenial  vet  up,  a»  kbown  In  figure  34,  die  anticipated  high  average  inicn 
kliy  levcli  could  lie  reduced  to  kafe  IcveU  A  50'4  duty  cycle  cliop|*cf  blade  with  10  ujvcns  wav 
med  iu  chop  the  Ikcam  yi  frecpienclck  fmin  50  to  1000  II/.  <  ’bopping  die  Ivum  piovided  a  way 
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Figure  52.  Measured  and  calculated  conversion  efficiency  as  a  function  of  the  deviation  about 
axis  E,  f  «=  6.05  cm,  at  the  phase  matched  angle,  0m.  Distance  from  the  laser  to  the  lens,  z7  -  50 
cm,  distance  from  the  lens  to  the  crystal,  z3  =  6.3  cm,  with  a  nominal  average  power  of  0.19  W. 
Error  bars  the  measurements  are  shown  for  each  of  the  data  points. 

to  maintain  the  peak  intensity  values  of  the  pulses,  and  therefore  the  peak  conversion  efficien¬ 
cies,  while  reducing  the  thermal  load  (and  average  intensities). 

An  interesting  trend  was  observed  that  was  not  accounted  for  by  these  predictions.  There 
was  a  marked  decrease  in  the  observed  conversion  efficiency  when  the  beam  was  chopped  as 
opposed  to  allowed  to  run  unhindered.  The  first  set  of  measurements  were  taken  as  a  function 
of  Zj  for  several  chopping  frequencies  with  the  blade  described  above  and  is  shown  in  Figure 
53.  The  averuge  power  of  the  non-chopped  beam  was  3.95  W  during  this  set  of  measurements 
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and  we  note  the  measurements  at  the  distances  further  from  the  focal  length  become  indepen¬ 
dent  of  the  chopping  frequency.  The  readings  at  the  lower  chopping  frequencies  were  discon¬ 
tinued  due  to  a  lack  of  time. 


Focal  Position  wrt  Crystal  Center  [cm] 


Figure  53.  Measured  and  calculated  conversion  efficiency  as  a  function  of  the  distance  from 
the  focussing  clement,  z3,  for  several  chopping  frequencies,  f  =  12.5  cm,  at  the  phase  matched 
distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with  a  nominal  average  power  of  3.95 
W.  Error  bars  for  a  typical  measurement  are  shown  for  one  of  the  data  points. 

Because  of  the  slow  variations  in  the  output  power  of  the  laser  as  a  function  of  time  (due  to 
cooling  water  temperature,  power  surges,  etc.),  a  probe  HeNe  laser  and  photovoltaic  detector 
were  incorporated  to  measure  the  exposure  time  of  the  shutter  mechanism  directly.  This 
allowed  us  to  back  out  the  average  power  in  the  beam  when  it  struck  the  crystal  and  track  this 
variation  without  incorporating  errors  in  the  mechanical  shutter  exposure  time. 

This  also  allowed  us  to  investigate  the  effect  of  exposure  duration  on  conversion  efficiency 
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for  the  chopped  and  non-chopped  beam.  These  results,  for  one  chopping  frequency  are  shown 
in  Figure  54.  The  average  power  of  the  laser  is  shown  as  the  series  of  unfilled  symbols  corre¬ 
sponding  to  the  y-axis  on  the  right  side  of  Figure  54.  The  average  power  fluctuations  are  rela¬ 
tively  minor.  It  is  interesting  to  note  the  trend  in  the  efficiency  falls  off  with  longer  exposures. 
This  can  be  expected  since  the  induced  thermal  gradients  will  generally  lead  to  greater 
dephasing  of  the  fundamental  beam  over  the  length  of  the  crystal. 
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Figure  54.  Measured  and  calculated  conversion  efficiency  as  a  function  of  the  exposure  dura¬ 
tion  at  a  fixed  distance  from  the  focussing  element,  z3,  for  0  Hz  (no  chopping)  and  500  Hz 
chopping  frequencies,  f  =  12.5  cm,  at  the  phase  matched  angle,  distance  from  the  laser  to 
the  lens,  z2  =  50  cm,  with  a  nominal  average  power  of  3.0  W.  Error  bars  for  the  measurements 
are  shown  for  each  of  the  data  points.  The  hollow  symbols  of  the  same  shape  as  the  data  points 
correspond  to  the  average  input  power  of  the  laser  during  those  exposures. 
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T.he  drop  off  near  the  shorter  exposures  is  due  to  the  limitations  of  the  mechanical  shutter 
reaction  time,  however,  the  difference  between  the  chopped  and  non-chopped  measurements  is 
not  as  easily  explained.  The  corrected  average  powers  indicate  that  the  peak  intensities  for 
each  is  approximately  the  same,  but  when  chopped,  the  efficiency  drops  significantly.  Thermal 
lensing  may  help  to  reduce  the  divergence  of  the  beam  inside  the  crystal  to  a  limited  degree  in 
the  non-chopped  mode,  leading  to  better  overall  conversion. 

Focussed  Second  Harmonic  Generation  with  ZnGeP , 

As  an  estimate  of  the  second  harmonic  power  we  expected  to  observe  using  the  configura¬ 
tion  specified  in  Chapter  V,  we  can  use  Eq.(4.56)  where  the  effective  nonlinear  coefficient  will 
be 

=  lpm/V)sin(2  •  85.16°)cos(2  •  90°)  =  1 8.7pm/V  (6-2> 

Using  the  same  configuration  as  in  the  previous  experiments,  where  z2  =  65cm  and  z3  =  1 3.6cm, 

and  an  input  power  of  approximately  2  Watts,  the  calculated  efficiency  is  1 .4%.  We  noted  in 
the  comparison  of  the  calculated  values  at  the  lower  efficiencies  with  Eq.(4.56)  were  generally 
much  higher  than  actually  observed. 

No  second  harmonic  power  was  observed  using  ZnGeP 2  with  either  the  calorimeters  or  the 
very  sensitive  HgCdTe  detectors.  This  may  be  due  to  several  sources.  At  these  oblique  angles 
of  incidence,  15°,  the  amount  of  energy  coupled  into  the  crystal  through  the  uncoated  surface  is 
decreased  by  25%,  while  the  dephasing  introduced  at  this  angle  was  not  considered  in  the 
estimates  of  the  efficiency.  is  still  relatively  small  at  this  angle,  and  the  crystal  is  now 
critically  phase  matched  which  wouid  lead  us  to  reconsider  the  accuracy  of  the  temperature 
dependent  Sellmeier  coefficients  (Table  3). 

Continued  experimental  work  may  achieve  very  low  efficiencies  at  best  if  the  crystal 


temperature  is  increased  even  further.  The  oven  used  in  this  experiment  cannot  achieve 
temperatures  much  above  250°C. 


Planned  as  the  final  stage  of  the  series  of  experiments  to  characterize  the  AgGaSe2  crystal, 
damage  occurred  at  power  levels  below  those  anticipated  from  published  values.  The  damage 
to  the  main  crystal  is  summarized  first  followed  by  the  results  from  the  exposures  of  the  wit¬ 
ness  samples. 

Damage  to  the  crystal  occurred  at  three  separate  places  on  the  crystal  face  after  a  series  of 
exposures  prior  to  the  breakdown  of  the  surface.  Bulk  damage  was  not  observed  on  any  of  the 
damage  points,  and  in  fact,  several  data  points  in  the  previous  analysis  were  taken  in  the  unaf¬ 
fected  areas  after  the  surface  had  been  marred.  The  first  of  the  damage  points  was  at  the  center 
of  the  crystal  face  and  occurred  after  approximately  seven  exposures  on  the  same  order  of 
magnitude  as  the  published  damage  threshold  of  19.7  kW/cmJ  for  repetitive  >y  pulsed  sources 
(12:10). 

Figure  55  shows  the  entire  face  of  the  crystal  and  the  effect  of  the  anti  reflection  coating 
flaking  away  near  the  damage  point.  Figure  55  is  a  close-up  of  the  damage  point  itself. 


Figure  55.  Damage  to  the  face  of  the  AgGaSe2  crystal  occurred  at  distance  from  the  focussing 
element,  z}  =  14.05  cm,  f  =  12.5  cm,  at  9„,  distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with 
a  nominal  average  power  of  5.45  W  at  a  PRF  of  30  KHz.  The  tick  marks  across  the  top  of  the 
photograph  correspond  to  increments  of  0.5  mm. 
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Figure  56.  Damage  to  the  face  of  the  AgGaSe,  crystal  occurred  at  distance  from  the  focussing 
element,  z3  =  14.05  cm,  f  =  12.5  cm,  at  0m,  distance  from  the  laser  to  the  lens,  z2  =  50  cm,  with 
a  nominal  average  power  of  5.45  W  at  a  PRF  of  30  KHz.  The  tick  marks  across  the  top  of  the 
photograph  correspond  to  increments  of  0.5  mm. 

The  results  of  these  damage  points  are  shown  in  Table  7  and  include  the  peak  and  average 
powers  at  the  damage  point,  the  pulse  repetition  frequency,  the  spot  size  of  the  beam  and  an 
estimate  of  the  radius  of  the  bum  spot. 

Table  6. 


Results  of  the  damage  points  to  the  AgGaSe2 
crystal  (Damage  point  1  is  shown  in  the  Figures  55  and  56). 


Damage  Point 

Average 

Intensity 

(kW/cra2) 

Peak 

Intensity 

PRF 

(kHz) 

Spot  size 
lum} 

Bum  spot 
Radius 

(urn) 

1  (Figures  55-56) 

9.06 

7794 

30 

138 

379 

2 

4.59 

1458 

100 

221 

288 

3 

4.92 

4240 

30 

155 

252 

The  witness  samples  were  mounted  perpendicular  the  input  beam,  backed  through  the  focus 
of  the  lens  and  subjected  to  a  series  of  one  second  exposures.  At  the  point  of  damage,  the  spot 
size  and  intensity  levels  were  recorded  for  the  coated  and  uncoated  sides  of  the  samples.  The 
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results  are  summarized  in  Table  8.  The  damage  points  for  the  coated  surfaces  were  restricted 
to  the  first  surface  and  the  results  are  shown  in  the  upper  portion  of  the  table.  The  damage  to 
the  uncoated  surfaces,  however,  reached  through  the  depth  of  the  witness  samples  to  form  a 
bum  spot  on  the  opposite  side  as  well.  It  appears  that  the  bulk  damage  limits  were  rarely  met 
due  to  the  low  threshold  values  of  the  antireflection  coatings. 

Table  7. 

Results  of  the  damage  studies  with  the 
AgGaSe 2  witness  samples. 


PRF 

Average 

Intensity 

IkW/crA 

Peak 

Intensity 

Spot  size 

Bum  spot 
Radius 

Damage  Point 

(kHz) 

(MW/cm2) 

(pm) 

(pm) 

Antireflection 

30 

3.1  ±0.3 

2.7  ±0.3 

209  ±21 

304  ±  53 

coated  surface 

50 

4.6  ±  1.0 

1.6  ±0.3 

193  ±  17 

447  ±  1 8 

100 

5.0  ±  1.1 

1.6  ±0.4 

229  ±  33 

661 ±  125 

Uncoated  surface 

30 

2.4  ±0.7 

2.1  ±0.7 

349  ±  30 

607  ±  30 

50 

10.8  ±9.6 

3.7  ±  3.4 

228  ±  99 

571  ±67 

100 

10.8  ±  1.3 

3.4  ±0.4 

153  ±  10 

571  ±30 

Opposing  coated 

30 

2.4  ±0.7 

2.0  ±  0.5 

315  ±  30 

178  ±30 

surface 

50 

7.7  ±6.1 

2.7  ±2.1 

208  ±  83 

154  ±39 

100 

13.3 ±  i.S 

4.2  ±0.6 

124  ±9 

125  ±  18 

An  interesting  trend  in  the  size  of  the  damage  points  to  the  coated  surfaces  is  shown  in 
Figure  57.  The  average  damage  threshold  of  the  coated  surface  was  4.85  i  1.51  k\Y/cm2  (aver¬ 
age  intensity).  In  Figure  57,  the  ratio  of  the  measured  to  the  calculated  spot  size  gives  a  rough 
measure  of  the  level  above  the  damage  threshold  was  achieved  during  that  exposure.  The 
minimum  threshold  for  damage  is  approximately  2.5  kW/cm2.  The  damage  threshold  for  the 
bulk  can  be  inferred  from  the  three  highest  intensity  levels  achieved  during  these  exposures  to 
be  16  ±  8  kW/cm2  which  is  within  the  error  bounds  of  the  previously  published  levels. 

The  surface  of  the  crystal  has  a  threshold  for  damage  which  may  vary  a  functio:.  of  the 
fatigue  of  the  surface  induced  by  several  exposures  over  time.  This  was  observed  in  the  varia- 
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Figure  57.  Average  intensity  levels  for  damage  versus  the  ratio  of  the  measured  bum  radius  to 
the  calculated  spot  sizes  at  the  face  of  the  crystal.  This  gives  a  rough  measure  of  minimum 
levels  for  damage  to  the  antireflection  coated  surface. 


tion  of  the  damage  points  to  the  main  crystal.  The  peak  intensity  levels  during  the  exposures  of 
the  witness  samples  averaged  to  a  value  of  2.0  ±  0.5  MW/cm2  which  is  well  below  the 
published  limits  and  it  is  unlikely  that  the  peak  intensity  was  mechanism  for  damage. 
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VII.  Summary 


The  results  of  this  study  showed  that  high-average  power  second  harmonic  generation  using 
AgGaSe2  provides  an  efficient  source  of  mid-infrared  radiation  at  5.3  pm  The  work  with 
ZnGeP2  revealed  that  this  material  may  be  more  appropriate  for  other  applications  such  as  OPO 
or  SHG  at  shorter  fundamental  wavelengths.  The  TEMqq  output  mode  of  the  C02  laser  source 
at  a  pulse  repetition  frequency  of  30  kHz  (Gaussian  temporal  profile  with  FWHM  =  36.4  ns) 
was  used  in  all  of  the  measurements. 

The  absorptivity,  a  =  0.027cm'1,  of  the  AgGaSe2  crystal  was  very  low  compared  to  those 
used  in  previous  studies,  (1 1 :786)  and  (12:1).  Measurements  were  made  using  a  low  power 
CW  beam  10.6  (im.  The  lead  and  thorium  fluoride  based  X74  anti-reflection  coatings 
reduced  surface  losses  to  less  than  0.5%  at  both  the  fundamental  and  second  harmonic.  There¬ 
fore,  all  efficiency  measurements  are  given  as  external  values. 

The  non-focussed  efficiency  measurements  near  the  output  port  of  the  laser  were  conducted 
using  two  different  methods.  The  first  involved  measurement  of  the  energy  produced  at  the 
second  harmonic  during  0.5  second  exposures  of  the  crystal.  The  far-field  divergence  half¬ 
angle  of  the  beam  inside  the  crystal,  60  =  1 ,85mrad,  with  an  average  input  power  of  9.53  watts 
and  u>=  2.35mm  at  the  face  of  the  crystal,  yielded  a  peak  efficiency  of  0.12%.  The  phase 
matched  angle,  0m  =  55.14°,  and  angular  acceptance,  A9  =  0.292°,  compared  very  well  with  the 
values,  0„  =  55.02°  and  A0  =  0.300°.  calculated  using  the  Sellmeier  coefficients  (Table  3). 

The  crystal  was  provided  from  the  manufacturer  cut  such  that  0  =  57.17°.  Since  0„  =  55°, 
the  fundamental  was  incident  on  the  crystal  face  at  5.25°  with  respect  to  the  normal  direction. 
The  theoretical  estimates  of  the  conversion  efficiency  using  Eq.(4.1)  and  Eq.(4.38)  were  on  the 
same  order  of  magnitude  as  the  measured  values.  The  most  exhaustive  computation,  Eq.(4.56), 
incorporating  the  radiaily  dependent  intensity  profile  of  the  Gaussian  beam  yielded  a  peak  effi¬ 
ciency  estimate  of  0.29%  which  is  in  fair  agreement  at  these  very  low  values.  Some  of  the 
deviation  may  be  attributed  to  the  oblique  angle  of  incidence  as  noted  above. 
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The  second  method  of  measuring  the  non-focussed  efficiency  involved  recording  the  aver¬ 
age  power  of  the  second  harmonic  produced  during  30  second  exposures  of  the  crystal.  With 
an  average  input  power  of  7.0  watts,  all  other  factors  the  same  as  above,  the  efficiency  at  0m 
was  0. 1 1%.  Here  we  begin  to  notice  the  effects  of  heating  on  the  crystal  with  an  increase  in 
the  angular  acceptance,  A0  =  0.333°.  Theoretical  estimates  of  the  efficiency  were  comparable 
to  those  above. 

Comparison  of  the  temporal  profiles  of  the  fundamental  and  second  harmonic  revealed  the 
(P[a °)2  dependence  of  the  power  produced  at  the  second  harmonic.  This  lends  support  to 
approximation  incorporated  in  the  development  of  Eq.(4.57)  which  allows  the  CW  theory  for 
focussed  SHG  to  be  applied  to  the  pulsed  case  (given  a  Gaussian  shaped  temporal  profile). 

The  focussed  SHG  experiments  were  conducted  using  a  1 "  diameter  ZnSc  lens  placed  50  cm 
from  the  output  port  of  the  laser  such  that  50  =  10.39  mrad  inside  the  crystal.  Measurements 
of  the  angular  acceptance  with  the  focus  of  the  beam  far  from  the  center  of  the  crystal  showed 
an  increase  with  this  divergent  input  beam  to  A0  =  0.656°. 

Efficiency  was  recorded  using  energy  measurements  where  the  exposure  duration  was  gen¬ 
erally  limited  to  one  second  to  avoid  possible  damage  to  the  crystal  or  coatings.  At  0m,  the 
distance  between  the  focussing  element  and  the  crystal  face  was  varied  in  order  of  track  the 
dependence  of  the  power  produced  at  the  second  harmonic  as  a  function  of  peak  and  average 
intensity  incident.  Plots  of  the  results  versus  the  position  of  the  beam  focus  inside  the  crystal 
with  respect  to  the  center  of  the  crystal  showed  that  the  peak  occurred  when  the  two  coincided. 

The  highest  efficiency  recorded  was  20.7%  and  occurred  with  an  average  input  power  of 
5.45  watts  with  the  focus  at  the  center  of  the  crystal.  Note,  this  value  is  an  average  external 
SHG  conversion  efficiency.  At  this  position,  the  average  output  power  observed  at  5.3  pm 
was  1.13  watts.  External  efficiencies  as  high  as  18.7%  have  been  previously  reported  for  SHG 
of  C02  output  with  an  antireflection  coated  samples  of  AgGaSe2  (12:1).  This  work  represents  a 
2%  increase  in  the  absolute  efficiency  and  a  10%  increase  relative  to  th;  published  values.  A 
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longer  crystal  length,  lower  absorptivity  and  larger  focussed  spot  size,  =  1 25  pm,  were  all 
contributing  factors  to  this  improvement. 

The  theoretical  estimates  of  the  efficiency  using  Eq.(4.56)  yielded  excellent  agreement  with 
the  measured  values  as  well  as  the  shaped  of  function  produced  when  moving  the  position  of 
the  focus.  Estimates  using  the  numerically  derived  approximation  (Pascal  program  in  Appen¬ 
dix  B)  agreed  with  the  peak  of  the  tuning  curves  as  well. 

An  unusual  drop  on  in  energy  produced  at  the  second  harmonic  was  observed  when  a  50% 
duty  cycle  chopper  was  placed  in  path  of  the  fundamental.  The  purpose  of  the  chopper  was  to 
reduce  the  average  intensity  of  the  beam  on  the  crystal  while  maintaining  the  peak  intensity 
values  of  the  pulses.  Since  conversion  efficiency  is  solely  a  function  of  the  instantaneous 
power,  it  should  not  have  been  effected  by  the  chopper  (or  a  drop  in  average  power).  How¬ 
ever,  given  an  average  input  power  of  3.0  watts,  with  the  focus  at  the  center  of  the  crystal,  a 
15%  drop  in  the  efficiency  was  observed  when  the  chopper  was  in  place  (operating  at  500  Hz). 
Further  studies  are  proposed  to  determine  the  cause  of  this  decrease. 

Damage  studies  revealed  that  the  antireflection  coatings  placed  on  the  crystal  significantly 
lowered  the  thresholds  for  catastrophic  breakdown  of  the  surface.  Tne  damage  threshold  for 
average  intensity  was  approximately  2.5 kW/cm2  and  that  for  the  peak  intensity  was 
2.0  MW/cm2.  These  are  much  lower  than  previously  reported  values  for  the  uncoated  (1 1:788) 
and  coated  crystal  (12:5). 


In  the  continuation  of  this  work  with  AgGaSe2,  better  overall  efficiencies  may  be  achieved 
by  using  higher  peak  intensities  in  conjunction  with  a  larger  focussed  spot  sizes.  This  will 
serve  to  increase  the  overall  output  power  and  avoid  the  damage  thresholds  of  the  material. 

The  crystal  must  first  be  recoated  if  the  same  sample  is  to  be  used  in  these  studies.  Developing 
antireflection  coatings  which  do  not  lower  the  damage  thresholds  for  the  material  is  of  primary 
concern.  As  noted  in  the  short  series  of  studies  above,  the  bulk  threshold  for  damage  was 
rarely  encountered  due  to  the  weakness  of  the  coating. 

The  computational  effort  should  be  continued  to  complete  the  codes  required  numerically 
integrate  the  indefinite  integral  shown  in  Eq.(4.57).  It  should  be  robust  enough  to  work  over 
all  limits  of  the  integration  encountered  for  variouc  crystal  lengths  and  orientations.  This 
would  provide  a  very  useful  100I  in  determining  the  optimum  experimental  arrangements  for 
focussed  second  harmonic  generation  in  a  variety  of  crystals. 

Of  particular  interest  are  the  transient  heating  effects  on  the  conversion  efficiencies  as  noted 
with  the  introduction  of  the  chopper  in  the  Figures  53  and  54.  This  unexpected  drop  in  the 
second  harmonic  output  at  lower  average  powers  may  lend  clues  to  the  problems  of  scaling 
these  devices  to  higher  levels. 

The  work  with  ZnGeP2  for  doubling  of  10.6)im  radiation  would  appear  to  be  of  interest  to 
the  research  community  only.  The  attraction  of  a  non-critically  phase  matched  material  is 
countered  by  very  poor  performance  (low  d^j)  at  these  wavelengths.  For  practical  purposes  it 
is  much  better  suited  for  optical  parametric  amplification  at  shorter  wavelengths  in  the  mid  to 
far  infrared  spectrum. 

The  use  of  AgGaSe2  as  a  doubling  material  for  C02  lasers  is  a  promising  candidate  for  an 
efficient  laser  source  in  the  5  ^im  range.  With  optimum  conversion  efficiencies  in  excess  of 
20%  and  routine  efficiencies  above  10%,  the  material  possesses  the  qualities  of  hardness, 
resistance  to  humidity  and  damage  necessary  for  use  in  several  applications. 
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In  an  anisotropic  medium,  we  can  characterize  the  speed  of  propagation  of  different  fre 
quencies  and  polarizations  of  light  as  the  solutions  to  the  Fresnel  equations  In  the  medium.  In 


developing  these  solutions,  the  relevant  aspects  of  crystal  symmetry  and  the  theory  of 
anisotropic  crystals  arc  presented  below, 

Crystal  Symmetry 

The  notation  of  crystal  symmetry  is  often  referred  to  in  classifying  various  nonlinear  male 
rial  and  a  review  is  useful,  A  ciystallinc  solid  is  com|>oscd  of  an  orderly  array  of  Identical 
building  blocks.  The  requirement  (hat  the  structure  be  repetitive  restricts  the  shape  of  basic 
units  to  those  which  can  fill  a  three  dimensional  space  without  leaving  gaps.  There  are  seven 
basic  systems  which  meet  this  requirement  and  arc  dcscubcd  by  a  set  of  three  characteristic 
length*  und  respective  angles  (Note  Figure  58). 


F’igutu  .  3.  Notation  for  angles  and  dimensions  within  a  unit  cell  (IH  *>), 

In  the  system  of  the  crystallographic  ascs  (a,  b  and  c),  directions  and  planes  may  Ire 
described  in  term*  of  the  number  basic  unit  lengths  to  travel  In  cadi  of  (lie  three  (not  ncctssar 


Table  8.  The  Seven  Crystal  Syiicmi 


Syitcin 

Charade  riiilc 

Unit  Cell 

Cry*t«l  Point 

Charge  ter  iitic 

Group 

Trie  link 

None 

a  *»  b  *  c 

l.T 

Monocllnic 

One  2  fold  rotation 
aaii 

a  *  b 

a  -  (J«  VO*  n  y 

2,  m,  2/m 

Orthorhombic 

Three  mutually 

perpendicular  2  fold 
roiailun  aici 

u  =  (i  «■  y  *»  v(r 

222,  mm2,  mmm 

'Iriragonal 

One  4  fold  rotation 
ail*  or  m  4  fold 
roiaiion  mvcnlon 

Bill 

a  *  b  1 1 
<( «  (1  -  y  *  VO" 

4, 4,  4/m, 

422,  4mm, 

42m,  4/mmm 

Cubic 

I’otii  3  fold  rotation 
a*i»  (diagonal*; 

U  m  fr*  < 
u  -  (1  -*  y  -  VO" 

23.  mV  4 >2, 

3.3m,  m3m 

lleaugoiial 

One  f>  fold  rotation 

mil 

a  *»  h  ?  c 
a  *.  jU  VO" 
y-  120* 

6,  6,  6/m, 

622,  6mm,  6m2,  6/nimm 

'Irlpona/ 

One  3  fold  rotation 

a  *  b  •*  r 

X\  32, 

(Kknnbohedrali 

ml* 

a  »  ft »  y  ^  VO* 
u  -  ft  <•  y  •*  120* 

3m,  3m 

Adapted  (rum  Mi  Kclvcy  ( IK  5)  and  Wood  (2 1  ;S5) 


lly  orthogonal)  dlrcc  tlmn  and  the  reciprocal  of  (hcnc  value*  are  referred  to  a*  the  Miller 
imiiici.  An  example  of  »  direction  in  Milln  imlicn  would  be  1 1 1 2 j.  Wc  take  the  rcciprocul  of 
tin*  ivt,  w*  2,2, 1  ,  and  llri*  would  Ueimbe  a  vccloi  mapjvd  from  the  origin  of  the 
lyMcm  to  a  point  two  baik  unlu  along  Nil*  a,  two  along  Bill  b  and  one  along  ml*  c.  Simi 
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larly,  the  description  of  the  plane  (112)  would  intersect  the  respective  axes  at  each  of  these  unit 
distances  and  if  the  system  were  cubic,  would  lie  perpendicular  to  the  direction  [112).  Note 
Figure  59  where  this  direction  and  plane  are  described  for  a  system  with  a  cubic  structure. 


C 


Figure  59.  Miller  indices  and  planes. 

Each  of  these  basic  units  possesses  certain  elements  of  symmetry  which  include  mirror 
planes,  rotational  axes  and  inversion  axes  (21 : 10).  A  mirror  plane  exists  if  the  system  can  be 
divided  along  a  direction  that  would  allow  one  half  to  be  reflected  through  a  mirror  to  repro¬ 
duce  the  other.  A  rotational  axis  is  an  axis  around  which  the  system  may  be  rotated  to  repro¬ 
duce  a  system  identical  to  that  with  which  one  started.  These  are  limited  to  1,  2,  3,  4  and 
6-fold  axes  (generally,  icccnt  evidence  of  5-fould  symmetry  is  still  under  investigation).  An 
inversion  axis  (often  referred  to  as  an  improper  rotation  axis)  is  present  if  the  system  may  be 
rotated  by  2jt/n  and  then  inverted  through  the  center  of  symmetry  (center  of  the  basic  system 
or  reflected  through  the  horizontal  plane  of  the  system)  to  reproduce  the  original  system.  The 
combination  of  these  elements  of  symmetry  and  the  seven  basic  units  leads  to  the  32  crystallo¬ 
graphic  point  groups. 

In  this  study,  both  crystals  have  a  tetragonal  basis  and  fall  into  the  42m  point  group.  These 
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possess  three  two-fold  rotation  axes,  two  diagonal  reflection  planes  and  two  inversion  axes. 

We  note  that  these  axes  are  useful  for  describing  physical  distances  within  the  crystal  but 
will  often  lack  general  utility  because  they  are  not  always  a  mutually  orthogonal  set  which  can 
cover  all  of  three  dimensional  space  within  the  crystal.  For  this  reason,  we  refer  to  the  polariz¬ 
ability  of  the  material  since  it  is  a  physical  property  which  will  lend  itself  to  designating  an 
appropriate  coordinate  system. 


The  induced  polarization  of  an  isotropic  medium  is  always  parallel  to  the  direction  of  an 
applied  electric  Field  and  is  related  to  it  by  a  constant  factor,  the  susceptibility,  x,  (23:69)  and  is 
a  measure  of  the  net  displacement  of  the  electrons  in  a  material  under  the  influence  of  an  opti¬ 
cal  E  field. 

P  =  £oX£  (A1) 

In  isotropic  materials,  the  crystallographic  axes  are  a  mutually  exclusive  set  of  orthogonal 

directions.  In  an  anisotropic  dielectric  material  (where  we  assume  the  conductivity  and  free 

charge  to  be  zero),  however,  this  is  not  the  case.  The  displacement  current,  D,  may  be 

expressed  as 

~D  -  e^E  4-  ~P  (A.2) 

where  the  polarization  will  not  be  parallel  to  the  electric  Field  expect  for  certain  cases. 

The  linear  polarization  will  depend  upon  the  direction  arid  magnitude  of  the  applied  Field 
and  will  be  related  by 

=  EoOCnE,  +  Xi^>  +  Xi3^»)  (A. 3) 

Py  ~  Eo(X:.£,  +  Xu^y  XiyE,) 

-  Eo(X3.£,  +  X-yfiy  +  X33^|) 

where  P  and  E  are  complex  Field  quantities  and  the  electric  susceptibility  becomes  a  second 
order  tensor,  x  (in  tensor  notation,  x.,  ).  From  Eq.(A.2)  and  Eq.(A.3)  above,  we  may  define 
the  relative  dielectric  tensor,  e* 
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D  =  eJ  +  e^E  =  eo(l  +  x)£  =  (A4) 

The  nine  components  make  up  efi  arc  assumed  to  be  real  and  the  tensor  is  required  by  the 

conservation  of  energy  to  be  Hermitian  and  therefore  symmetric.  This  means  that  there  will  at 
most  be  six  independent  terms  that  make  up  the  second  order  tensor. 

i . •) 

e*  = 

Recall  from  matrix  theory  that  if  a  given  matrix,  B,  is  real  and  symmetric,  there  exists  an 
orthogonal  matrix,  A,  such  that 


£ll 

E12 

^13 

eI2 

&22 

^23 

^3 

£j3 

ABAt  = 


|%  0  ox 
0  V  0 

0  0  X, 


(A.6) 


which  is  a  diagonal  matrix  and  the  Vs  are  the  eigenvalues  of  matrix  B  and  the  row  vectors  of 
A  are  the  orthogonal  eigendirections  of  B. 

Since  the  relative  dielectric  tensor  is  real  and  symmetric,  a  coordinate  transformation  may 
be  found  which  allows  the  system  to  be  described  by  three  independent  components.  The 
reference  coordinate  system  is  established  by  chosing  the  axes  parallel  to  these  eigendirections, 
x,  y  and  z,  such  that  the  off-axis  components  vanish.  This  axis  system  is  called  the  principal 
coordinate  system  and  is  not  necessarily  parallel  to  the  crystallographic  axes.  In  the  principal 
coordinate  system  the  dielectric  tensor  appears  as 

(A.7) 


fe 

0 

0" 

V 

0 

o N 

II 

0 

0 

= 

0 

A 

0 

0 

c° 

0 

A, 

where  the  diagonal  components  are  related  to  the  indices  of  refraction  of  the  material  by 


e,  =  nl 
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With  this,  we  are  able  to  define  the  optic  axes,  although  their  significance  will  not  be 
established  until  the  analytic  solution  to  the  Fresnel  equations  is  shown  below.  When 

*  n2  *  the  two  optic  axes  lie  in  a  plane  containing  two  of  the  eigendirections  and 
the  crystal  is  biaxial.  Similarly,  when  nt  =  n 2  *  n3,  the  optic  axis  is  parallel  to  one  of  the 
eigendirections  and  the  crystal  is  uniaxial.  If  /»|  =  =  /i3,  the  crystal  is  isotropic. 

Waves  in  an  Anisotropic  Medium 

Now  that  a  description  of  the  medium  is  available,  it  is  possible  to  derive  its  effect  on  inci¬ 
dent  radiation  by  considering  the  solution  to  Maxwell’s  equations  given  here  in  their  general 
form 

V  •  D  =  p  (A. 8) 

-  (A.9) 

Vx£  =  — — 

dt 


V  •  B  =  0 


(A. 10) 


VxH 


(A.l  1) 


where  in  a  dielectric  medium,  p  =  0,  a  =  0  and  p  =  Mo.  We  assume  a  monochromatic  plane 
wave  of  frequency  w  propagating  in  the  anisotropic  medium  with  an  electric  field 

£g!'  <«-*•?))  (A.  12) 


and  a  magnetic  field 


(A.  13) 

Here  k  is  the  propagation  or  wave  vector,  k  =  (o )/c)n£,  where  n  is  the  index  of  refraction  to  be 


solved  for  and  /  is  the  unit  vector  in  the  direction  of  propagation.  Substituting  these  forms  into 
Maxwell’s  equations  for  an  anisotropic  dielectric  medium 

k  •  D  =0  (A. 14) 

kxE  =  (opj/  (A.  ID) 
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k'B  =0  (A.16) 

k  xH  =  -<S  (A. 17) 

where  the  geometric  relations  between  the  various  Held  quantities  is  shown  in  Figure  60.  Note 
5,  the  Poynting  vector,  is  defined  by  S  =E  xH,  and  is  the  direction  of  energy  propagation,  not 
the  unit  propagation  vector  f . 


Figure  60.  Representation  of  relations  between  the  field  vectors. 

To  obtain  a  relationship  between  the  E  and  D,  we  can  eliminate  H  from  Eq.(A.17)  by  sub¬ 
stitution  from  Eq.(A.15)  to  obtain 

k*k*E  =  -o )2p0D  (A.  18) 

By  the  use  of  the  vector  identity  A  x~B  *C  -B{A  •  C)-C  (A  •  B),  and  solving  forD 

—  1  --  -  (A. 19) 

D  =  —  [|ifc|  £  -  *(*•£)] 
to 

which  is  expressed  in  terms  of  the  propagation  vector  explicitly.  We  can  eliminate  the 
magnitude  of  this  component  and  express  D  in  terms  of  the  direction  of  propagation  alone, 
where  k  =  [Jt„  k Jk3],  by  the  relations 

->  to  k  (A. 20) 

\k\  =  -n  and  —  =  £  =  [s,,  s7,  j,1 
c  1*1  123 
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By  substitution  into  Eq.(A.19),  the  displacement  vector  can  be  expressed  as 

D  =  f^n\E  -  f  (f  •  £)]  (A  21) 

where  s  is  the  unit  vector  in  the  direction  of  propagation  of  the  plane  wave  and  n  is  the  index 

of  refraction  or  the  wave  property  in  the  f  direction. 

At  this  point  there  are  several  approaches  available  for  finding  the  solutions  to  this  equation 
(23:77),  the  n's  associated  with  the  two  possible  orthogonal  polarisations  of  the  incident  plane 
represented  by  E.  An  analytic  approach  involving  the  simultaneous  solution  of  three  equations 
followed  by  the  geometric  representation  is  the  most  convenient  of  these  methods. 

By  the  linearity  assumption,  D  =  £<jE.RE,  Eq.(A.21)  may  be  readily  expressed  in  terms  of 
the  material  property  er. 


e£RE  =  £^n2[E  -  /(/•£))  (A-22) 

Assuming  that  we  are  working  in  the  principal  coordinate  system  where  the  dielectric  tensor  is 

a  diagonal  matrix,  we  are  able  to  sum  like  vectorial  components  on  both  sides  of  this  equation 
and  represent  it  in  the  general  form  below  where 
5V  =  0  when  i  *  j  and  8 tJ  =  1  when  i  =  j, 


I  n2(8„  -  V,)]  E}  =  0 

j 

This  represents  a  set  of  three  homogeneous  equations  for  the  values  in  E,,  E2, 


(A. 23) 

and  £3,  one 


for  each  value  of  i  =  1,  2,  and  3.  These  may  be  written  in  full  form,  where  the  indices 

of  refraction  have  been  substituted  for  the  eigenvalues  £„  £5,  and  £3, 

[n,2-/i2(l -5,2)]£,  +  [rt2J,jrJ£2  +  [n2$,$3]£3  =  0 

(nVjE,  +  K2-V(l-s2)]£2  +  I/iV3]E3  =  0  (A,24) 

[n25,s3]£,  +  [n^jJEj  +  [n2-rt2(l  -532)]£3  =  0 


For  the  nontrivial  solution  to  this  set  of  homogeneous  equations  the  determinant  must  equal 
zero.  This  is  the  Fresnel  equation  (note  this  must  be  in  the  principal  coordinate  system) 
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0 


(A. 25) 


n^-n2{\  -s2) 
n2sts2 
n2sxs3 


n2sxs2  n2sxs  3 

nf  -  n2(\  ~ s2)  n 2s2si 

«V3  rt32-«2(l-5j) 


Given  a  particular  direction  of  propagation  with  respect  to  the  principal  axes,  the  determinant 
will  reduce  to  the  general  form  An6  +  Rn*  +  Cn2  +  D  =0  where  A  =0.  Thus  the  solution  is 
quadratic  in  n2  and  will  yield  at  most  two  independent  solutions,  the  eigenindices,  n1  and  n”. 
Substitution  of  each  of  these  solutions  into  Eqs.(2.24)  will  yield  the  eigenpolarization  direc¬ 
tions,  E?  and  Eu,  which  are  the  orthogonal  linearly  polarized  eigenmodes  for  this  direction  of 
propagation.  The  crystal  will  break  up  the  incident  radiation  into  these  two  orthogonal 
polarizations  with  respect  to  the  principal  coordinate  system.  Each  polarization  will  travel  at  a 
speed  dictated  by  the  respective  index.  The  direction  of  propagation  for  both  is  considered  to 
be  the  same,  the  direction  of  energy  flow  may  be  different. 

Recall  from  Figure  60  and  that  Maxwell’s  equations  imply  that  D ,  E  and  H  are  related 
by 


D 


x  H  and  H  =  —  /  x  E 
c  \ic 


(A. 26) 


According  to  these  equations,  D  and  E  are  both  perpendicular  to  the  direction  of  propagation  f 

and  the  resulting  direction  of  the  energy  flow,  the  Poynting  vector,  E  x  H  is  not,  in 
general,  parallel  to  f  (23:75). 

The  Index  Surface  Representation 

This  last  statement  is  easier  to  visualize  once  an  understanding  of  the  fact  that  the  different 
polarizations  of  light  within  the  media  travel  at  different  speeds  has  been  developed.  Consider 
that  the  index  of  refraction  for  each  of  the  two  eigenpolarizations  can  be  represented  as  a  three 
dimensional  surface  as  depicted  in  Figure  61.  Given  the  general  case  that  nx<n2<  n3,  these 
will  be  two  separate  index  surfaces  or  ellipsoids  which  are  the  normalized  k  surfaces  (recall 
Eq.(A.20)). 


108 


*2 


Figure  61.  First  quadrant  of  a  biaxial  index  surfaces. 

This  may  be  interpreted  best  by  considering  the  example  polarizations  shown  as  the  dashed 
lines  with  the  arrow  heads.  For  both  of  these  polarizations,  the  propagation  is  in  the  x1  direc¬ 
tion.  The  polarization  parallel  to  the  x2  direction  is  subject  to  the  n2  index.  Similarly,  the  ort¬ 
hogonal  polarization  (parallel  to  x3)  sees  the  n3  index.  This  is  the  most  general  case  of  a  biaxial 
crystal  which  has  two  optic  axes  as  defined  previously  and  shown  in  the  diagram  as  that  axis 
along  which  the  two  surfaces  meet.  For  a  biaxial  crystal,  there  are  only  four  such  points  which 
define  the  two  axes. 

In  the  course  of  this  study,  only  uniaxial  crystals  were  characterized,  which  is  the  case  when 
n\  =  rh*  ny  The  index  surfaces  then  share  a  common  boundary  in  the  x,-x2  plane.  In  addition, 
the  surfaces  will  be  ellipsoid*  of  revolution  about  the  x3  axis  (the  optic  axis)  and  the  surface 
defined  by  nt  and  n2  will  be  a  sphere.  and  n2  are  called  the  ordinary  indices  of  refraction,  n0, 
and  n3  is  called  the  extraordinary  index,  nt. 

It  can  be  shown  that  the  direction  of  energy  propagation  will  be  perpendicular  to  the  index 
surface  normal  defined  by  the  direction  of  the  propagation  through  the  media.  In  a  uniaxial 
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material,  Figure  62  depicts  the  two  index  surfaces  present  and  the  manner  in  which  the  ortho¬ 
gonally  polarized  components  of  the  incident  ray  are  broken  up  with  respect  to  optic  axis,  x3. 
Note  that  the  angle  of  incidence  with  the  crystal  interface  is  normal. 


WdkcK  1 
angle 


Figure  62.  Propagation  of  two  polarizations  of  an  incident  ray  in  a  uniaxial  material  where  the 
orientation  of  the  index  surfaces  produces  walkoff  of  the  polarization  in  the  plane  of  incidence 

The  deviation  of  the  energy  flow  of  one  of  the  polarizations  with  respect  to  the  other  is 
called  walkoff,  and  is  a  result  of  the  relations  between  the  Poynting  vector  and  propagation 
vector  in  an  anisotropic  media.  Keep  in  mind  that  the  normal  boundary  conditions,  Snell’s 
law,  is  still  applicable,  with  the  result  that  the  two  different  eigenpolarizations  will  be  deviated 
at  different  angles  due  to  the  two  eigenindices  as  shown  in  Figure  63.  This  is  an  example  of 
double  refraction. 

This  review  of  the  properties  of  anisotropic  material  properties  provides  the  basis  upon 
which  the  theory  of  nonlinear  optics  may  be  developed. 
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Figure  63.  Propagation  in  a  uniaxial  crystal  showing  the  effect  of  double-refraction  at  the 
interface  and  walkoff  of  the  polarization  in  the  plane  of  incidence  (in  the  plane  of  the  optic  axis 
and  the  propagation  vector). 
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Appendix  B;  Numerical  Computation  of  SHG  Efficiency 


This  program  follows  the  development  outlined  in  Chapter  IV  beginning  with  Eq.(4.57).  Tne 
indefinite  and  definite  integrals  are  evaluated  using  Romberg  integration  routines  given  the  user 
inputs.  The  values  calculated  agree  with  those  previously  published  (5:3600). 

The  program  was  written  in  Turbo  Pascal,  Version  5.5,  for  MS  DOS  based  machines.  The 
source  code  follows  a  brief  summary  of  the  required  inputs  and  operation  of  the  program  in  its 
present  form.  Reference  Chapter  FV  for  the  definitions  of  the  variables. 

Upon  compilation  and  execution  of  the  source  code,  you  have  the  option  to  specify  the  inputs 
for  o,  p,  K,  ^  and  p  in  a  number  of  ways.  When  the  program  is  initiated,  the  first  column 
of  characters  on  the  lower  left  side  of  the  screen,  as  shown  in  Figure  64,  will  accept  one  of  four 
characters  which  determine  how  the  values  are  calculated.  These  character  responses  are: 

ChauKULL  Description _ 

C  Calculate  the  value  directly  from  inputs  to  the  right  (note  definition  of  the  vari¬ 

ables). 

-  The  value  to  the  right  in  immediately  calculated  and  the  cursor  skips  down  to 
the  next  field. 

I  Input  the  value  directly  (be  wise  in  this  choice  due  to  the  interplay  between 

some  of  the  variables). 

-  The  cursor  moves  to  the  right  and  requests  the  desired  value. 

O  Optimum  value  for  the  parameter  (specified  in  Boyd  (5:3602)). 

-  Only  useful  for  o  at  this  point. 

R  Repeat  the  computation  over  a  range  of  values  for  the  specific  variable  (all  oth¬ 

ers  fixed). 

-  Only  use  this  option  for  one  of  the  five  variables  at  a  time! 

-  Must  be  aware  of  the  interplay  between  variables  in  this  case  as  well. 

-  When  you  have  completed  all  other  inputs,  the  program  will  ask  for  the  lower 
and  upper  limits  of  the  variable  inge  and  the  number  of  points  to  take  between 
*hose  limits. 

-  After  the  series  has  run,  the  program  asks  for  a  filename  under  which  to  save 
the  data  and  allows  you  to  continue. 

The  input  screen  appears  as  shown  in  Figure  64.  Upon  completion  of  the  computation,  the 
value  of  a,  P,  tc,  £  and  p.  and  the  conversion  efficiency  are  displayed  with  the  options  to 
continue. 
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Efficiency,  n,  as 

a  =  W  b  5k 

P  =  B /Vs 

§  =  L/b 
H  -  (L-2f)/L 
k  *  lfl(«l  -  l in  «2)b 
K  =  (2^2  w‘d‘/cA3 


function  of  the  parameters  defined  below  : 

b  *  confocal  parameter  ■  w*kl  dephseing  »  5k 

kl  “  prop,  oonet.  inside  w  «  minimum  spot  size 
B  ■  1/3  p  V  (L  kl)  p  ■  walkoff  angle  in  radians 

focussing  parameter 
f  ■  focal  position  (Mote  Fig.  1) 

«<1  e  abs.coef.  9  vl  «2  ■  abs.coef.  0  w2 

n*€l)R  P  m  Peak  power  ■  Bavg  /  (Vk  w  X) 

▼  >  pulse  rep.  freq.  t  ■  char,  time  of  pulse 


may  be  written  as  : 

n(0,P,k,§,n)  -  K  P  L  exp(-«'L)  h«J,P(k,S,H> 
where  the  Boyd  &  Kleinman  focussing  factor  is  defined  by  : 

h(C,P,k,§,)i)  -  (*,’/§)  exp  (2p.k$)  (2 /Vk)  f  «  |H(o  +  4ps,k,§,p)|*  exp(-4a>) 

Enter  all  parameters  in  8Z  units  (MKS)  :  Input  the  values  7  nr/Nj  Y 
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Lamb 

«  1 . IE-0005 

L  b  0.02900 

p  b 

0.01169 

1  I 

P  = 

0.00000 

nl 

■  2.59150 

5k  -  0.00000 

Favg  m 

5.000 

:  C 

s  « 

1.20246 

n2 

«  2.61600 

w  b  1.3E-0004 

V  ■ 

30000.0 

!  R 

p  = 

0.00000 

-1 

■  2.70000 

t  b  0.01450 

X  « 

2 . 2E-000  8 

;  c 

k  = 

0.01628 
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V _ 

Figure  64.  Appearance  of  the  input  screen  to  numerical  calculate  the  conversion  efficiency  fol¬ 
lowing  the  developments  in  Boyd  (5:3600)  and  Eckardt  (12:5). 
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PROGRAM  BKJ\.Factor;  (*  Borland  Turbo  Patcal  Vcriion  5.5  *) 

(*  Program  to  numerically  integrate  the  the  double  Integral  required  to 
estimate  the  SHG  power  a*  outlined  in  the  article: 

Boyd,  G.D.  and  D.A.  Kleinman,  "Parametric  Interaction  of  Focussed  Gauiilan 
Light  Beams, "  Journal  of  Applied  Physics,  Vol,  39.  No.  8;  3597-3639 
(July  1968). 


This  theory  takes  into  account  the  focussed  Gaussian  nature  of  (he  beam, 
wulkoff,  double  refraction,  position  of  the  focus  and  absorption  of  the 
fundamental  and  integrates  the  far  field  irradiance  pattern  to  evaluate 
the  total  second  harmonic  cower  generated.  The  original  program  was 
written  in  Pascal  for  the  Macintosh  by  Mr,  Robert  Eckardt  of  Standford 
University  using  a  combination  of  Romberg  and  Gaussian  quadrature  routines. 
It  converges  well  for  values  of  the  double  refraction  parameter  less  than 
two.  I  hope  to  improve  the  model  by  adding  other  types  of  integration  to 
extend  the  range  of  convergence. 

*) 


Uses  PR  INTER, CRY, 


CONST  PEEP  -  07; 

MAXNUMVAR  »  24; 

TYPE  STR 1  NO  1 5  -  STR I NO|  1 5 J; 

STR1NG60  »  $TRIN(}|60|; 

SLTOI  CHAR  -  SET  OF  CHAR; 

PARAMETERS  »  ARRAY  1 1., MAXNUMVAR)  OF  HEAL; 

VAR  Continue  :  Char; 

. . J 

PROCEDURE  RombH2(slgma, beta, *1, mu, kappa, s, lower, upper, tol'.real, 

q:intcgcr;  var  h;rcal);  FORWARD; 

(• . . - . *) 

JUNCTION  Or  tchar  (X,Y:INTF.aEK;QUICK,DFAl)LT;HOOLHAN; 

OKS  ET:SE'lt)l ’CHAR;  DIX’IIAR;  UI.INH:STKINOI5)  :CIIAR; 

(*  Rr’ads  the  response  at  the  screen  coordinates  X,Y,  QUICK  toggles  whether 
it  is  u  single  character  response  or  one  followed  by  a  CR,  DFAUi.T  toggles 
the  default  response  OF,  and  the  function  checks  to  see  if  the  input  Is 
in  the  acceptable,  set  of  characters,  SE'JX)FCIIAK.  *) 

VAR  CH  :  CHAR; 

SIIRT:. STRING  121; 

OOOD :  BOOLEAN; 

BEGIN  (•  (1ETCHAR  *) 

REPEAT 

GOTOXY  (X.Y);  WRITE  (F.UNE);  GOTOX Y  (X,Y), 

WRITE  (OF);  G07 OX Y  (X,Y); 

IE  QUICK  77 iEN  BEGIN 

CH  ReadKey;  GOTOXY  (X,Y);  WRITE  (CH)  END 
ELSE  BEGIN 
RE  ADEN  (SIIRT); 

IF  LENGTH  (SIIRT)  <  1  THEN  CH  >  '#’ 

ELSE  Cl  I  >  SIIR'l'l  !J; 
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END; 

CH  :*  UPCASE  (CH); 

IF  ((CH  IN  |CHR(  1 3),*#’J)  AND  (DFAULT  «  TRUE)}  THEN  BEGIN 
GOOD  :«  TRUE;  CH DP  END 
ELSE  GOOD CH  IN  OKSET; 

IP  NOT  GOOD  THEN  WRITE  (CHR(BEEP)); 

UNTIL  GOOD; 

GETCHAR  > CH; 

END,  (♦  GETCHAR  •) 

C . •) 

PROCEDURE  GetNumber  (X,Y:INTEGER;  VAR  DF:REAL); 

VAR  HOLD :  StrinfillO); 

IIOLDKcftl  :renl; 

GOOD :  BOOLEAN; 

IChir ;  CHAR; 

Numbers,  Arrow#,  Direct#*,  Option# ;  SETOPCHAR; 

Code  ;  Imcgri; 

BEGIN  (*  OETNumber  *) 

TcxtBackgroundtCyan); 

TcxtColor(Hhie); 

GOTOXY  (X.Y);  WRITE  C  *); 

GOTOXY  (X,Y);  WRITE  (DP:  10:5);  GOTOXY  (X,Y); 

RHADLN  (HOLD); 

IP  LENOni(HOLD)  <>  0  THEN  BEGIN 
V  AL(J  IOLD.I  lOLDRcal.Codc); 

IP  Code  -  0  THEN  DP  >  HOLD  Real; 

END: 

GOTOXY  (X.Y);  WRITE  (DP;10;5);  CKVOXY  (X.Y); 

Tcx.(  Huckgrounjf  Blnck ); 

TextColin  (While); 

END;  (*  GETN umber  •) 

(■'•••• . ♦) 

FUNCTION  fx2(x:rcnl;q:integeriigm*, beta, xl, my, kappa, c:n:al);rcal; 

(*  Evaluate#  the  rex)  and  imaginary  pam  ul  HA2  depending  upon 
the  toggle  q  *i 

VAR  Holdvttr,  JI2,KI»,Ih,  algr.taP.R.I :  real; 

BEGIN  (•  fx2  *) 

(IgniuP  :-  alama  <■  4  *  beta  *  s; 

IP  q  *  2  THEN  BEGIN 

R  ;-cxp(-kippa*xj/(l4-x*x)  *  (x*iln(iigmaP*x)'fco#(ilgmaP*x)); 
Moldvar :»  R  END 
ELSE  BEGIN 

I  ;•*  cxp(-kunpn*x)/(l+x*x)  *  (#in{iigmaP*x)-x*co»(#igfnaP*x)); 
HoldvHi  :» I  END; 
fx2  w  lluldvar; 

PND,  0  fx 2  *) 


FUNCTION  h.l(x:rcal;q:integrf;iiiigina,bctj,xi,mu(kappa:real):rcal; 
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(*  Function  that  retuins  the  value  of  the  function  F  and  makes 
two  calls  to  the  second  Romberg  integration  routine  which  separately 
evaluates  the  real  an  imaginary  parts  of  HA2  *) 

VAR  Holdvar,  H2,Rh,Ih,  sigmaP,R,I :  real; 

BEGIN  (*  fxl  *) 

(*  Calls  to  the  second  Romberg  integration  routine  RombH2  *) 
RombH2(sigma,beta,xi,mu,kappa,x,-xi*(l-mu),xi*(l+mu),lE-05,2,Rh); 
RombH2(sigma, beta, xi, mu, kappa, x,-xi*(l -mu), xi*(l+mu),lE-05,3,Ih); 
H2  :•  (Rh  *  Rh  +  Ih  *  Ih)/(4  *  pi  *pi); 

Holdvai  ;=  H2  *  exp(-4  *  x  *  x); 
fxl  :=  Holdvar, 


END;  (*  fxl  *) 

(• . J . *) 


PROCEDURE  ErrorMessage(q:integer); 

BEGIN  (*  ErrorMessage  *) 
TextColor(Yellow); 

GOTOXY(42, 16);  WritcC  ’); 

GOTOXY(42,16); 

CASE  q  OF 

1  :  write  (’Error:  Loop  for  exp(-4sA2)’); 

2  :  write  (’Error:  REAL  pan  of  HA2’); 

3  :  write  (’Error:  IMAGINARY  pan  of  HA2’); 
END; 

TextColor(White); 

Readln;  Exit 


END;  (*  En-orMessage  *) 

(* . *) 


PROCEDURE  RombH2;  (*  (sigma, beta, xi, mu, kappa, lower, upper, tokreal; 

q:integer;  varh.real);  *) 

(*  Romberg  integration  of  HA2  components,  declared  above  by  a 
FORWARD  declaration.  Both  of  the  Romberg  integration  routines, 
RombH2  and  Romb  are  identical  expect  for  the  calls  the  inner 
and  outer  functions.  They  were  taken  from: 

Alan  R.  Miller,  Pascal  Programs  for  Scientists  and  Engineers, 

Sybex  Inc.  (1981) 

The  only  change  is  in  the  convergence  criteria,  as  noted  below.  *) 

VAR  NX  :  anayll..  16]  of  integer; 

T  :  anay[1..136]  of  real; 
done,  error :  boolean; 

steps,  NT,  I,  II,  N,  NN,  L,  NTRA,  K,  M,  J  :  integer; 
delta_x,  C,  sum,  fotom,  X  :  real; 

BEGIN  (*  RombH2  *) 
done  :=  false;  error  :=  false; 
steps  :=1;  NX(1]  :=  1; 

delta_x  :=  (upper-lower)/steps; 

C:=  (fx2(lower,q,  sigma,  beta,  xi.mujcappa.s) 

+  fx2(upper,q,  sigma,  beta,  xi,  mu,  kappa,  s))*0.5; 

T[l]  :=  delu_x*C; 

N  :=  1;  NN  :=2;  sum  :=  C; 
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REPEAT 
N:=  N+l; 
fotom  :=  4.0; 

NX[N]:=  NN; 
steps  :=  steps*2; 

L:-  steps- 1; 

delta_x:=(upper-lower)/steps; 

FOR  D  :=  1  TO  (L+l)  div  2  EX)  BEGIN 
I  :=  II  *  2  - 1; 

X:=  lower  +  I*delta_x; 

sum  :=  sum+fx2(X,q,sigma,beta,xi,mu, kappa, s)  END; 

T[NN]  :=  delta  x*sum; 

NTRA  :=  NX[N-1]; 

K:=N-1; 

FOR  M:=  1  TO  K  DO  BEGIN 
J:=  NN  +  M; 

NT  :=  NX[N-1]+M-1; 

T[J]:=(fotom*T[J- l]-T[NT])/(fotom- 1 .0); 
fotom  :=  fotom*4.0  END; 

IF  N>4  THEN  BEGIN 
EFT[NN+1]  oO.OTHEN 

IF  (abs(T[NTRA+l]-T[NN+l])<=abs(T[NN+l]*tol)) 

OR  (abs(T[ N N- 1  ] -T[ J] )  <=  tol)  THEN  done  :=-•  true 

(*  In  the  source  for  this  code  segment,  the  tolerance  is  checked  according  to 
OR  (abs(T(NN-l]-T(J])  <=  abs(T[J]*tol))  THEN  done  :=  true 
instead  of  the  above.  I  found  that  the  routines  would  not  converge  for 
imaginary  part  of  the  HA2  because  it  is  often  very  close  to  zero  to 
begin  with.  Multiplying  T[ J]  by  tol  drives  the  value  so  close  to 
zero  that  we  start  seeing  limits  in  the  REAL  variable  type  (16  bits). 

The  results  with  this  converge  test  works  great.  *) 

ELSE  IF  N  >  15  THEN  BEGIN 
done:=true;  error  :=  true 
END; 

END;  (‘if*) 

NN  :=  J+l; 

UNTIL  done; 
h  :=T[J]; 

IF  error  THEN  ErrorMessage(q); 

END;  (*  RombH2  *) 

(* .  *) 

PROCEDURE  Romb  (sigma,  beta,  xi,mu,kappa,lower.upper,tol:real; 

qrinteger;  var  h:real); 

VAR  NX  :  array[1..16]  of  integer; 

T  :  array[1..136]  of  real; 
done,  error :  boolean; 

steps,  NT,  I,  II,  N,  NN,  L,  NFRA,  K,  M,  J  :  integer; 
delta_x,  C,  sum,  fotom,  X  *  real; 


117 


BEGIN  (*  Romb  *) 
done  :=  false;  error  :=  false; 
steps  :=  1;  NX[1]  :=  1; 

delta_x  :=  (upper-lower)/steps; 

C:=  (fxl(lower,q, sigma, beta, xi,mu,kappa) 

+  fxl(upper,q,sigma,beta,xi,mu,kappa))*0.5; 

T(l]  :=delta_x*C; 

N  :=  1;  NN:=2;  sum:=C; 

REPEAT 
N:=  N+l; 
fotom  :=  4.0; 

NX[N]:=  NN; 
steps  :=  steps*2; 

L:=  steps- 1; 

delta_x:=(upper-lower)/steps; 

FOR  D  :=  1  TO  (L+l)  div  2  DO  BEGIN 
I  :=  13  *  2  - 1; 

X:=  lower  +  I*delta_x; 

sum  :=  sum+fxl(X,q, sigma, beta, xi, mu, kappa)  END; 

T[NN]  :=  delta  x*sum; 

NTRA  :=  NX[N-1]; 

K:=N-1; 

FOR  M:=  1  TO  K  DO  BEGIN 
J:=  NN  +  M; 

NT  :=  NX[N-1]+M-1; 

T[J]:=(fotom*T[J- 1  ]-T[NTJ)/(fotom- 1 .0); 
fotom  :=  fotom*4.0  END; 

IF  N>4  THEN  BEGIN 
IF  T[NN+1]  oO.O  THEN 

IF  (abs(T[NTRA+ 1  ]-T[NN+l  ])<=abs(T[NN+ 1  ]*tol)) 

OR  (abs(T[NN-l]-T[J])  <=  abs(T[J]*tol))  THEN  done  :=  tme 
ELSE  IF  N  >  15  THEN  BEGIN 
done:=true;  error  :=  true 
END; 

END;  (*if*) 

NN  :=  J+l; 

UNTIL  done; 
h  :=T[J]; 

IF  error  THEN  ErrorMessage(q); 

END;  (*  Romb  *) 

(* . *) 

PROCEDURE  InputParms(VAR  sigma:real;  VAR  beta:real;  VAR  xi:real; 
VAR  mu:real;  VAR  kappa:real; 

VAR  FwdConst:real;  VAR  DFV PARAMETERS); 

(*  Allows  you  to  input  some  or  all  of  the  variables  used  to  calculate 
the  efficiency  *) 

VAR  Entervalues  :  CHAR; 

DONE :  BOOLEAN; 

P,  K,  Cst :  real; 

(*  . *) 
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PROCEDURE  DfaultVals; 

BEGIN  (*  DfaultVals  *) 

DFV[1]:=  1.0;  DFV[2]:=  1.0;  DFV[3]:=  1.0;  DFV[4]:=  1.0;  DFV[5]:=  1.0; 
DFV[6]  :=  0.57;  DFV[7]  :=  0.0;  DFV[8]  :=  2.84;  DFV[9]  :=  0.0; 

DFV[10]  :=0.0;  DFV[11J  :=  10.6E-6;  DFV[12]  ;=  2.5915;  DFV[  13]  :=  2.616; 
DFVJ14]  :=  2.7;  DFV[15j  :=  2.7;  DFV[16]  :=  0.029;  DFV[17]  :=  0.0; 

DFV[18]  :=  125.3E-6;  DFV[19]  :=  DFV[16]/2; 

DFV[20]  :=  (43c- 12)  *  Sin(55.02*Pi/180);  DFV[21]  :=  0.67*Pi/180; 

DFV[22]  :=  5.0;  DFV[23]  :=  30000.0;  DFV[24]:=21.86E-9; 

END;  (*  DfaultVals  *) 

(* . *) 

PROCEDURE  TextBox(UpleftX,  UpLeftY,  LoRightX,  LoRightY  :  integer); 
const 

UpLfCh  =#201;  LoLfCh  =#200; 

UpRtCh  =  #187  ;  LoRtCh  =#188  ; 

HorizCh  =  #205  ;  VertCh  =#186; 

var 

X,  Y,  NumToDoX,  CursorX,  CursorY;  integer; 
begin 

if  (UpLeftX  >=  LoRightX)  or  (UpLeftY  >=  LoRightY)  or 
(UpLeftX  <  1)  or  (LoRightX  >  mem[0:$44A])  or 
(UpLeftY  <  1)  or  (LoRightY  >  25)  then 

exit; 

CursorX  :=  wherex; 

CursorY  :=  wherey; 
gotoxy(UpLeftX,  UpLeftY); 
write(UpLfCh); 

NumToDoX  :=  LoRightX  -  UpLeftX  -  1; 
for  X  :=  1  to  NumToDoX  do 
write(HorizCh); 
write(UpRtCh); 

for  y:=  UpLeftY  +  1  to  LoRightY  -  1  do 
begin 

gotoxy(LoRightX,  Y); 
write(VertCh); 
gotoxy(UpLeftX,Y); 
write(VertCh); 
end; 

gotoxy(UpLcftX,  LoRightY); 
write(jjoLfCh); 
for  x  :=  1  to  NumToDoX  do 
write(HorizCh); 
write(LoRtCh); 
gotoxy(CursorX,  CursorY) 
end; 


{ 


} 


PROCEDURE  Write Equs; 

BEGIN  (*  WriteEqus  *) 

Clrscr, 

Writeln;  Textcolor(Green); 

Writeln(’EfTiciency,  \CHR(252),\  as  function  of  the  parameters  defined  below 
TextColor(White); 

Writeln(CHR(229), ’  =  ’,CHR(171),’  b  ’,CHR(235),’k’, 

’  b  =  confocal  parameter  =  w’,CHR(253), 

’kl  dephasing  =  ’,CHR(235),’k’); 

Writelnf  kl  =  prop,  const,  inside  w  =  minimum  spot  size’); 

Writeln(CHR(225),’  =  B/’,CHR(251),CHR(21),’  B  =  ’,CHR(171), 

’  p  \CHR(251),’(L  kl)  p  =  walkoff  angle  in  radians’); 

Writeln(CHR(21),’  =  L/b  focussing  parameter’); 

Writeln(CHR(230),’  =  (L-2f)/L  f  =  focal  position  (Note  Fig.  1)’); 

Writeln(’k  =  \CHR(171),’(’,CHR(224),’l  -  ’,CHR(171),’  ’,CHR(224),’2)b’, 

’  \CHR(224),’l  =  abs.coef.  @  wl  ’,CHR(224),’2  =  abs.coef.  @  w2’); 

WritelnfK  =  (2’,CHR(251),’2  w\CHR(253),’d’,CHR(253).7cA3  n’,CHR(253), 
CHR(238),’1)R  P  =  Peak  power  =  Pavg  /  (’,CHR(251),CHR(227),’  v  \CHR(231)/)’); 
Writeln(’  v  =  pulse  rep.  freq.  \ 

CHR(23l),’  =  char,  time  of  pulse’); 

TextColor(Cyan);  Writeln(’may  be  written  as :  ’);  TextColor(Green); 

Writeln(’  ’lCHR(252),’(’,CHR(229),’,’,CHR(225),’,k,’,CHR(21)t,,',CHR(230), 

’)  =  K  P  L  exp(-’,CHR(224),’”L)  h(’,CHR(229),’,’,CHR(225),’.k.’,CHR(21)> 
’,’,CHR(230),’)’); 

TextColor(Cyan); 

Writeln(’where  the  Boyd  &  Kleinman  focussing  factor  is  defined  by  :’); 

Textcolor(Green); 

WritclnC  h(’,CHR(229),’,’,CHR(225),’,k,’,CHR(21),’,’,CHR(230),’)  =  ’, 
’(’,CHR(227),CHR(253),7’,CHR(21),’)exp(2’,CHR(230),’k’,CHR(21),’)(2/’t 
CHR(25 1  ),CHR(227),’ )  ’,CHR(244),’  ’,CHR(236),’  ’.CHR(179),’H(’.CHR(229). 

’  +  4’,CHR(225),’s»k,’,CHR(21),’,’fCHR(230),’)’,CHR(179),CHR(253), 

’  exp(-4s’,CHR(253),’)  ds’); 

Writeln(’  ’,CHR(245).’-',CHR(236)); 

TextColor(White); 

Writeln(’Enter  all  parameters  in  SI  units  (MKS) ;’); 

END;  (*  WriteEqus  *) 

{  . ) 

PROCEDURE  Write  Values; 

VAR  I :  integer, 

(* . *) 

PROCEDURE  LineupChar; 

BEGIN  (*  LineupChar  *) 

GOTOXY  (7,18);  Write(CHR(229),’  =  ’);  GOTOXY(24,18);  Write(’Lamb  =  ’); 
GOTOXY  (43,18);  Write(’  L  =  ’);  GOTOXY(62,18);  Write(’  p  = ’); 

GOTOXY  (7,19);  Write(CHR(225),’  =  ’);  GOTOXY(24,19);  Write(’  nl  =  ’); 

GOTOXY  (43,19);  Write(CHR(235),’k  =  ’);  GOTOXY(62,19);  Write(’Pavg  =  ’); 
GOTOXY  (7,20);  Write(CHR(21),’  =  ’);  GOTOXY(24,20);  Writef  n2  =  ’); 

GOTOXY  (43,20);  Write(’  w  =  ’);  GOTOXY(62,20);  Write(’  v  =  ’); 

GOTOXY  (7,21);  Write(CHRv230),’  =  ’);  GOTOXY(24,21);  Write(’  \CHR(224),’l  =  ’); 
GOTOXY  (43,21);  Write(’  f  =  ’);  GOTOXY(62,21);  Write(’  ’,CHR(231),’  =  ’); 
GOTOXY  (7,22);  Write(’k  =  ’);  GOTOXY(24,22);  Write(’  ’,CHR(224),’2  =  ’); 
GOTOXY(43,22);  Write(’deff  = ’); 


120 


GOTOXY(3,24);  Writc(’K  =  ’);  GOTOXY(  19,24);  Write(’P  =  ’); 
GOTOXY(35,24);  Writc(’C  =  ’); 

END;  {*  LineupChar  *) 

(* . . *) 

PROCEDURE  Calc  Vais  (VAR  Peak:Real;  VAR  KConst:Real;  VAR  LConst:Real); 
VAR  enter :  integer, 

Db,  b,  alphal,  alpha2,  rho  :  real; 

BEGIN 

b  :=  Sqr(DFV[18])*2*Pi*DFV[12]/DFV[l  1]; 
alphal  :=DFV[14]-DFV[15]/2; 
alpha2  :»  DFV[14]  +  DFV[15]/2; 

Db  :=  DFV[21]*Sqrt(DFV[16)*2*Pi*DFV[12]/DFV[l  l])/2; 

KConst  :=  2*Sqrt(2)*Sqr((2*Pi*3E8/DFV[l  1])*DFV[20]/(DFV{  13]*3E8)); 
KConst  :=  KConst/(8.8542E- 1 2*3E8*DFV[  11]); 

L Const  :=  DFV[16]*exp(-alpha2*DFV[16]); 

Peak  :=  DFV(22]/(Sqrt(Pi)*DFV[23]*DFV[24]); 

FOR  enter  :=  1  TO  5  DO  BEGIN 
IF  DFV[cnter]  <>  1 .0  THEN  BEGIN 
CASE  enter  OF 

1  :  DFV[6]  :=  b*DFV[17]/2; 

2  :  DFV[7]  :=  Db/Sqrt(DFV[8]); 

3  ;  DFV[8]  :=  DFV[  16j/b; 

4  :  DFV[9]  :=  (DFV[16]  -  2*DFV(19])/DFV(16]; 

5  :  DFV[10]  :=  alpha l*b/2; 

END; 

END; 

END; 

IFDFV[1]  =  2.0  THEN 

DFV[6]  :=  1/DFV[8]  *  (cxp(-Sqrt(DFV[8])/2)  *  arctan(DFV[8])  + 

(1  •  exp(-Sqrt(DFV[8])/2))  *  2  *  arctan(DFV[8]/2)); 


END; 

(* . *) 

BEGIN  (*  WriteValues  *) 


Calcvals(P,K,Cst); 

FOR  I  :=  1  TO  7  DO  BEGIN 
GOTOXY(  1,16+1);  ClrEol  End; 

TextColor(Red);  TextBox(l,17,80,23);  TextColor(White); 
LineupChar; 

TextBackground(Black);  TextColor(White); 

FOR  I  :=1  TO  5  DO  BEGIN 
GOTOXY(3,17+I); 

CASE  TRUNC(DFV[I])  OF 
0 :  Write(’C'); 

1  :  BEGIN 

TextBackground(Cyan);  TextColor(Blue); 

Write(T);  GOTOXY  (13,17+1);  Write(DFV[5+I]:10:5); 
TcxtBackground(Black);  TextColorOV'hite); 

END; 

2  :  Write(’O’); 

3  :  Write(’R’); 

END; 

GOTOXY(  13,17+1); 
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IFDFV[I]  o  1.0  THEN  Write(DFV[5+I]:10:5); 

END; 

TcxtBackground(Cyan);  TcxtColor(Bluc); 

GOTOXY(30,18);  Writc(DFV[ll]:10);  GOTOXY  (49,18);  Write(DFV[16]:10:5); 
GOTOXY(68,18);  Write(DFV[21]:lO;5); 

GOTOXY(30,19);  Write(DFV[12]:lO:5);  GOTOXY  (49,19);  Write(DFV[17]:10:5); 
GOTOXY (68, 19);  Writc(DFV[22]:10:3); 

GOTOXY (30,20);  Write(DFV[13]:10:5);  GOTOXY  (49,20);  Write(DFV[18]:10); 
GOTOX Y(68 ,20);  Write(DFV[23] :  1 0: 1 ); 

GOTOXY(30,21);  Write(DFV(14]:l0:5);  GOTOXY  (49,21);  Write(DFV[19]:10:5); 
GOTOXY(68,2 1 );  Write(DFV[24] :  1 0); 

GOTOXY (30,22);  Write(DFV[15]:10:5);  GOTOXY  (49,22);  Write(DFV[20]:10); 
TextBackground(Black);  TextColor(White); 

GOTOXY(7,24);  Write(K:10); 

GOTOX Y(23, 24);  Write(P:10); 

GOTOXY  (39,24);  Wriie(Cst:10); 

END;  (*  WriteValues  *) 

{ . ) 


PROCEDURE  GetKeys(VAR  Complete:BOOLEAN); 

VAR  Arrows,  Numbers,  Direct,  Characters :  SETOFCHAR; 
TEMP :  CHAR; 

X,Y,I :  Integer; 

BEGIN  (*  GetKeys  *) 

DONE False; 

Characters  :=  (T.’i’.’C’/c’.’O’.V.’R’.V]; 

FOR  I  :=  1  TO  5  DO  BEGIN 
CASE  TRUNC(DFV[I])  OF 


0  :  TEMP 

1  :  TEMP 

2  :  TEMP 

3  :  TEMP 
END; 


=  ’C’; 
=  T; 
=  ’O’; 
=  ’R’; 


TEMP  :=  GETCHAR(3,17+I,Tnie,True, Characters, TEMP, TEMP); 
CASE  TEMP  OF 


’C’  :  DFV[I]  :=  0.0; 
’I’ :  BEGIN 


DFV[I]  :=  1.0; 

GetNumber(13,17+I,DFV{5+IJ); 

END; 

’O’  :  DFV[I]  :=  2.0; 

’R’  :  DFV[I]  :=  3.0; 

END; 

WriteValues; 

END; 

FOR  I  :=  1  TO  5  DO  BEGIN 


GetNumbcr(30, 1 7+I,DFV[  10+1]); 
WriteValues; 

END; 

FOR  I  :=  1  TO  5  DO  BEGIN 


GeiNumber(49,17+I,DFV(15+I]); 

WriteValues; 

END; 

FOR  I  ;=  1  TO  4  DO  BEGIN 
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GetN  umber(68 , 1 7 +1  ,DFV  [20+1]); 

WritcValucs; 

END; 

TextColor(Y  ellow); 

GOTOXY(51,24);  Write(’Done  ?  [Y/N]  ’); 

TEMP  ;=  GetChar(69,24,Truc,True,[,Y,,y ,’N\ Vl/YW); 
IFTEMP  =  ’Y’  THEN  Done:=True; 

TextColor(White); 

END;  (*  GetKeys  *) 

{ . ) 

BEGIN  (*  InputParms  *) 

TcxtBackground(Black);  TextColor(White); 

DfaultVals; 

WriteEqus; 

Write  Values; 

TextColor(Y  ellow); 

GOTOXY(42,16);  Writcflnput  the  values  ?  [Y/N]  *); 

EnterValues  :=  GETCHAR(70,16,True,True,rY7y7N\V],’Y7Y’) 
TextColor(White); 

IF  EnterValues  =  ’Y’  THEN  BEGIN 
REPEAT 
Getkeys(DONE) 

UNTIL  DONE; 

END; 

FwdConst  P*K*Cst; 

Sigma  ;=  DFV[6]; 

Beta  :=  DFV[7]; 
xi  DFV[8]; 
mu  ;=  DFV[9]; 

kappa  :=DFV[10]; 

END;  (*  InputParms  *) 

(* . *) 

PROCEDURE  Calllnteg; 

(*  Calls  the  numerical  integration  routines  given  the  user  inputs  *) 

VAR  sigma,  beta,  B,  xi,  mu,  kappa,  h,  FwdConst  :  real; 
ok  :  char, 

DFV  :  PARAMETERS; 

Cycle  :  BOOLEAN; 
ent ;  integer, 

(* . *) 

PROCEDURE  CycleThru; 

VAR  cyclevar,numcycles,cyclecount :  integer; 
upper,  lower,  tnum,  range,  spacing,  cyclevalue,  EfFtc  :  real; 
resultants ;  ARRAY[1..200]  OF  REAL; 

(* . *) 
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PROCEDURE  Getrange; 

VAR  j :  integer; 

BEGIN  (*  GetRange  *) 

FOR  J  :=  1  TO  200  DO  Resultams[J]  :=  0.0; 
TextColor(LightRed); 

GOTOXY(42, 16);  WriteC  *); 

GOTOXY(42,16); 

CASE  cyclevar  OF 

1  :  Write  ('Input  range  of  values  for  ’,CHR(229)); 

2  :  Write  (’Input  range  of  values  for  ’,CHR(225)); 

3  :  Write  (’Input  range  of  values  for  \CHR(21»; 

4  :  Write  (’Input  range  of  values  for  ’,CHR(230)); 

5  :  Write  (’Input  range  of  values  for  k’); 

END; 

Readln; 

numcycles  :=  1; 

GOTOXY(42, 16);  Write(’  ’); 

TextColor(LightRe<i), 

GOTOXY(42,16);  write(’Lower  =  ’);  GetNumber  (50, 16, Lower); 
GOTOXY(42, 16);  Writer  ’); 

TextColor(LightRed); 

GOTOXY(42,16);  write(’Upper  =  ’);  GetNumber  C-0,16,Upper); 
GOTOXY(42, 16);  Writer  ’); 

TextColor(LightRed); 

GOTOXY(42,16);  write(’Number  of  points  =  ’); 

GetNumber  (62,16,tnum); 

IF  ((tnum  >  0.0)  AND  (tnum  <  200.0))  THEN  BEGIN 
numcycles  :=  Round(tnum); 
spacing  :=  (Upper  -  Lower)/numcycles;  END 
ELSE  spacing  (Upper  -  Lower); 

GOTOXY(42, 16);  Writer  ’); 

TextColor(Y  ellow); 

GOTCXY^Z.lb);  wTite(’Beginning  the  calculations  ....  ’); 
TextColor(Whitej; 

END;  (*  GetRange  *) 

(* . *) 

PROCEDURE  Saveresults; 

VAR  Saveit :  Char; 

Savename.Tempstr  :  String60; 

DestFile  :  Text; 
j.multdelay :  integer; 

BEGIN  (*  Saveresults  *) 

TextColor(LightRed); 

C«OTOXY(42, 16),  Writer  ’); 

GOTOXY(42,16);  write(’Save  the  results  ?  [Y/NJ  ’); 

Saveit  :=  GetChariTO.l^True.True.CY’.’y’/N’/n’J.’Y’.’Y’); 

IF  Saveit  =  ’Y’  THEN  BEGIN 
GOTOXY(42, 16);  Writer  ’); 

GOTOXY(42,16);  writer  Filename;  ’); 

Readln(Savename); 

Assign(DestFile.Savename); 

Rewrite(DesiFiie); 
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CASE  Cyclevar  OF 

1  :  TempStr  :=  ’sigma  varied’; 

2  :  Tempstr  :=  ’beta  varied’; 

3  :  Tempstr  :=  ’kappa  varied’; 

4  :  Tempstr  :=  ’xi  varied’; 

5  :  Tempstr  :=  ’mu  varied’; 

END; 

Writeln(DestFilc,TempStr); 

Str(Lowen  10:5, Tempstr);  Writeln(DestFile, ’Lower  =  ’  +  TempStr); 
Str(Upper:  10:5, Tempstr);  Writeln(DestFile, ’Upper  =  ’  +  TempStr); 
Str(Spacing:10:5,Tempstr);  Writeln(DestFile,’Spacing  =  ’  +  TempStr); 
Str(numcycles, Tempstr);  Writeln(DestFile,’Num.Pt.=  ’  +  TempStr); 
Str(FwdConst:  10:6, Tempstr);  Writeln(DestFile,’ Const  =  ’  +  TempStr); 
FOR  j  :=  1  TO  numcycles  DO  BEGIN 
Str(Resultams(j]:  10:6, TempStr); 

Writeln(DestFile,TempStr); 

END; 

Close(DestFile); 

END; 

Te  x  tColor(Wh  i  te ) ; 

END;  (*  Saveresults  *) 

(* . *) 

BEGIN  (*  CycieThru  *) 

Lower  :=  0.0;  Upper  :=  1.0;  tnum  :=  1.0; 

FOR  cyclevar  :=  1  to  5  DO 
IF  DFV (cyclevar]  =  3.0  THEN  BEGIN 
TextColor(LightRed); 

GetRange; 

cyclevalue  :=  Lower, 

FOR  cyclecount  :=  1  TO  numcycles  DO  BEGIN 
CASE  cyclevar  OF 

1  :  sigma  :=  cyclevalue; 

2  :  beta  :=  cyclevalue; 

3  :  xi  :=  cyclevalue; 

4  :  mu  :=  cyclevalue; 

5  :  kappa  :=  cyclevalue; 

END; 

Romb(sigma,beta,xi,mu,kappa,-5,5, 1 E-04, 1  ,h); 
cyclevalue  :=  cyclevalue  +  spacing; 
h  :=  (2/Sqrt(pi))*h*(pi  *  piyxi)*exp(2*mu*kappa*xi); 

Effic  :=  FwdConst*h; 

Effic  :=  Effic/O+Effic);  (*  Accounts  for  Rsat  *) 
Resultants[cyc!e^:ount]  :=  Effic; 

GOTOXY(42,16);  Write(’  ’); 

GOTOXY(42,16);  write(’h  =  ’,n:10:6); 

GOTOXY(58,16);  write(’Eff.  =  ’.Effic:  10:6); 

END; 

Saveresults; 

END; 

END;  (*  CycieThru  *) 

(* . *) 
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BEGIN  (*  Calllnteg  *) 

Cycle  :=  False; 

InputPamis(sigma,beta,xi,mu,kappa,FwdConst,DFV); 

FOR  cnt  :=  1  TO  5  DO 
IF  DFV[cnt]  =  3.0  THEN  Cycle  :=  Tme; 

GOTOXY(5 1,24);  WriteC  ’); 

TextColor(Y  ellow); 

GOTOXY(51,24);  Write(’Continue  ?  [Y/N]  ’); 

ok  :=  GetChar(73, 24 .True, True, [’YYyYNYn’J.’YYY’); 

TextColor  (White); 

IF  ok  =  ’Y’  THEN 
IF  Cycle  THEN  CycleThru 
ELSE  BEGIN 

Romb(sigma,beta,xi,mu,kappa,-S,5,  IE-04, 1  ,h); 
h  ;=  (2/Sqrt(pi))*h'*(pi  *  pi/xi)*exp(2*mu*kappa*xi); 
GOTOXY(42, 16);  WriteC  ’); 

write(CHR(BEEP)); 

TextColor(LightRed); 

GOTOXY(42,16);  write(  h  =  ’,h:10:6); 

GOTOXY(58,16);  write(’Eff.  =  \FwdConst*h:10:6); 
TextColor(White); 

END; 

END;  (*  Calllnteg  *) 

(* . *) 

BEGIN  (*  MAIN  *) 

REPEAT 

ClrScr; 

Calllnteg; 

GOTOXY(5 1,24);  Writer  ’); 

TextColor(Yellow); 

GOTOXY(5 1 ,24);  Wnte(’Run  again  ?  [Y/N]  ’); 

Continue  :=  GETCHAR(73, 24, True, True, [’Y’.’y’/N’.’n’l.’Y’.’Y’) 
TextColor(Whitc); 

UNTIL  Continue  »  ’N’; 

END.  (*  MAIN  *) 

(* . . . *) 
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